Bài báo này trình bày kết quả nghiên cứu mô phỏng ảnh hưởng của việc bổ sung hydro trên đường ống nạp tới sự hình thành muội than trong động cơ diesel có trang bị hệ thống luân hồi khí thải.
HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 ẢNH HƯỞNG CỦA HYDRO TỚI HÌNH THÀNH MUỘI THAN CỦA ĐỘNG CƠ DIESEL CĨ LN HỒI KHÍ THẢI EFFECT OF HYDROGEN ON SOOT FORMATION OF DIESEL ENGINE INTERGRATED EXHAUST GAS RECIRCULATION TRỊNH XN PHONG1*, HỒNG ĐÌNH LONG1, PHẠM NGỌC ANH2 Viện Cơ khí động lực, Trường Đại học Bách khoa Hà Nội Khoa Công nghệ ô tô, Trường Cao đẳng Cơ khí nơng nghiệp *Email liên hệ: txphong@nute.edu.vn Tóm tắt Bài báo trình bày kết nghiên cứu mô ảnh hưởng việc bổ sung hydro đường ống nạp tới hình thành muội than động diesel có trang bị hệ thống luân hồi khí thải Động nghiên cứu động diesel R180, kỳ, xi lanh Phần mềm mô phần mềm AVL Boost Chế độ nghiên cứu chế độ 2600 (v/ph) Kết nghiên cứu cho thấy hydro có vai trị tích cực việc giảm hình thành muội than Cụ thể với tỷ lệ hydro tối ưu khối lượng muội than giảm trung bình 22% so với động nguyên Từ khóa: Động diesel, bổ sung hydro, ln hồi khí thải, giảm hình thành muội than Abstract This paper presents the results of simulation study on the effect of hydrogen addition to the intake manifold on soot formation of diesel engines equipped with exhaust gas recirculation system The research engine is a single-cylinder R180, 4stroke diesel engine The simulation software is AVL Boost software The study mode is the 2600rpm mode Research results show that hydrogen has play an important role in reducing soot formation Specifically, with the optimal hydrogen ratio, the amount of soot reduced on average by 22% compared to the original engine Keywords: Diesel engine, hydrogen addition, EGR, reduction soot formation Giới thiệu Động diesel nguồn động lực vô quan trọng ngành giao thông vận tải Tuy nhiên khí thải lại chứa nhiều chất độc hại muội than, NOx chất gây ô nhiễm khác CO, HC, vấn đề đáng lo ngại toàn cầu Nồng độ chất nhiễm phụ thuộc vào đặc tính động 132 thông số vận hành Để xử lý tất thành phần phát thải độc hại cần trang bị hệ thống xử lý khí thải gồm xử lý DOC (Diesel Oxidation Catalysts) để xử lý CO HC, xử lý DPF (Diesel Particulate Filter) để lọc bụi, xử lý SCR (Selective Catalytic Reduction) kết hợp DEF (Diesel Exhaust Fluid) để giảm NOx Việc trang bị xử lý khí thải động diesel giúp giảm đáng kể thành phần phát thải động Tuy nhiên, kết cấu hệ thống xử lý khí thải cồng kềnh, chi phí tốn sử dụng phức tạp Phương pháp luân hồi khí thải EGR (Exhaust Gas Recirculation) xem phương pháp hiệu để giảm phát thải NOx động diesel sử dụng từ lâu Tuy nhiên, công nghệ lại làm giảm hiệu suất động tăng hàm lượng CO, HC PM khí thải nên cần sử dụng kết hợp với biện pháp xử lý khác Nhiên liệu hydro coi nguồn nhiên liệu đầy tiềm với nhiệt trị khối lượng cao, tốc độ cháy nhanh hy vọng loại nhiên liệu khả thi cho việc giảm phát thải độc hại cho động diesel có ln hồi khí thải Theo Bose [2] cộng cung cấp hydro với khối lượng 0,15kg/h cho động có EGR0%, EGR10% EGR20% , hệ số nạp giảm thực EGR bổ sung hydro Khi thực EGR10% EGR20% phát thải khói, CO, CO HC tăng Nếu bổ sung hydro cho trường hợp luân hồi phát thải kể giảm so với động nguyên Giảm nhiều 42% (khói), 45,8% (CO 80% tải), 40,5% (CO2 80% tải ), 58% (HC 80% tải) với EGR0% Khi thực EGR10% 20% phát thải NOx giảm Việc bổ sung hydro làm gia tăng phát thải NO x thấp động nguyên Yadav [3] cộng nghiên cứu bổ sung hydro với khối lượng 40mg/hr động xi lanh, công suất định mức 3,5kW kết cho thấy BTE tăng 1,83%; CO2 giảm 10%; CO giảm 21,7%; HC giảm 5,13% SỐ ĐẶC BIỆT (10-2021) HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 Một nghiên cứu khác Việt Nam tác giả Anh Tuan Hoang [4] cộng nghiên cứu bổ sung hydro với lưu lượng 25 l/p với EGR10% EGR20% động D243 Kết cho thấy: Hiệu suất nhiệt động tăng lên 5%, tiêu hao nhiên liệu giảm đến 8% Lượng khí thải CO HC giảm mạnh khí thải NOx tăng so với động nguyên Có thể thấy việc cung cấp hydro đường ống nạp cho động diesel có ln hồi khí thải chia làm trường hợp Trường hợp 1: Thay phần diesel hydro cho đảm bảo tổng lượng đầu vào cung cấp không đổi Trường hợp 2: Giữ nguyên lượng nhiên liệu diesel cung cấp bổ sung hydro với lưu lượng không đổi Các kết khẳng định phát thải muội than động diesel giảm nguyên nhân: Một lượng diesel cung cấp đầu vào giảm, hai hydro cải thiện trình cháy giúp trình cháy trở nên hồn tồn từ góp phần giảm phát thải Tuy nhiên chưa có cơng trình nghiên cứu ảnh hưởng hydro tới tốc độ hình thành muội than khối lượng hình thành muội than động R180, loại máy phổ biến dùng nông nghiệp Việt Nam Chính báo trình bày kết nghiên cứu mơ ảnh hưởng việc bổ sung hydro đường ống nạp tới tốc độ hình thành muội than tổng khối lượng muội than hình thành động diesel có ln hồi khí thải Nghiên cứu mơ 2.1 Động nghiên cứu Động nghiên cứu động diesel R180, bốn kỳ, xi lanh Thông số kỹ thuật động nghiên cứu trình bày Bảng Bảng Thông số kỹ thuật động R180 [1] Thơng số kỹ thuật Mơ hình mơ trình bày Hình Các thơng số nhập vào mơ hình bao gồm đường kính xi lanh, hành trình piston, kiểu loại động cơ, Mơ hình cho q trình cháy động lưỡng nhiên liệu mơ hình cháy AVL MCC Mơ hình sử dụng cho nhiên liệu diesel đa nhiên liệu diesel loại nhiên liệu khác Mơ hình có khả dự đốn tốc độ giải phóng nhiệt động hình thành hỗn hợp khơng đồng nhất, ảnh hưởng q trình ln hồi dự đốn hình thành NOx, CO muội than [5] d mc u dQ dm dV dQF = - pc + - å w - hBB BB da da da da da (1) Trong đó: d mc u da - pc - Biến đổi nội bên xi lanh; dV da dQ F da - Cơng chu trình thực hiện; - Nhiệt lượng cấp vào; - Tổn thất nhiệt qua vách; a hBB dmBB da - Tổn thất enthalpy lọt khí; mc - Khối lượng môi chất bên xi lanh; u - Nội năng; pc - Áp suất bên xi lanh; V - Thể tích xi lanh; a - Góc quay trục khuỷu; hBB - Trị số enthalpy; - Biến thiên khối lượng dòng chảy Giá trị Mã động R180 Kiểu động Diesel, kỳ, xi lanh, Công suất định mức 5,15kW 2600v/ph Tiêu hao nhiên liệu 278 (g/kWh) Tỉ số nén 20±1 Dung tích cơng tác 0,402 (lít) SxD 80 x 80 (mm) Kiểu làm mát Nước, đối lưu tự nhiên Thời điểm phun 22 ± Độ (BTDC) Áp suất phun 140 bar SỐ ĐẶC BIỆT (10-2021) 2.2 Mơ hình mơ Hình Mơ hình mơ AVL Boost C1 động diesel R180; CO1 làm mát khí luân hồi; TH1 van EGR;I1 vòi phun H2; SB1, SB2: điều kiện biên; MP1đến MP6: điểm đo;J1, J2: điểm nối đường ống; E1: tử đặt thông số chung 133 HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 Quá trình truyền nhiệt từ buồng cháy qua thành buồng cháy nắp xi lanh, piston, lót xi lanh tính dựa vào phương trình truyền nhiệt sau: Qwi = Ai a w Tc - Twi (2) Trong đó: Qwi - Nhiệt lượng truyền cho thành (nắp xi lanh, pít-tơng, lót xi lanh); Ai - Diện tích truyền nhiệt (nắp xi lanh, pít tơng, lót xi lanh); aw - Hệ số truyền nhiệt; Tc - Nhiệt độ môi chất xi lanh; t char - Hệ số hòa trộn (-); n1, n2, n3 - Hằng số mơ hình (-) Theo quan điểm Hiroyasu, q trình hình thành phát thải muội than phụ thuộc trực tiếp vào nồng độ nhiên liệu Q trình xy hóa định lượng muội than thực tế xi lanh hàm lượng xy cịn dư thừa vùng phản ứng Hai phản ứng (5) (6) diễn đồng thời lượng muội than thực tế sinh phụ thuộc vào tốc độ ô xy hóa tốc độ hình thành muội than xi lanh: dms , f Twi - Nhiệt độ thành (nắp xi lanh, pít tơng, lót xi dt lanh) dt = As2 dms , f dt t char (msoot ) n ( PO2 , ref n1 = As , f E so PO2 ) n e Rm T sf dm fb, diff æ p Rm T ữ e ỗ dt ỗố pref ữứ E (3) (4) Trong đó: dms, f dt dm s ,o dt Es , f (5) Rm T dt - s ,o po dms ,o = Cso msoot p1.8 e Rm T dt p (6) Trong đó: dms, f dt dms ,o dt - Tốc độ hình thành muội than, (kg/s); - Tốc độ xy hóa muội than, (kg/s); po2 - Áp suất thực/áp suất tham chiếu, (Pa); Esf, Eso - Năng lượng kích hoạt q trình hình thành xy hóa, (J/mol); msoot - Khối lượng muội than, (kg); Cs,f - Hệ số hình thành muội than (-); Cs,o - Hệ số xy hóa muội than (-) 2.3 Trình tự mơ - Thực nghiệm động NB để xác định số liệu xây dựng mơ hình mơ phỏng; - Tốc độ hình thành muội than, (kg/s); - Tốc độ xy hóa muội than, (kg/s); dm fb ,diff - Tốc độ trình cháy khuếch tán, (kg/s); po2, po2,ref - Áp suất thực/áp suất tham chiếu, (Pa); Esf, Eso - Năng lượng kích hoạt q trình hình thành xy hóa, (J/mol); msoot - Khối lượng muội than, (kg); As,f - số hình thành muội than (-); As,o - Hệ số xy hóa muội than (-); 134 - E Cơ chế hình thành phát thải muội than mơ AVL-Boost sử dụng mơ hình Schubiger [6] coi q trình cháy gồm hai vùng mơ hình Hiroyasu [7, 8] coi trình cháy đa vùng Schubiger quan tâm đến hai phản ứng định đến trình hình thành phát thải muội than theo chế phản ứng động Phản ứng hình thành muội than phụ thuộc vào tốc độ cháy khuếch tán xi lanh Q trình xy hóa định lượng muội than thực tế xylanh hàm lượng ô xy cịn dư thừa vùng phản ứng Thơng số τchar liên quan tới tổng nhiệt tỏa theo hai phương trình dms2 = Cs , f m fv p 0.5 e - Mô đánh giá tính xác mơ hình động NB; - Mô động NB cung cấp nhiên liệu diesel với giá trị nhiên liệu 0,2; 0,4; 0,6; 0,8; 1,0 (kg/h) ký hiệu từ Gnl1 đến Gnl5 tương đương với gia tăng tải trọng động tốc độ 2600 (v/ph); - Mơ ln hồi khí thải với tỷ lệ luân hồi 10% 20% động NB; - Mơ cung cấp hydro có động NB có trang bị ln hồi khí thải với tỷ lệ hydro bổ sung 4%; 6% 8% khối lượng diesel cung cấp; - Đánh giá hiệu giảm phát thải muội than hydro cho động diesel có trang bị hệ thống luân hồi khí thải SỐ ĐẶC BIỆT (10-2021) HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 Gnl khác Lượng hydro xác định theo công thức (8): Kết nghiên cứu 3.1 Đánh giá độ tin cậy mơ hình Trong đó: 700 mH - Khối lượng hydro bổ sung; 600 500 400 300 200 100 Ne_tn k - Tỷ lệ thay 3.3 Ảnh hưởng hydro bổ sung tới hệ số lambda 4.5 Tốc độ động (v/ph) Hình Đồ thị đánh giá độ tin cậy mơ hình Để có thơng số đầu vào nhập mơ hình phần mềm ALV Boost cần phải nghiên cứu thực nghiệm Hình trình bày diễn biến công suất suất tiêu hao nhiên liệu đặc tính ngồi động diesel thực nghiệm mô Kết mô thực nghiệm sai lệch khơng q 5% kết mơ từ mơ hình sử dụng làm số liệu tham chiếu nghiên cứu ảnh hưởng việc bổ sung hydro đến hình thành muội than động R180 với tỷ lệ luân hồi khác tốc độ 2600 (v/ph) Trong nghiên cứu, tỷ lệ EGR tính cơng thức (7): megr ´ 100 (7) mi Trong đó: megr ,mi lượng khơng khí ln hồi, lượng khơng khí nạp lấy từ AVL Boost 3.2 Lượng nhiên liệu hydro bổ sung 0.6 Hydro bổ sung (g/p) 0.4 0.3 ố ệ ữ ỷ ệ lưu lượ 8% ổ Hình trình bày mối quan hệ tỷ lệ hydro bổ sung lưu lượng hydro bổ sung cho trường hợp SỐ ĐẶC BIỆT (10-2021) E20 E10H4 E10H6 E10H8 2.5 D20H4 E20H6 3.5 1.5 Gnl1 Gnl2 Gnl3 Gnl4 Gnl5 Gnl (kg/h) Hình Ảnh hưởng tỷ lệ hydro bổ sung tới hệ số lambda Hình trình bày ảnh hưởng việc bổ sung hydro động R180 với tỷ lệ khí thải luân hồi khác tốc độ 2600 (v/ph) tới hệ số lambda Kết mô cho thấy việc bổ sung hydro làm cho lamba giảm Điều lý giải hydro chiếm chỗ khơng khí buồng đốt Trong nghiên cứu mô lấy giá trị lambda lớn 1,2 để đảm bảo đủ khơng khí đốt cháy hoàn toàn diesel hydro Theo quan điểm nên thực đồng thời ln hồi bổ sung hydro trường hợp Gnl1, Gnl2 Gnl3 Có thể thấy thực bổ sung hydro động có trang bị hệ thống luân hồi khí thải thời điểm bắt đầu hình thành muội than ln có xu hướng muộn so với trường hợp không bổ sung hydro 0.1 6% Tỷ lệ hydro bổ sung E10 Hình đến Hình 10 trình bày kết mơ tốc độ hình thành muội than trường hợp Gnl1 đến Gnl3 0.2 4% NB 3.4 Ảnh hưởng hydro tới tốc độ hình thành muội than Gnl1 Gnl2 Gnl3 Gnl4 Gnl5 0.5 Hệ số Lambda 1300 1500 1800 2000 2200 2400 2600 EGR (%) = (8) Gnl - Lượng nhiên liệu diesel cung cấp; ge(g/kWh) Ne (kW) mH = Gnl ´ k 800 Tại trường hợp ln có tỷ lệ hydro cung cấp tối ưu Tại tỷ lệ tối ưu ln hình thành vùng diện tích tạo trục hoành đường cong tốc độ hình thành muội than nhỏ so với động NB Các tỷ lệ tối ưu liệt kê Bảng 135 HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 Tốc độ hình thành soot (g/s) Gnl1, E10 0.0008 0.0007 NB 0.0006 E10 0.0005 E20 E20H4% 0.0006 0.0004 E10H6% 0.0004 E10H8% 0.0002 E20H6% E20H8% 0.0002 0.0001 0 Góc quay trục khuỷu (độ) Góc quay trục khuỷu (độ) Hình Ảnh hưởng hydro tới tốc độ hình thành Hình Ảnh hưởng hydro tới tốc độ hình thành muội than Gn1, E10 muội than Gnl1, E20 0.016 Gnl 2, E10 0.012 NB 0.01 E10 0.008 E10H4% 0.006 E10H6% 0.004 E10H8% 0.002 Tốc độ hình thành soot (g/s) Tốc độ hình thành soot (g/s) 0.014 Gnl2, E20 0.014 0.012 0.01 NB E20 0.008 E20H4% 0.006 E20H6% 0.004 E20H8% 0.002 0 Góc quay trục khuỷu (độ) Góc quay trục khuỷu (độ) Hình Ảnh hưởng hydro tới tốc độ hình thành Hình Ảnh hưởng hydro tới tốc độ hình thành muội than Gnl2, E10 muội than Gnl2, E20 0.08 Gnl3, E20 0.08 Tốc độ hình thành soot (g/s) Tốc độ hình thành soot (g/s) 0.09 NB 0.07 0.06 E20 0.05 0.04 E20H4% 0.03 E20H6% 0.02 E20H8% 0.01 Gnl3, E10 0.07 0.06 0.05 NB E10 0.04 E10H4% 0.03 E10H6% 0.02 E10H8% 0.01 -0.01 Góc quay trục khuỷu (độ) Góc quay trục khuỷu (độ) Hình Ảnh hưởng hydro tới tốc độ hình thành Hình 10 Ảnh hưởng hydro tới tốc độ hình muội than Gnl3, E10 thành muội than Gnl3, E20 Bảng Tỷ lệ hydro tối ưu cho động R180 có luân hồi khí thải 136 Gnl1, E20 NB 0.001 0.0008 E10H4% 0.0003 Tốc độ hình thành soot (g/s) 0.0012 0.0009 Gnl1 Gnl1 Gnl2 Gnl2 Gnl1 Gnl2 E10 E20 E10 E20 E10 E20 H8% H6% H6% H6% H6% H4% 3.5 Ảnh hưởng hydro bổ sung tới khối lượng muội than hình thành Hình 11 trình bày kết mơ ảnh hưởng việc bổ sung hydro tới khối lượng muội than hình thành động R180 có trang bị luân hồi khí thải với tỷ lệ hydro cung cấp tối ưu theo Bảng Kết mô cho thấy: - Khi thực luân hồi khí thải khối lượng muội than có xu hướng tăng theo gia tăng tỷ lệ EGR; SỐ ĐẶC BIỆT (10-2021) HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 - Khi thực bổ sung hydro cho động R180 có ln hồi phát thải muội than giảm mạnh Cụ thể Gnl1 giảm 31,5%; Gnl2 giảm 25,4%; Gnl3 giảm 22,5% so với động NB; Soot (g/kWh) NB E10 E20 EH Gnl1 Gnl2 Gnl3 Gnl4 Gnl5 Lượng diesel cung cấp (kg/h) Hình 11 Ảnh hưởng hydro tới khối lượng muội than hình thành động R180 Việc phát thải muội than động diesel giảm mạnh cung cấp hydro giải thích theo bước sau: - Q trình thành muội than động diesel: nhiên liệu diesel cháy buồng đốt, phần nhiên liệu vùng có nhiệt độ cao tạo thành muội than Cơ chế hình thành muội than bao gồm trình: Nhiệt phân, tạo mầm, phát triển hạt, kết tụ, thiêu kết oxi hóa; - Q trình thành muội than thực luân hồi khí thải: Khi thực luân hồi khí thải, phần khí thải quay lại buồng đốt điều dẫn tới thiếu hụt xy khí ln hồi chiếm chỗ yếu tố dẫn đến gia tăng hình thành muội than đặc biệt chế độ tải trọng lớn; - Giảm muội than bổ sung hydro: nói trên, hydro với khả khuếch tán tốt bổ sung đường ống nạp nạp vào xi lanh giúp cho hỗn hợp khơng khí-hydro đồng Hỗn hợp điều kiện để giúp nhiên liệu diesel cháy cách triệt để tránh tạo vùng cục khó cháy, khơng cịn nhiên liệu thừa hình thành muội than giảm Thêm hydro có tốc độ cháy nhanh giúp lan tràn màng lửa để đốt cháy nhiên liệu nơi xa nguồn lửa trước xu-páp thải mở Tất yếu tố cải thiện trình cháy động cơ, nhờ phát thải muội than giảm đáng kể Kết luận Kết nghiên cứu mô ảnh hưởng việc bổ sung hydro đường ống nạp tới tốc độ khối lượng muội than hình thành động diesel R180 SỐ ĐẶC BIỆT (10-2021) - Hệ số lambda động giảm hydro tích chiếm chỗ lớn làm giảm lượng khí nạp cung cấp cho động cơ; - Tốc độ hình thành muội than khối lượng muội than hình thành giảm mạnh Cụ thể khối lượng muội than hình thành giảm trung bình 22% tỷ lệ hyrdo bổ sung tối ưu vùng lưu lượng nhiên liệu diesel cung cấp từ 0,2kg/h đến 0,6kg/h; có trang bị hệ thống luân hồi khí thải chế độ 2600 (v/ph) khẳng định: - Ở vùng lưu lượng nhiên liệu lớn 0,8kg/h không nên thực bổ sung hydro lúc lambda thấp; Như bổ sung hydro đường ống nạp phương án khả thi để khắc phục tăng hình thành muội than động diesel có ln hồi khí thải TÀI LIỆU THAM KHẢO [1] R180 Diesel engine operation manual, ed, 2018 [2] P K Bose and D Maji, An experimental investigation on engine performance and emissions of a single cylinder diesel engine using hydrogen as inducted fuel and diesel as injected fuel with exhaust gas recirculation, International Journal of Hydrogen Energy, Vol.34, pp.4847-4854, 2009 [3] V SinghYadav, S L Soni, and D Sharma, Performance and emission studies of direct injection C.I engine in duel fuel mode (hydrogendiesel) with EGR, International Journal of Hydrogen Energy, Vol.37, pp.3807-3817, 2012 [4] A T Hoang and V V Pham, A study on a solution to reduce emissions by using hydrogen as an alternative fuel for a diesel engine integrated exhaust gas recirculation, AIP Conference Proceedings, Vol.2235, p.020035, 2020 [5] A Boost, AVL Boost Combustion Model Austria, 2011 [6] R A Schubiger, K Boulouchous, and M K Eberle, Rußbildung und Oxidation bei der dieselmotorischen Verbrennung, MTZ Motortechnische Zeitschrift, Vol.63, pp.342-353, 2002/05/01, 2002 [7] H Hiroyasu, T Kadota, and M Arai, Development and Use of a Spray Combustion Modeling to Predict Diesel Engine Efficiency and Pollutant Emissions: Part Combustion Modeling, Bulletin of JSME, Vol.26, pp.569-575, 1983 137 HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 [8] H a T K Hiroyasu, Models for Combustionand Formation of Nitric Oxide and Soot in Direct Injection Diesel Engines., SAE paper 760129, 1976 138 Ngày nhận bài: Ngày nhận sửa: Ngày duyệt đăng: 28/6/2021 06/8/2021 28/8/2021 SỐ ĐẶC BIỆT (10-2021) ... (độ) Hình Ảnh hưởng hydro tới tốc độ hình thành Hình 10 Ảnh hưởng hydro tới tốc độ hình muội than Gnl3, E10 thành muội than Gnl3, E20 Bảng Tỷ lệ hydro tối ưu cho động R180 có ln hồi khí thải. .. Hình Ảnh hưởng hydro tới tốc độ hình thành Hình Ảnh hưởng hydro tới tốc độ hình thành muội than Gnl2, E10 muội than Gnl2, E20 0.08 Gnl3, E20 0.08 Tốc độ hình thành soot (g/s) Tốc độ hình thành. .. bổ sung hydro đường ống nạp tới tốc độ hình thành muội than tổng khối lượng muội than hình thành động diesel có ln hồi khí thải Nghiên cứu mô 2.1 Động nghiên cứu Động nghiên cứu động diesel R180,