Nghiên cứu ảnh hưởng quá trình phun nhiều giai đoạn khí thải động cơ Vikyno RV 125-2 bằng phương pháp mô phỏng

ĐẠI HỌC QUOÁC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA HỒ THANH BÌNH NGHIÊN CỨU ẢNH HƯỞNG Q TRÌNH PHUN NHIỀU GIAI ĐOẠN ĐẾN KHÍ THẢI ĐỘNG CƠ VIKYNO RV125-2 BẰNG PHƯƠNG PHÁP MÔ PHỎNG Chuyên Ngành: Kỹ Thuật Ô Tô - Máy Kéo Mã số: 605235 LUẬN VĂN THẠC SĨ TP.HỒ CHÍ MINH, tháng 06 năm 2014 CƠNG TRÌNH ĐƯỢC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA-ĐẠI HỌC QUỐC GIA TP.HCM Cán hướng dẫn khoa học: TS.NGUYỄN LÊ DUY KHẢI (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán chấm nhận xét 1: TS.NGUYỄN NGỌC DŨNG (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán nhận xét 2: PGS.TS.ĐẶNG THÀNH TRUNG (Ghi rõ họ, tên, học hàm, học vị chữ ký) Luận văn thạc sĩ ñược bảo vệ Trường Đại Học Bách Khoa Tp.HCM, ĐHQG Tp.HCM, Ngày 28 Tháng 07 Năm 2014 Thành phần Hội Đồng ñánh giá luận văn thạc sĩ gồm: (Ghi rõ họ, tên, học hàm, học vị Hội ñồng chấm bảo vệ luận văn thạc sĩ) TS.HUỲNH THANH CÔNG TS.PHẠM TUẤN ANH TS.NGUYỄN NGỌC DŨNG PGS.TS.ĐẶNG THÀNH TRUNG TS.NGUYỄN LÊ DUY KHẢI Xác nhận Chủ tịch Hội ñồng ñánh giá LV Trưởng Khoa quản lý chuyên ngành sau luận văn ñã ñược chỉnh sữa chữa (nếu có) CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA KỸ THUẬT GIAO THÔNG Độc Lập – Tự Do – Hạnh Phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: HỒ THANH BÌNH MSHV: 10130459 Ngày, tháng, năm sinh: 20/03/1985 Nơi sinh: Long An Chuyên ngành: Kỹ thuật ô tô – máy kéo Mã số: 605235 1.TÊN ĐỀ TÀI: NGHIÊN CỨU ẢNH HƯỞNG QUÁ TRÌNH PHUN NHIỀU GIAI ĐOẠN ĐẾN KHÍ THẢI CỦA ĐỘNG CƠ VIKYNO RV125-2 BẰNG PHƯƠNG PHÁP MÔ PHỎNG NHIỆM VỤ CỦA LUẬN VĂN: 2.1 Tổng quan tình hình nghiên cứu tác động q trình phun nhiên liệu ñến khí thải ñộng Diesel phun trực tiếp Việt Nam Thế Giới 2.2 Nghiên cứu ứng dụng phần mềm mơ KIVA3V vào động cụ thể ( động Diesel Vikyno RV125-2 phịng thí nghiệm ñộng ñốt trong, Trường Đại Học Bách Khoa Tp.HCM) 2.3 Nghiên cứu ảnh hưởng “ lượng phun ban ñầu” ñến khí thải ñộng Diesel phun trực tiếp 2.4 Nghiên cứu ảnh hưởng “ thời gian nghỉ hai lần phun” đến khí thải động Diesel phun trực tiếp 2.5 Tổng kết ñánh giá NGÀY GIAO NHIỆM VỤ: 03/09/2013 NGÀY HOÀN THÀNH NHIỆM VỤ: 20/06/2014 HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN (Ghi ñầy ñủ học hàm, học vị, họ, tên): TS NGUYỄN LÊ DUY KHẢI Tp.HCM, ngày tháng năm 2014 CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN ĐÀO TẠO (Họ tên chữ ký) (Họ tên chữ ký) TRƯỞNG KHOA (Họ tên chữ ký) LỜI CẢM ƠN Qua thời gian thực luận văn tốt nghiệp chuyên ngành kỹ thuật tơ – máy kéo, em nhận giúp đỡ tận tình q thầy bạn bè nên hồn thành thời gian, ñạt yêu cầu Em xin chân thành cảm ơn! Toàn thể thầy trường Đại Học Bách Khoa Tp.HCM tận tình dạy dỗ chúng em thời gian học tập trường Em xin chân thành cảm ơn thầy cô khoa Kỹ Thuật Giao Thông thầy cơ, cộng tác viên phịng thí nghiệm Đã giúp đỡ em suốt q trình thực nghiệm Đặc biệt em xin chân thành cảm ơn Thầy T.s Nguyễn Lê Duy Khải giáo viên hướng dẫn em ñã tạo ñiều kiện thuận lợi, hướng dẫn giúp ñỡ tận tình em trình thực luận văn Xin chân thành bạn học viên lớp ñã tham gia giúp ñỡ thời gian thi công mơ hình thực nghiệm Trong q trình thực luận văn cịn nhiều khiếm khuyết, thiếu sót kính mong thầy cô thông cảm lượng thứ! Cuối em xin kính chúc q thầy Khoa Kỹ Thuật Giao Thông Trường Đại Học Bách Khoa Tp.HCM sức khỏe, hạnh phúc ! Học viên thực Hồ Thanh Bình -i- LỜI CAM ĐOAN Tơi xin cam đoan ñề tài nghiên cứu thực Các số liệu kết luận nghiên cứu trình bày luận văn chưa cơng bố nghiên cứu khác Tơi xin chịu trách nhiệm nghiên cứu Học viên Hồ Thanh Bình TĨM TẮT Một vấn ñề quan trọng giảm chất phát thải gây nhiễm từ động Diesel bồ hóng NOx Phun nhiều giai ñoạn phương pháp hiệu để giảm khí thải động Diesel Luận văn nghiên cứu ảnh hưởng phun nhiên liệu hai lần chu trình để điều tra q trình hình thành khí thải động Diesel phun trực tiếp Vikyno RV125-2 Nghiên cứu thông số phun tách, bao gồm lượng phun ban ñầu thời gian nghỉ hai lần phun ñược thực hiện, phần mềm mô (KIVA3V) Trong luận văn, nghiên cứu với lượng phun ban ñầu 90, 80, 75, 35% thời gian nghỉ hai lần phun ñược thay ñổi ñể ñánh giá ảnh hưởng phun tách đến khí thải bồ hóng NOx Kết phun tách với lượng phun ban ñầu 80%, lượng phun sau 20% thời gian nghỉ hai lần phun 4ºCA NOx giảm 62,27%, bồ hóng 53,15% tăng so với phun đơn ABSTRACT One of the important problems in reducing of pollutant emission from DI diesel engines is trade-off between soot and NOx Multiple injections is one of the most powerful method that makes the chance to shift the trade-off curve closer to origin This thesis introduces the research into influence of two times the fuel injection within a cycle to investigate emissions formation processes in a direct injection Vikyno RV125-2 diesel engine The study of split injection parameters, including the delay dwell and the fuel quantity injected between injection pulses are carried out, using simulation (KIVA-3V) In this thesis, the five different split injection schemes, in which 90, 80, 75, 50 and 35% of total fuel injected in the first pulse, have been considered And the delay dwell between injections was aslo changed to appreciate the effect of split injection on the soot and NOx emissions The results show that the split injection of scheme of 80.4.20, in which 80% of total fuel injected in the first pulse is accompanied with the 4ºCA delay dwell between injections, reduces 62,27%, NOx, increase 53,15% soot emissions in comparison to the single injection strategy -ii- MỤC LỤC Trang Lời cảm ơn -i- Tóm tắt -ii- Mục lục -iiiDanh mục hình ảnh -vi- Danh mục bảng biểu -viiiDanh mục chữ viết tắt -ix- Danh mục phụ lục -xi- CHƯƠNG I: TỔNG QUAN 1.1 Lý chọn ñề tài .1 1.2 Nghiên cứu tổng quan nước 1.3 Mục tiêu ñề tài ñối tượng nghiên cứu 1.3 Nội dung phạm vi nghiên cứu 1.4 Phương pháp nghiên cứu 1.5 Nghiên cứu tổng quan .5 1.6 Ý nghĩa ñề tài CHƯƠNG II: CƠ SỞ LÝ THUYẾT 2.1 Khái qt chất gây nhiễm khí thải ñộng Diesel ·················7 2.1.1 Oxit carbon (CO) Carbonic (CO2) ·····················································7 2.1.2 Oxit nitoric (NOx) ··················································································7 2.1.3 Hydrocarbon (HC) 2.1.4 Bồ hóng .8 2.1.5 Oxit lưu huỳnh (SO2) 2.2 Cơ chế hình thành chất ô nhiễm 2.2.1 Cơ chế hình thành bồ hóng 2.2.2 Cơ chế hình thành NOx 2.3 Các biện pháp giảm ô nhiễm ñộng Diesel .12 -iii- 2.4 Những vấn đề q trình phun nhiên liệu 12 2.4.1.Yêu cầu ñối với hệ thống nhiên liệu diesel .12 2.4.2 Phun nhiều giai ñoạn 13 2.4.3 Vấn ñề phun nhiên liệu 14 2.4.4 Cấu trúc tổng quát tia nhiên liệu 15 2.4.5 Góc mở tia phun 18 2.4.6 Độ xuyên thâu tia nhiên liệu buồng cháy 20 2.4.7 Sự phân bố kích thước hạt .23 2.5 Các mơ hình tốn sử dụng mơ 25 2.5.1 Mơ hình phân rã tia phun “Kelvin-Helmholtz Rayleigh-Taylor” .25 2.5.2 Mơ hình bay hạt nhiên liệu 27 2.5.3 Mơ hình dịng chảy rối “Renormalized Group k-epsilon” (RNG k - ε) 29 2.5.4 Mơ hình hình thành NOx 30 2.5.5 Mơ hình hình thành bồ hóng 31 CHƯƠNG III: TẠO LƯỚI VÀ THIẾT LẬP CÁC THÔNG SỐ ĐỘNG CƠ NGHIÊN CỨU 3.1 Cấu trúc chương trình mơ KIVA3V-ERC .35 3.2 Sơ lược động mơ .36 3.3 Xây dựng mơ hình động Vikyno RV125-2 .39 3.4 Sơ ñồ tiến hành thí nghiệm động RV125-2 hiệu chỉnh đồng thông số thực nghiệm mô 41 3.4.1 Giới thiệu sơ ñồ thiết bị 41 3.4.2 Điều kiện bước thực thí nghiệm .42 3.4.3 Các hiệu chỉnh ñể ñồng thông số mô thực nghiệm 42 3.4.3 Đồng áp suất xylanh không phun nhiên liệu .42 3.4.3.2 Đồng áp suất xylanh phun nhiên liệu 44 CHƯƠNG IV: ẢNH HƯỞNG LƯỢNG PHUN BAN ĐẦU ĐẾN KHÍ THẢI 4.1 Điều kiện mơ thơng số thay đổi 4.2 Kết mơ bình luận 47 51 -iv- 4.2.1 Ảnh hưởng lượng phun ban ñầu ñến áp suất nhiệt ñộ xylanh 51 4.2.2 Ảnh hưởng lượng phun ban đầu đến hình thành NOx bồ hóng 4.2.2.1 Ảnh hưởng lượng phun ban đầu ñến hình thành NOx 54 54 4.2.2.2 Ảnh hưởng lượng phun ban đầu đến hình thành bồ hóng 57 4.3 Kết luận 63 CHƯƠNG V: ẢNH HƯỞNG THỜI GIAN NGHỈ GIỮA HAI LẦN PHUN ĐẾN KHÍ THẢI 5.1 Điều kiện mơ thơng số thay đổi 64 5.2 Kết mơ bình luận 65 5.2.1 Ảnh hưởng thời gian nghỉ ñến áp suất nhiệt ñộ xylanh .65 5.2.2 Ảnh hưởng lượng phun ban đầu đến hình thành NOx 68 5.2.3 Ảnh hưởng lượng phun ban đầu đến hình thành bồ hóng .71 5.3 Kết luận .75 CHƯƠNG VI: KẾT LUẬN 6.1 Kết luận 77 6.2 Đề xuất hướng phát triển 77 Phụ lục -xiiTài liệu tham khảo -xlv- -v- DANH MỤC BẢNG BIỂU Số hiệu bảng biểu Trang 2-1 Hằng số mơ hình KH-RT 27 2-2 Các tốc độ phản ứng sử dụng mơ hình bồ hóng “8 bước” Foster 33 3-1 Thơng số kỹ thuật động Vikyno RV 125-2 37 4-1 Thơng số q trình mơ 49 4-2 Áp suất nhiệt độ trung bình cực đại trường hợp mô 53 4-3 Khối lượng NOx phát thải thay ñổi lượng phun ban ñầu 57 4-4 Khối lượng bồ hóng phát thải ứng với lượng phun ban ñầu khác 58 4-5 Giá trị bồ hóng, NOx, cơng suất theo lượng phun ban ñầu 61 4-6 Tỉ lệ tăng/giảm NOx, Bồ hóng, cơng suất phun đơn so với phun tách thay ñổi lượng phun ñầu khác 63 5-1 Nhiệt ñộ áp suất xylanh ứng thay ñổi thời gian nghỉ 66 5-2 Khối lượng NOx ứng với thời gian nghỉ khác 68 5-3 Khối lượng NOx, Bồ hóng, cơng suất ứng với thời gian nghỉ khác 73 5-4 Tỉ lệ tăng/giảm NOx, Soot, cơng suất phun đơn so với phun tách thay ñổi thời gian nghỉ khác 75 -viii- numinj numvel t1inj 8.12500E-03 tdinj 5.62500E-04 ca1inj -1.80000E+01 " Start Of First Injection Timing" cadinj 8.10000E+00 " Duration Of First Injection" tspmas 2.46600E-02 " 90%Spray Mass To Be Injected, In Grams" tnparc 4.00000E+03 t1inj 8.96528E-03 tdinj 6.25000E-05 ca1inj -5.90000E+00 " Start Of Second Injection Timing" cadinj 9.00000E-01 " Duration Of Second Injection" tspmas 2.74000E-03 " 10% Spray Mass To Be Injected, In Grams" tnparc 4.00000E+03 pulse 2.00000 injdist kolide tpi 3.23000E+02 turb 1.00000E+00 breakup 1.00000E+00 evapp 1.00000E+00 drnoz 1.20000E-02 dznoz 1.10700E+01 dthnoz 4.50000E+01 tiltxy 4.50000E+01 tiltxz 7.00000E+01 cone 1.00000E+00 dcone 1.00000E+00 anoz 2.26800E-04 smr 1.65000E-02 amp0 0.00000E+00 nsp 18 species data c14h30 mw 1.98394E+02 htf -2.55308E+12 o2 mw 3.20000E+01 htf 0.00000E+00 n2 mw 2.80160E+01 htf 0.00000E+00 co2 mw 4.40110E+01 htf -3.93150E+12 h2o mw 1.80160E+01 htf -2.38919E+12 h mw 1.00800E+00 htf 2.16024E+12 h2 mw 2.01600E+00 htf 0.00000E+00 o mw 1.60000E+01 htf 2.46810E+12 n mw 1.40080E+01 htf 4.70784E+12 oh mw 1.70080E+01 htf 3.88652E+11 co mw 2.80110E+01 htf -1.13805E+12 no mw 3.00080E+01 htf 8.97719E+11 c2h2 mw 2.60000E+01 htf 2.26731E+12 rn mw 6.00000E+02 htf 7.16677E+12 pn mw 1.20000E+01 htf 7.16677E+12 r mw 1.15000E+02 htf 0.00000E+00 b mw 2.30000E+02 htf 0.00000E+00 q mw 2.30000E+02 htf 0.00000E+00 stoifuel 2.00000E+00 stoio2 4.30000E+01 nreg presi 990000.000000000 -xxx- tempi 313.000000000000 tkei 0.100000000000000 scli 0.000000000000000E+000 er 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.232900000000000 mfrac 0.767100000000000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 mfrac 0.000000000000000E+000 for region the following will be used: equivalence ratio= 0.0000000000E+00 average molecular weight of the non-fuel species= 2.8852611774E+01 total density= 1.0976156694E-03 total moles per unit volume (= sum (rho/mw))= 3.8042159859E-05 pressure= 9.9000000000E+05 temperature= 3.1300000000E+02 species densities and mass fractions for this region are: 0.0000000000E+00 0.0000000000E+00 2.5563468941E-04 2.3290000000E-01 8.4198098002E-04 7.6710000000E-01 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 10 0.0000000000E+00 0.0000000000E+00 11 0.0000000000E+00 0.0000000000E+00 12 0.0000000000E+00 0.0000000000E+00 13 0.0000000000E+00 0.0000000000E+00 14 0.0000000000E+00 0.0000000000E+00 15 0.0000000000E+00 0.0000000000E+00 16 0.0000000000E+00 0.0000000000E+00 17 0.0000000000E+00 0.0000000000E+00 18 0.0000000000E+00 0.0000000000E+00 v( 1) = 46795 nrk nvalves isoot tcuts 1.20000E+03 drni 2.24800E-08 -xxxi- dnuci 1.25000E-07 drncr 1.25000E-07 2.00000E-08 surf 2.00000E+04 rhorn 2.00000E+00 rhopn 2.00000E+00 xnc 1.20000E+01 ee1 1.16000E+02 ee2 3.00000E+01 ee3 8.00000E+00 ee4 0.00000E+00 eka 3.00000E+01 ekb 1.52000E+01 ekt 9.70000E+01 ekz 4.10000E+00 fka 2.00000E+01 fkb 4.46000E-03 fkt 1.51000E+05 fkz 2.13000E+01 -********************************* -grams of each species in entire mesh: c14h30 2.385926894693E-03 o2 9.962417567406E-02 n2 5.036821284488E-01 co2 3.575381572020E-02 h2o 1.700187960183E-02 h 0.000000000000E+00 h2 3.384651788951E-04 o 0.000000000000E+00 n 0.000000000000E+00 oh 0.000000000000E+00 co 6.182466027474E-03 no 4.618738106000E-04 c2h2 3.438667789936E-04 rn 1.404176617561E-10 pn 4.837587505968E-05 r 1.088150820857E-05 b 8.443514956197E-07 q 1.184414274713E-02 particle mass distribution: entrained= 0.00000E+00, and on surfaces numbered -1 to nmoving: 0.000000000000000E+000 6.996224910056666E-006 parts/million by volume of each species in region 1: c14h30ppm 5.154250297838E+02 o2ppm 1.334294415457E+05 n2ppm 7.705260369115E+05 co2ppm 3.481753866983E+04 h2oppm 4.044598989197E+04 hppm 0.000000000000E+00 h2ppm 7.195489989474E+03 oppm 0.000000000000E+00 nppm 0.000000000000E+00 ohppm 0.000000000000E+00 coppm 9.459540279764E+03 noppm 6.596646198654E+02 c2h2ppm 5.668312409181E+02 rnppm 1.003014917588E-05 pnppm 1.727764290097E+02 rppm 4.055348292361E+00 bppm 1.573375367772E-01 qppm 2.207052696346E+03 crank= 130.12 fluxin= 0.00000E+00 fluxout= 0.00000E+00 wheat= 1.08889E+09 swirl ratio= -8.60588E-01 crenb= 130.12 0.00000E+00 2.65833E+09 1.08889E+09 -3.14673E+09 0.00000E+00 cremi= 130.12 3.48175E+04 9.45954E+03 6.59665E+02 Indicated Gross Work= 307.722831285856 Nm Indicated power= 6.15445662571711 kW ISFC= 80.1370502700607 g/kWh Execution time= 11715.5 s (= 195.26 m, = 3.25 h) Total time from beginning= 11713.2 s (= 195.22 m, = 3.25 h) -xxxii- Phụ Lục 6: Otape12.txt_Phun Tách 80%4CA20% -********************************* -grams of each species in entire mesh: c14h30 1.727884762759E-03 o2 9.924541946623E-02 n2 5.036818048265E-01 co2 3.562497567273E-02 h2o 1.683499661918E-02 h 0.000000000000E+00 h2 3.999856184527E-04 o 0.000000000000E+00 n 0.000000000000E+00 oh 0.000000000000E+00 co 6.841948592626E-03 no 4.414156920862E-04 c2h2 3.263930702825E-04 rn 7.481314905788E-11 pn 3.697506724328E-05 r 1.108016723516E-05 b 6.034725107756E-07 q 1.217845731920E-02 particle mass distribution: entrained= 0.00000E+00, and on surfaces numbered -1 to nmoving: 0.000000000000000E+000 6.696464304013395E-006 parts/million by volume of each species in region 1: c14h30ppm 3.728568849753E+02 o2ppm 1.327750333372E+05 n2ppm 7.696726593088E+05 co2ppm 3.465367242128E+04 h2oppm 4.000466022668E+04 hppm 0.000000000000E+00 h2ppm 8.493951454131E+03 oppm 0.000000000000E+00 nppm 0.000000000000E+00 ohppm 0.000000000000E+00 coppm 1.045700026560E+04 noppm 6.297477763910E+02 c2h2ppm 5.374319822003E+02 rnppm 5.338049675853E-06 pnppm 1.319118017342E+02 rppm 4.124814295201E+00 bppm 1.123273677363E-01 qppm 2.266837394128E+03 crank= 130.03 fluxin= 0.00000E+00 fluxout= 0.00000E+00 wheat= 1.11492E+09 swirl ratio= -8.66779E-01 crenb= 130.03 0.00000E+00 2.66578E+09 1.11492E+09 -3.12120E+09 0.00000E+00 cremi= 130.03 3.46537E+04 1.04570E+04 6.29748E+02 Indicated Gross Work= 305.269860609765 Nm Indicated power= 6.10539721219530 kW ISFC= 80.7809848988779 g/kWh Execution time= 9936.3 s (= 165.60 m, = 2.76 h) Total time from beginning= 9934.0 s (= 165.57 m, = 2.76 h) -xxxiii- Phụ Lục 7: Otape12.txt_Phun Tách 75%4CA25% -********************************* -grams of each species in entire mesh: c14h30 1.779179610354E-03 o2 9.963980795966E-02 n2 5.036817670821E-01 co2 3.501840293143E-02 h2o 1.658588075590E-02 h 0.000000000000E+00 h2 4.413566995672E-04 o 0.000000000000E+00 n 0.000000000000E+00 oh 0.000000000000E+00 co 7.321536267037E-03 no 4.272964558454E-04 c2h2 4.237243624169E-04 rn 2.071323284390E-10 pn 5.698554183557E-05 r 1.083598835796E-05 b 6.212571482979E-07 q 1.217042640845E-02 particle mass distribution: entrained= 0.00000E+00, and on surfaces numbered -1 to nmoving: 0.000000000000000E+000 1.757102910726001E-005 parts/million by volume of each species in region 1: c14h30ppm 3.834739324780E+02 o2ppm 1.331458061453E+05 n2ppm 7.687669232610E+05 co2ppm 3.402355479943E+04 h2oppm 3.936631411463E+04 hppm 0.000000000000E+00 h2ppm 9.361464258881E+03 oppm 0.000000000000E+00 nppm 0.000000000000E+00 ohppm 0.000000000000E+00 coppm 1.117681834415E+04 noppm 6.088871762233E+02 c2h2ppm 6.968746632825E+02 rnppm 1.476186484820E-05 pnppm 2.030617029262E+02 rppm 4.029167084315E+00 bppm 1.155016400041E-01 qppm 2.262676918207E+03 crank= 130.03 fluxin= 0.00000E+00 fluxout= 0.00000E+00 wheat= 1.11054E+09 swirl ratio= -8.61897E-01 crenb= 130.03 0.00000E+00 2.66328E+09 1.11054E+09 -3.04316E+09 0.00000E+00 cremi= 130.03 3.40236E+04 1.11768E+04 6.08887E+02 Indicated Gross Work= 297.556342722382 Nm Indicated power= 5.95112685444764 kW ISFC= 82.8750608183393 g/kWh Execution time= 14369.3 s (= 239.49 m, = 3.99 h) Total time from beginning= 14366.8 s (= 239.45 m, = 3.99 h) -xxxiv- Phụ Lục 8: Otape12.txt_Phun Tách 50%4CA50% -********************************* -grams of each species in entire mesh: c14h30 7.780374627207E-04 o2 9.447984508869E-02 n2 5.036821987820E-01 co2 4.152142858651E-02 h2o 1.937372968193E-02 h 0.000000000000E+00 h2 2.508925670668E-04 o 0.000000000000E+00 n 0.000000000000E+00 oh 0.000000000000E+00 co 4.861271881138E-03 no 3.427908726371E-04 c2h2 4.859465666957E-04 rn 1.788770185041E-10 pn 8.198049616017E-05 r 1.003992805504E-05 b 6.737863777220E-07 q 1.129625727784E-02 particle mass distribution: entrained= 0.00000E+00, and on surfaces numbered -1 to nmoving: 0.000000000000000E+000 4.243421823555173E-006 parts/million by volume of each species in region 1: c14h30ppm 1.680402107876E+02 o2ppm 1.265116257364E+05 n2ppm 7.703564372054E+05 co2ppm 4.042521081085E+04 h2oppm 4.607826362187E+04 hppm 0.000000000000E+00 h2ppm 5.332593878569E+03 oppm 0.000000000000E+00 nppm 0.000000000000E+00 ohppm 0.000000000000E+00 coppm 7.436396536120E+03 noppm 4.894782976359E+02 c2h2ppm 8.008596492443E+02 rnppm 1.277451849988E-05 pnppm 2.927322284285E+02 rppm 3.740881944159E+00 bppm 1.255265615861E-01 qppm 2.104495403532E+03 crank= 130.04 fluxin= 0.00000E+00 fluxout= 0.00000E+00 wheat= 1.17148E+09 swirl ratio= -8.63352E-01 crenb= 130.04 0.00000E+00 3.17098E+09 1.17148E+09 -3.27348E+09 0.00000E+00 cremi= 130.04 4.04252E+04 7.43640E+03 4.89478E+02 Indicated Gross Work= 321.088234605818 Nm Indicated power= 6.42176469211636 kW ISFC= 76.8013192083282 g/kWh Execution time= 14742.9 s (= 245.71 m, = 4.10 h) Total time from beginning= 14740.6 s (= 245.68 m, = 4.09 h) -xxxv- Phụ Lục 9: Otape12.txt_Phun Tách 35%4CA65% -********************************* -grams of each species in entire mesh: c14h30 1.454879867738E-03 o2 9.657234887287E-02 n2 5.036822934349E-01 co2 3.932767868085E-02 h2o 1.826000585102E-02 h 0.000000000000E+00 h2 3.193340976709E-04 o 0.000000000000E+00 n 0.000000000000E+00 oh 0.000000000000E+00 co 5.522556311414E-03 no 3.470485318436E-04 c2h2 4.925581758718E-04 rn 2.447081919849E-10 pn 6.621801953020E-05 r 1.015933247955E-05 b 7.205628147871E-07 q 1.145191246715E-02 parts/million by volume of each species in region 1: c14h30ppm 3.140297896562E+02 o2ppm 1.292334874050E+05 n2ppm 7.698795837719E+05 co2ppm 3.826567019785E+04 h2oppm 4.340250411536E+04 hppm 0.000000000000E+00 h2ppm 6.783081165867E+03 oppm 0.000000000000E+00 nppm 0.000000000000E+00 ohppm 0.000000000000E+00 coppm 8.442747232541E+03 noppm 4.952510535682E+02 c2h2ppm 8.112532170251E+02 rnppm 1.746503682314E-05 pnppm 2.363018867635E+02 rppm 3.783028219243E+00 bppm 1.341579019864E-01 qppm 2.132172960908E+03 crank= 130.05 fluxin= 0.00000E+00 fluxout= 0.00000E+00 wheat= 1.12102E+09 swirl ratio= -8.59546E-01 crenb= 130.05 0.00000E+00 2.99565E+09 1.12102E+09 -3.20952E+09 0.00000E+00 cremi= 130.05 3.82657E+04 8.44275E+03 4.95251E+02 Indicated Gross Work= 314.629206574485 Nm Indicated power= 6.29258413148969 kW ISFC= 78.3779747229603 g/kWh Execution time= 19232.5 s (= 320.54 m, = 5.34 h) Total time from beginning= 19229.9 s (= 320.50 m, = 5.34 h) -xxxvi- Phụ Lục 10: Otape12.txt_Phun Tách 80%2CA20% -********************************* -grams of each species in entire mesh: c14h30 1.669854983414E-03 o2 9.960448893028E-02 n2 5.036812842149E-01 co2 3.640314929666E-02 h2o 1.689010394762E-02 h 0.000000000000E+00 h2 3.161202886770E-04 o 0.000000000000E+00 n 0.000000000000E+00 oh 0.000000000000E+00 co 5.324817822370E-03 no 4.380215587979E-04 c2h2 3.721145432134E-04 rn 1.711568365475E-10 pn 4.480411755915E-05 r 1.092908885521E-05 b 5.959414050516E-07 q 1.203234166664E-02 particle mass distribution: entrained= 0.00000E+00, and on surfaces numbered -1 to nmoving: 1.775558545108192E-006 2.362282967706724E-005 parts/million by volume of each species in region 1: c14h30ppm 3.613007176648E+02 o2ppm 1.336126365982E+05 n2ppm 7.717351579341E+05 co2ppm 3.550555627258E+04 h2oppm 4.024320399728E+04 hppm 0.000000000000E+00 h2ppm 6.731013212249E+03 oppm 0.000000000000E+00 nppm 0.000000000000E+00 ohppm 0.000000000000E+00 coppm 8.160086102259E+03 noppm 6.265807326320E+02 c2h2ppm 6.143585358356E+02 rnppm 1.224507928562E-05 pnppm 1.602711241049E+02 rppm 4.079479146563E+00 bppm 1.112229284018E-01 qppm 2.245644058881E+03 crank= 130.06 fluxin= 0.00000E+00 fluxout= 0.00000E+00 wheat= 1.07222E+09 swirl ratio= -8.65184E-01 crenb= 130.06 0.00000E+00 2.73605E+09 1.07222E+09 -3.09202E+09 0.00000E+00 cremi= 130.06 3.55056E+04 8.16009E+03 6.26581E+02 Indicated Gross Work= 302.479559207401 Nm Indicated power= 6.08259118414802 kW ISFC= 81.5261701141643 g/kWh Execution time= 12569.2 s (= 209.49 m, = 3.49 h) Total time from beginning= 12566.8 s (= 209.45 m, = 3.49 h) -xxxvii- Phụ Lục 11: Otape12.txt_Phun Tách 80%6CA20% -********************************* -grams of each species in entire mesh: c14h30 1.442118210991E-03 o2 9.811790215952E-02 n2 5.036818922985E-01 co2 3.606723410240E-02 h2o 1.749353975416E-02 h 0.000000000000E+00 h2 3.765095095959E-04 o 0.000000000000E+00 n 0.000000000000E+00 oh 0.000000000000E+00 co 7.216136057179E-03 no 4.378689530573E-04 c2h2 3.303329660103E-04 rn 7.584915453051E-11 pn 7.120834202771E-05 r 1.110455641210E-05 b 7.201954034071E-07 q 1.218832059086E-02 parts/million by volume of each species in region 1: c14h30ppm 3.109978470985E+02 o2ppm 1.311847436395E+05 n2ppm 7.691928956834E+05 co2ppm 3.506199773817E+04 h2oppm 4.154362394029E+04 hppm 0.000000000000E+00 h2ppm 7.990435997023E+03 oppm 0.000000000000E+00 nppm 0.000000000000E+00 ohppm 0.000000000000E+00 coppm 1.102201898419E+04 noppm 6.242983064857E+02 c2h2ppm 5.435801936928E+02 rnppm 5.408596079267E-06 pnppm 2.538836206455E+02 rppm 4.131316138508E+00 bppm 1.339700021576E-01 qppm 2.267258758012E+03 crank= 130.14 fluxin= 0.00000E+00 fluxout= 0.00000E+00 wheat= 1.14875E+09 swirl ratio= -8.65371E-01 crenb= 130.14 0.00000E+00 2.69681E+09 1.14875E+09 -3.19197E+09 0.00000E+00 cremi= 130.14 3.50620E+04 1.10220E+04 6.24298E+02 Indicated Gross Work= 312.271471326140 Nm Indicated power= 6.24542942652281 kW ISFC= 78.9697499271228 g/kWh Execution time= 10444.3 s (= 174.07 m, = 2.90 h) Total time from beginning= 10441.6 s (= 174.03 m, = 2.90 h) -xxxviii- Phụ Lục 12: Otape12.txt_Phun Tách 80%8CA20% -********************************* -grams of each species in entire mesh: c14h30 1.489499754095E-03 o2 9.906820292339E-02 n2 5.036818671576E-01 co2 3.490883061001E-02 h2o 1.708681705210E-02 h 0.000000000000E+00 h2 3.938895886742E-04 o 0.000000000000E+00 n 0.000000000000E+00 oh 0.000000000000E+00 co 7.486592039778E-03 no 4.556541197601E-04 c2h2 4.111175301735E-04 rn 9.556235207628E-11 pn 9.055009595866E-05 r 1.114100932280E-05 b 7.936948662706E-07 q 1.232351897104E-02 parts/million by volume of each species in region 1: c14h30ppm 3.211441907012E+02 o2ppm 1.324257602558E+05 n2ppm 7.690212767175E+05 co2ppm 3.392831038554E+04 h2oppm 4.056868816072E+04 hppm 0.000000000000E+00 h2ppm 8.357418344421E+03 oppm 0.000000000000E+00 nppm 0.000000000000E+00 ohppm 0.000000000000E+00 coppm 1.143256614619E+04 noppm 6.495108619381E+02 c2h2ppm 6.763645169528E+02 rnppm 6.812770328375E-06 pnppm 3.227719879091E+02 rppm 4.143953423539E+00 bppm 1.476093620888E-01 qppm 2.291896862772E+03 crank= 130.03 fluxin= 0.00000E+00 fluxout= 0.00000E+00 wheat= 1.11332E+09 swirl ratio= -8.65774E-01 crenb= 130.03 0.00000E+00 2.65135E+09 1.11332E+09 -3.11902E+09 0.00000E+00 cremi= 130.03 3.39283E+04 1.14326E+04 6.49511E+02 Indicated Gross Work= 304.985469982541 Nm Indicated power= 6.09970939965083 kW ISFC= 80.8563109626555 g/kWh Execution time= 10575.4 s (= 176.26 m, = 2.94 h) Total time from beginning= 10573.2 s (= 176.22 m, = 2.94 h) -xxxix- Phụ Lục 13: Otape12.txt_Phun Tách 80%10CA20% -********************************* -grams of each species in entire mesh: c14h30 1.370577556943E-03 o2 9.774575043128E-02 n2 5.036818675945E-01 co2 3.742066985052E-02 h2o 1.783968575936E-02 h 0.000000000000E+00 h2 3.053952945915E-04 o 0.000000000000E+00 n 0.000000000000E+00 oh 0.000000000000E+00 co 5.860679510525E-03 no 5.347457754255E-04 c2h2 3.367215267285E-04 rn 9.630900886990E-11 pn 6.072649048530E-05 r 1.054276295794E-05 b 9.645603065300E-07 q 1.199013786603E-02 parts/million by volume of each species in region 1: c14h30ppm 2.961267552428E+02 o2ppm 1.309334081427E+05 n2ppm 7.706421405193E+05 co2ppm 3.644625307141E+04 h2oppm 4.244547446807E+04 hppm 0.000000000000E+00 h2ppm 6.493433230532E+03 oppm 0.000000000000E+00 nppm 0.000000000000E+00 ohppm 0.000000000000E+00 coppm 8.968543622995E+03 noppm 7.638584315312E+02 c2h2ppm 5.551369046898E+02 rnppm 6.880471951249E-06 pnppm 2.169199534835E+02 rppm 3.929697922857E+00 bppm 1.797645763337E-01 qppm 2.234595430764E+03 crank= 130.12 fluxin= 0.00000E+00 fluxout= 0.00000E+00 wheat= 1.16435E+09 swirl ratio= -8.63725E-01 crenb= 130.12 0.00000E+00 2.72665E+09 1.16435E+09 -3.23925E+09 0.00000E+00 cremi= 130.12 3.64463E+04 8.96854E+03 7.63858E+02 Indicated Gross Work= 316.925347075262 Nm Indicated power= 6.33850694150524 kW ISFC= 77.8101222498427 g/kWh Execution time= 9772.2 s (= 162.87 m, = 2.71 h) Total time from beginning= 9769.8 s (= 162.83 m, = 2.71 h) -xl- Phụ Lục 14: Áp Suất Xylanh (Không Phun)_Thực Nghiệm Pressure Angle PCYL1 (Deg) (bar) -134 0.73941 -133 0.74505 -132 0.75163 -131 0.75673 -130 0.76287 -129 0.76922 -128 0.77578 -127 0.78158 -126 0.78856 -125 0.79577 -124 0.80322 -123 0.81092 -122 0.81888 -121 0.82711 -120 0.83561 -119 0.8444 -118 0.85349 -117 0.86289 -116 0.8726 -115 0.88264 -114 0.89454 -113 0.90377 -112 0.91487 -111 0.92637 -110 0.93826 -109 0.94877 -108 0.96144 -107 0.97456 -106 0.98813 -105 1.00219 -104 1.01674 -103 1.03182 -102 1.04743 -101 1.06361 -100 1.08038 -99 1.09777 -98 1.11579 -97 1.13449 -96 1.15388 -95 1.17401 -94 1.19489 -93 1.21657 -92 1.23909 -91 -90 -89 -88 -87 -86 -85 -84 -83 -82 -81 -80 -79 -78 -77 -76 -75 -74 -73 -72 -71 -70 -69 -68 -67 -66 -65 -64 -63 -62 -61 -60 -59 -58 -57 -56 -55 -54 -53 -52 -51 -50 -49 -48 -47 -46 1.25908 1.28325 1.30837 1.33449 1.36164 1.3899 1.4193 1.4499 1.48177 1.51497 1.54957 1.58563 1.62324 1.66247 1.7034 1.74606 1.79052 1.83695 1.87841 1.92884 1.98159 2.0368 2.09461 2.15517 2.21864 2.2852 2.35502 2.42831 2.50527 2.58614 2.67115 2.76056 2.85466 2.95375 3.05813 3.16816 3.26725 3.38877 3.51707 3.6526 3.79578 3.94668 4.10629 4.27485 4.45338 4.64256 -xli- -45 -44 -43 -42 -41 -40 -39 -38 -37 -36 -35 -34 -33 -32 -31 -30 -29 -28 -27 -26 -25 -24 -23 -22 -21 -20 -19 -18 -17 -16 -15 -14 -13 -12 -11 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 4.84315 5.05595 5.28178 5.52157 5.77631 6.04706 6.33493 6.64114 6.91914 7.26271 7.6274 8.01348 8.42405 8.86099 9.32595 9.8206 10.3344 10.8836 11.4758 12.1043 12.7707 13.4758 14.2175 14.9601 15.9494 16.5948 17.5338 18.462 19.6089 20.4698 21.4897 22.6703 23.5757 24.5431 25.6474 26.5847 27.7821 28.5252 29.4252 30.1973 30.8872 31.4731 31.9085 32.2593 32.4481 32.4838 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 32.3375 32.0569 31.6363 31.0849 30.4152 29.6409 28.777 27.8394 26.8365 25.8084 24.7396 23.6529 22.5645 21.4867 20.4665 19.3306 18.3332 17.3723 16.4506 15.6913 14.7722 14.0202 13.2598 12.5416 11.8643 11.2263 10.6262 10.062 9.53181 9.08517 8.56095 8.12171 7.7096 7.37696 7.01146 6.66875 6.34731 6.04568 5.76251 5.51034 5.25267 5.01734 4.79586 4.5874 4.39107 4.20609 4.03176 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 3.86745 3.7125 3.56627 3.4282 3.29776 3.17444 3.05779 2.96279 2.85741 2.75756 2.66287 2.57302 2.48772 2.40669 2.32965 2.25635 2.18658 2.12015 2.05687 1.99656 1.93905 1.8842 1.83185 1.78189 1.73418 1.69498 1.65112 1.60918 1.56906 1.53067 1.4939 1.45869 1.42694 1.39056 1.35962 1.33411 1.30545 1.27794 1.25153 1.22615 1.20176 1.17831 1.15712 1.13283 1.11201 1.09479 1.0754 -xlii- 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 1.05672 1.03873 1.02139 1.00469 0.98858 0.97306 0.95808 0.94364 0.9297 0.91626 0.90328 0.89074 0.87864 0.86696 0.85567 0.84477 0.83572 0.82548 0.8156 0.80604 0.79681 0.78788 0.77925 0.7709 0.76283 0.75503 0.74749 0.74019 0.73314 0.72632 0.71973 0.71336 0.7072 0.70124 0.69629 0.6907 Phụ Lục 15: Áp Suất Xylanh (Có Phun)_Thực Nghiệm Pressure Angle PCYL1 (Deg) (Bar) -134 0.99588 -133 1.00273 -132 1.01084 -131 1.01851 -130 1.02678 -129 1.03496 -128 1.04327 -127 1.05323 -126 1.06166 -125 1.07185 -124 1.08182 -123 1.09213 -122 1.10278 -121 1.11379 -120 1.12517 -119 1.13693 -118 1.14732 -117 1.15983 -116 1.17276 -115 1.18612 -114 1.19994 -113 1.21423 -112 1.22901 -111 1.2443 -110 1.26011 -109 1.27647 -108 1.29339 -107 1.31091 -106 1.32904 -105 1.34781 -104 1.36724 -103 1.38736 -102 1.4082 -101 1.4298 -100 1.44892 -99 1.47199 -98 1.4959 -97 1.52069 -96 1.5464 -95 1.57308 -94 1.60075 -93 1.62947 -92 1.65929 -91 -90 -89 -88 -87 -86 -85 -84 -83 -82 -81 -80 -79 -78 -77 -76 -75 -74 -73 -72 -71 -70 -69 -68 -67 -66 -65 -64 -63 -62 -61 -60 -59 -58 -57 -56 -55 -54 -53 -52 -51 -50 -49 -48 -47 -46 1.69025 1.72241 1.75583 1.79055 1.82665 1.8642 1.90325 1.9439 1.9862 2.02385 2.06947 2.11701 2.16656 2.21822 2.27209 2.32827 2.38675 2.44776 2.51148 2.57808 2.64771 2.72055 2.79676 2.87654 2.96011 3.04767 3.13946 3.22169 3.32202 3.42734 3.53795 3.65418 3.77637 3.9049 4.04015 4.18256 4.33259 4.4907 4.65742 4.83332 5.01898 5.21468 5.41857 5.63597 5.86586 6.10902 -xliii- -45 -44 -43 -42 -41 -40 -39 -38 -37 -36 -35 -34 -33 -32 -31 -30 -29 -28 -27 -26 -25 -24 -23 -22 -21 -20 -19 -18 -17 -16 -15 -14 -13 -12 -11 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 6.32871 6.59899 6.88529 7.18866 7.51026 7.85128 8.2129 8.59658 9.00374 9.43587 9.89453 10.3793 10.8913 11.4341 12.0097 12.6198 13.2535 13.937 14.66 15.4239 16.23 16.9552 17.8417 18.7631 19.952 20.7572 21.8927 22.9045 24.202 25.0368 26.0726 27.3462 28.0956 29.0404 29.9226 30.7391 31.708 32.3802 33.1544 33.8763 34.5154 35.0532 35.4874 35.7987 35.9939 36.0692 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 36.4091 39.4658 46.7951 52.0042 56.9927 61.9578 65.2157 67.6706 69.0629 69.6142 68.9052 67.4938 65.5776 63.2791 61.0504 58.5675 56.1411 53.7613 51.27 49.1432 46.6315 44.4627 42.5514 40.5199 38.6133 36.9318 35.1568 33.5193 31.9607 30.6092 29.086 27.7702 26.5535 25.433 24.2382 23.2605 22.2089 21.3369 20.3849 19.628 18.7471 17.9907 17.2758 16.6002 15.9614 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 15.441 14.8634 14.3175 13.7972 13.3028 12.8335 12.3882 11.9659 11.5651 11.1859 10.8244 10.4802 10.1522 9.83949 9.54099 9.25613 8.98516 8.72625 8.5139 8.27677 8.05043 7.83421 7.62746 7.42959 7.24055 7.05886 6.88455 6.71745 6.5574 6.40415 6.25739 6.11674 5.9821 5.85252 5.74416 5.61253 5.50063 5.38929 5.28219 5.17922 5.08024 4.98529 4.89379 4.80572 4.73267 -xliv- 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 4.64814 4.56446 4.48811 4.41439 4.34321 4.27453 4.20828 4.15352 4.09159 4.03183 3.97418 3.91854 3.86481 3.8129 3.76273 3.71421 3.66725 3.62183 3.57788 3.53549 3.49438 3.45462 3.41616 3.37899 3.34305 3.31319 3.27941 3.24675 3.21513 3.18451 3.15488 3.12619 3.09841 3.07152 3.04561 3.02045 2.99612 2.9726 2.94987 2.92791 ... TÀI: NGHIÊN CỨU ẢNH HƯỞNG Q TRÌNH PHUN NHIỀU GIAI ĐOẠN ĐẾN KHÍ THẢI CỦA ĐỘNG CƠ VIKYNO RV1 25-2 BẰNG PHƯƠNG PHÁP MÔ PHỎNG NHIỆM VỤ CỦA LUẬN VĂN: 2.1 Tổng quan tình hình nghiên cứu tác động q trình. .. tượng nghiên cứu Mục tiêu ñề tài: nghiên cứu ảnh hưởng trình phun nhiều giai đoạn đến khí thải động VIKYNO RV1 25-2 phương pháp mơ theo định hướng giảm nhiễm mơi trường - Xác định lưu lượng phun. .. buồng cháy, cải tiến q trình phun Trong đó, cải tiến q trình phun, cụ thể phun nhiều giai ñoạn biện pháp khả thi cho động Vì việc nghiên cứu ảnh hưởng q trình phun nhiều giai đoạn nhằm làm cho hệ