Tổng quan về luận án
Sự cạn kiệt nhanh chóng của trữ lượng nhiên liệu hóa thạch cùng áp lực giảm thiểu phát thải khí nhà kính theo Thỏa thuận Paris (COP21) đã đặt ngành công nghiệp ô tô toàn cầu trước yêu cầu cấp thiết về việc tìm kiếm năng lượng tái tạo thay thế. Mặc dù xu hướng điện hóa phương tiện đang phát triển, động cơ đốt trong (ICE) được dự báo vẫn là nguồn động lực chiếm ưu thế đến năm 2040 (BP, 2016). Trong bối cảnh đó, ethanol sinh học nổi lên như giải pháp hàng đầu nhờ trị số octan cao (RON = 108,6 so với RON = 92 của xăng khoáng), nhiệt ẩn hóa hơi lớn ($L = 840\text{ kJ/kg}$ so với $305\text{ kJ/kg}$) và chu trình cacbon trung tính (IPCC, 2013). Tuy nhiên, việc ứng dụng ethanol tỷ lệ cao trên động cơ đánh lửa cưỡng bức thương mại đang gặp rào cản lớn về sự suy giảm công suất ở hỗn hợp nghèo, hiện tượng phân tầng nhiên liệu khi hòa trộn tĩnh và khả năng bay hơi kém khi khởi động nguội.
Khoảng trống nghiên cứu (Research Gap) cốt lõi được xác định là: Hầu hết các công trình trước đây chỉ tập trung khảo sát hỗn hợp xăng - ethanol hòa trộn sẵn (blends) ở tỷ lệ thấp (E5 - E20) trên động cơ nguyên bản hoặc can thiệp ECU hạn chế (Mustafa Koç và cộng sự, 2009; Costagliola và cộng sự, 2013). Chưa có nghiên cứu nào tích hợp toàn diện mô hình 3D-CFD đa pha để làm sáng tỏ bản chất nhiệt động lực học, quá trình bay hơi cục bộ và cơ chế lan truyền màng lửa rối của các cấu hình phun riêng rẽ (Dual Injection: Gasoline Port Injection - Ethanol Direct Injection/Ethanol Port Injection) nhằm giải phóng giới hạn tỷ lệ ethanol lên tới E100 theo chế độ tải.
Luận án thiết lập 4 câu hỏi nghiên cứu ($RQ$) và 4 giả thuyết khoa học ($H$):
- $RQ_1$: Giới hạn tỷ lệ phối trộn ethanol tối ưu trên động cơ phun xăng điện tử thương mại không thay đổi phần cứng/phần mềm là bao nhiêu?
- $H_1$: Ở tỷ số nén $\varepsilon = 9{,}5$, tỷ lệ ethanol trong hỗn hợp hòa trộn sẵn bị giới hạn ở mức $\le\text{E20}$; vượt quá ngưỡng này, công suất và mô-men xoắn suy giảm do hòa khí bị nhạt hóa cục bộ ($f < 1$).
- $RQ_2$: Nhiệt độ môi chất và thời điểm phun tác động như thế nào đến tốc độ bay hơi và phân bố nồng độ hơi ethanol trong đường nạp và buồng cháy?
- $H_2$: Do nhiệt ẩn hóa hơi cao và áp suất hơi bão hòa thấp, ethanol đòi hỏi nhiệt độ môi chất $T \ge 345\text{ K}$ hoặc thời điểm phun sớm để đạt độ đồng nhất hòa khí tương đương xăng.
- $RQ_3$: Cấu hình phun riêng rẽ (GPI-EDI và GPI-EPI) có khắc phục được hiện tượng hòa trộn kém và sụt giảm áp suất cháy của hỗn hợp ethanol tỷ lệ cao không?
- $H_3$: Phun riêng rẽ xăng trên đường nạp kết hợp phun trực tiếp ethanol (GPI-EDI) tạo ra sự phân tầng octan và làm mát buồng cháy hiệu quả, cho phép tăng góc đánh lửa sớm $\theta_s$, tối ưu hóa áp suất chỉ thị cực đại ($p_{\max}$) và giảm phát thải.
- $RQ_4$: Cơ chế hình thành phát thải $\text{NO}_x$ và $\text{CO/HC}$ biến thiên như thế nào theo tỷ lệ ethanol và cấu hình phun?
- $H_4$: Hàm lượng oxy nội tại trong ethanol thúc đẩy phản ứng oxy hóa hoàn toàn giúp giảm mạnh $\text{CO}$ và $\text{HC}$, trong khi phát thải $\text{NO}_x$ được kiểm soát nhờ hiệu ứng làm mát nạp từ nhiệt ẩn hóa hơi của ethanol.
Khung lý thuyết của luận án xây dựng trên nền tảng cơ học chất lưu tính toán (3D-CFD), kết hợp hệ phương trình Navier-Stokes trung bình Reynolds (RANS), mô hình rối $k\text{-}\varepsilon$ tiêu chuẩn (Launder & Spalding, 1974), mô hình ngọn lửa rối Zimont (Zimont, 2000), mô hình bay hơi Dukowicz (1979), mô hình phân rã giọt TAB (O'Rourke & Amsden, 1987) và cơ chế động học hình thành $\text{NO}_x$ nhiệt mở rộng Zeldovich (1946).
Đóng góp đột phá của nghiên cứu bao gồm việc định lượng hóa chính xác giới hạn vận hành của động cơ Daewoo A16DMS (1,6 lít, 4 xilanh, $\varepsilon = 9{,}5$) với xăng sinh học, đồng thời chứng minh giải pháp phun riêng rẽ GPI-EDI cho phép vận hành linh hoạt đến E100, cải thiện hiệu suất nhiệt có ích ($\eta_e$) lên đến $4{,}2%$, giảm phát thải $\text{CO}$ từ $15\text{--}35%$ và $\text{HC}$ từ $20\text{--}40%$ so với xăng khoáng E0. Nghiên cứu cung cấp luận cứ khoa học thực chứng phục vụ Lộ trình áp dụng tỷ lệ phối trộn nhiên liệu sinh học theo Quyết định số 53/2012/QĐ-TTg của Thủ tướng Chính phủ.
Literature Review và Positioning
Nghiên cứu ứng dụng cồn sinh học trên động cơ đánh lửa cưỡng bức (SI) đã trải qua ba làn sóng học thuật chính:
Dòng nghiên cứu 1: Đánh giá tính năng hỗn hợp hòa trộn sẵn (Blends) trên động cơ truyền thống. J. Mantilla và cộng sự (2010) khi nghiên cứu trên động cơ sử dụng bộ chế hòa khí đã chỉ ra rằng tốc độ cháy đạt cực đại tại tỷ lệ E25, tương ứng với hệ số tương đương $f \approx 1{,}0$, nhưng suy giảm nhanh khi tăng hàm lượng ethanol do hòa khí bị loãng. Đồng thuận với kết luận này, M. Al-Hasan (2003) trên động cơ TOYOTA TERCEL-3A khẳng định công suất và mô-men có ích tăng dần từ E0 đến E20 nhưng giảm mạnh ở tỷ lệ $> \text{E20}$ nếu không tăng lưu lượng cấp nhiên liệu. Mustafa Koç và cộng sự (2009) khi thử nghiệm trên động cơ tỷ số nén biến thiên ($\varepsilon = 10\text{--}11$) chứng minh rằng hỗn hợp E50 và E85 giúp tăng mô-men xoắn nhờ nhiệt ẩn hóa hơi lớn và chống kích nổ tốt, song suất tiêu hao nhiên liệu có ích ($g_e$) tăng tương ứng do nhiệt trị thấp của ethanol ($Q_t = 26{,}8\text{ MJ/kg}$ so với $42{,}5\text{ MJ/kg}$ của xăng).
Dòng nghiên cứu 2: Can thiệp hệ thống điều khiển điện tử (ECU). Costagliola và cộng sự (2013) phát hiện rằng trên động cơ phun xăng điện tử (MPI), với tỷ lệ pha dưới E20, ECU nguyên bản có thể tự thích ứng bằng cách tăng thời gian phun; tuy nhiên, ở chế độ toàn tải ($bmep = 6{,}6\text{ bar}$, $3000\text{ rpm}$), việc sử dụng E85 đòi hỏi góc mở vòi phun kéo dài tới $180^\circ\text{CA}$, làm vượt ngưỡng giới hạn điều khiển của ECU gốc. Tại Việt Nam, nhóm nghiên cứu của Phạm Minh Tuấn, Trần Anh Trung và Nguyễn Khánh Tùng (2015) đã phát triển bộ ECU phụ điều khiển tự động hệ số dư lượng không khí $\lambda = 1$ và tối ưu hóa góc đánh lửa sớm, cho phép động cơ Toyota 1NZ-FE chạy ổn định đến E100 nhưng vẫn tồn tại nguy cơ phân tầng nhiên liệu và khó khởi động ở nhiệt độ thấp.
Dòng nghiên cứu 3: Phun riêng rẽ và phân tầng nhiên liệu (Dual Injection & Octane On-Demand). Nhằm khắc phục hạn chế của việc trộn tĩnh, Di Iorio và cộng sự (2014) cùng Catapano và cộng sự (2013) đã tiên phong nghiên cứu kết hợp phun xăng trên đường nạp (GPI) và phun ethanol trực tiếp (EDI). Tuy nhiên, các tác giả ghi nhận áp suất buồng cháy của cấu hình DI ethanol thuần túy đôi khi thấp hơn PI do thời gian bay hơi ngắn và hiệu ứng làm mát quá mức gây dập tắt màng lửa cục bộ. Ngược lại, Zhuang và cộng sự (2014) chứng minh rằng khi tỷ lệ năng lượng ethanol (EER) đạt $15\text{--}35%$, cứ mỗi $2\text{--}3%$ tăng EER cho phép tăng góc đánh lửa sớm thêm $2^\circ\text{CA}$ mà không xảy ra kích nổ, tối ưu hóa quá trình giải phóng nhiệt.
┌─────────────────────────────────────────────────────────┐
│ TIẾN TRÌNH HỌC THUẬT │
└────────────────────────────┬────────────────────────────┘
│
┌───────────────────────────────┴───────────────────────────────┐
▼ ▼
┌───────────────────────────────┐ ┌───────────────────────────────┐
│ HÒA TRỘN TĨNH (BLENDS) │ │ ECU PHỤ & DUAL-FUEL │
│ - Mantilla (2010): E25 max │ │ - Costagliola (2013): │
│ - Al-Hasan (2003): E20 limit │ │ E85 vượt ngưỡng ECU gốc │
│ - Koç (2009): E50/E85 tăng ge│ │ - Nguyễn Khánh Tùng (2015): │
│ - Hsieh (2002): Giảm CO/HC │ │ ECU phụ cho E100 │
└──────────────┬────────────────┘ └───────────────┬───────────────┘
│ │
└───────────────────────────────┬──────────────────────────────┘
│
▼
┌──────────────────────────────────┐
│ RESEARCH GAP & POSITIONING │
│ - Thiếu mô hình 3D-CFD đa pha │
│ - Hiện tượng bay hơi/làm mát │
│ - Cấu hình GPI-EDI vs GPI-EPI │
│ - Tối ưu động học màng lửa rối │
└────────────────┬─────────────────┘
│
▼
┌──────────────────────────────────┐
│ ĐÓNG GÓP CỦA LUẬN ÁN │
│ Giải pháp 3D-CFD Ansys Fluent │
│ kết hợp thực nghiệm APA204: │
│ Phun riêng rẽ tối ưu E0 - E100 │
└──────────────────────────────────┘
Tranh luận học thuật cốt lõi tồn tại giữa hai quan điểm: (1) Nhóm nghiên cứu thực nghiệm truyền thống cho rằng ethanol tỷ lệ cao gây suy giảm công suất do nhiệt trị thấp và hòa khí loãng; (2) Nhóm nghiên cứu nhiệt động học động lực cho rằng nếu tận dụng tốt tỷ số nén cao và phân tầng nồng độ qua phun riêng rẽ, hiệu suất chỉ thị sẽ vượt trội nhờ nhiệt ẩn hóa hơi và tốc độ cháy của ethanol.
Luận án này định vị chính xác tại giao điểm của hai trường phái: Sử dụng công cụ 3D-CFD chuyên sâu kết hợp kiểm chứng thực nghiệm nghiêm ngặt để chứng minh rằng: Sự suy giảm công suất không phải do bản chất của ethanol mà do giới hạn của phương thức cung cấp nhiên liệu truyền thống; giải pháp phun riêng rẽ linh hoạt (GPI-EDI/GPI-EPI) là chìa khóa hiện thực hóa toàn bộ tiềm năng nhiệt động lực học của ethanol.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đóng góp quan trọng vào lý thuyết quá trình cháy động cơ đốt trong thông qua việc mở rộng và tinh chỉnh các mô hình động lực học nhiệt:
- Mở rộng mô hình ngọn lửa rối Zimont cho hòa khí đa cấu tử không đồng nhất: Xác lập lại tương quan giữa vận tốc ngọn lửa chảy tầng ($U_l$) và cấu trúc màng lửa nếp nhăn trong điều kiện hòa khí phân tầng xăng/ethanol, chứng minh sự hiện diện của ethanol làm tăng cường gradient nồng độ cục bộ và gia tốc tốc độ lan truyền màng lửa rối $U_t = A \cdot u' \cdot (U_l / \chi)^{1/2} \cdot \tau_t^{1/4}$.
- Lý thuyết hóa hiện tượng tương tác màng tia phun - thành buồng cháy (Film-Wall interaction): Ứng dụng mô hình Weber điều chỉnh với chuẩn số va đập $K = \text{Oh} \cdot \text{Re}_D^{1{,}25}$, xác định chính xác ngưỡng chuyển pha từ lắng đọng màng chất lỏng ($K < 57{,}7$) sang phân rã phản xạ ($K > 57{,}7$) đối với giọt ethanol có độ nhớt và sức căng bề mặt khác biệt so với hydrocarbon mạch thẳng.
- Lý thuyết phân tầng trị số octan động (Dynamic Octane Stratification): Thiết lập các mệnh đề khoa học về việc điều khiển tỷ lệ ethanol theo thời gian thực ($SOC, EOC$) để triệt tiêu hiện tượng kích nổ tự kích mà không cần tăng chỉ số octan của toàn bộ khối nhiên liệu nạp.
Khung phân tích độc đáo
Khung phân tích tích hợp ba trục tiếp cận lý thuyết nền tảng:
- Cơ học chất lưu tính toán đa pha (Eulerian-Lagrangian Framework): Pha khí được giải bằng hệ phương trình RANS trên lưới động cơ biến dạng theo chuyển động piston/xúp-páp; pha lỏng được theo dõi qua các hạt phân tán rời rạc (Discrete Phase Model - DPM).
- Lý thuyết bay hơi và phân rã giọt tương hỗ: Tích hợp đồng thời mô hình Dukowicz cho truyền nhiệt/khối đối lưu và mô hình TAB cho quá trình phân rã thứ cấp dưới tác dụng của lực khí động học bên ngoài, sức căng bề mặt và lực cản nhớt.
- Động học hóa học phản ứng khử và oxy hóa mở rộng: Mô phỏng cơ chế cháy qua phương trình Arrhenius kết hợp mô hình Laminar Finite-Rate và cơ chế Zeldovich mở rộng với 3 phản ứng thuận nghịch then chốt hình thành $\text{NO}$ nhiệt.
┌──────────────────────────────────────────────────────────────┐
│ KHUNG PHÂN TÍCH 3D-CFD │
└──────────────────────────────┬───────────────────────────────┘
│
┌─────────────────────────────────┼─────────────────────────────────┐
▼ ▼ ▼
┌──────────────────────┐ ┌──────────────────────┐ ┌──────────────────────┐
│ PHA KHÍ LIÊN TỤC │ │ PHA RỜI RẠC DPM │ │ ĐỘNG HỌC HÓA HỌC │
│ - RANS k-ε Std │ │ - Phun tia DPM │ │ - Finite-Rate │
│ - Bảo toàn khối │ │ - Bay hơi Dukowicz │ │ Arrhenius │
│ lượng, động lượng │ │ - Phân rã giọt TAB │ │ - Ngọn lửa rối │
│ - Phương trình │ │ - Film-Wall Weber │ │ Zimont │
│ năng lượng │ │ (K = 57.7) │ │ - Zeldovich NOx │
└──────────┬───────────┘ └──────────┬───────────┘ └──────────┬───────────┘
│ │ │
└─────────────────────────────┼─────────────────────────────┘
│
▼
┌──────────────────────────────────────────────────────────────┐
│ HIỆU CHỈNH & KIỂM CHỨNG (VALIDATION) │
│ Băng thử APA204 + Động cơ Daewoo A16DMS + KEG-500 │
│ Độ lệch áp suất buồng cháy p_max giữa CFD & EXP < 3.5% │
└──────────────────────────────────────────────────────────────┘
Điều kiện biên xác định: Động cơ Daewoo A16DMS; đường kính xilanh $D = 79\text{ mm}$; hành trình piston $S = 81{,}5\text{ mm}$; dung tích công tác $1598\text{ cm}^3$; tỷ số nén $\varepsilon = 9{,}5$; nhiệt độ khí nạp $T_{kn} = 300\text{--}345\text{ K}$; dải tốc độ $n = 1000\text{--}4250\text{ rpm}$; góc mở bướm ga $\text{THA} = 10%,\ 30%,\ 50%,\ 70%$.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu tuân thủ thế giới quan thực chứng (Positivism), kết hợp chặt chẽ giữa mô phỏng số học 3D-CFD độ phân giải cao và thực nghiệm kiểm chứng trên băng thử công nghiệp. Thiết kế đa tầng bao gồm:
- Cấp độ vi mô: Mô phỏng hành vi động học giọt nhiên liệu, tốc độ bay hơi $E_r$, phân bố nồng độ hơi xăng ($G_v$) và hơi ethanol ($E_v$) trong trường dòng rối.
- Cấp độ vĩ mô: Phân tích diễn biến góc quay trục khuỷu ($\theta$), áp suất buồng cháy trong xilanh ($p$), tỷ lệ hòa khí cháy (Mass Fraction Burn - MFB) và sự hình thành các chất ô nhiễm.
- Cấp độ hệ thống: Đánh giá đặc tính kinh tế - kỹ thuật toàn diện ($N_e, M_e, g_e, \eta_e$) và nồng độ phát thải đầu ra ($\text{CO}, \text{HC}, \text{CO}_2, \text{NO}_x$).
┌────────────────────────────────────────────────────────────────────┐
│ SƠ ĐỒ PHƯƠNG PHÁP NGHIÊN CỨU │
└─────────────────────────────────┬──────────────────────────────────┘
│
┌───────────────────────┴───────────────────────┐
▼ ▼
┌───────────────────────────────┐ ┌───────────────────────────────┐
│ NGHIÊN CỨU MÔ PHỎNG │ │ NGHIÊN CỨU THỰC NGHIỆM │
│ - Phần mềm Ansys Fluent │ │ - Động cơ Daewoo A16DMS │
│ - Xây dựng hình học 3D │ │ - Băng thử APA204/E/0943 │
│ (Port 1 phía & 2 phía) │ │ - Phân tích khí thải KEG-500 │
│ - Lưới động dynamic mesh │ │ - Cảm biến áp suất Kistler │
│ - Khai báo DPM, Zimont, NOx │ │ - Quy trình Stationary Step │
└───────────────┬───────────────┘ └───────────────┬───────────────┘
│ │
└───────────────────────┬───────────────────────┘
│
▼
┌───────────────────────────────────────────────┐
│ SO SÁNH & HIỆU CHỈNH MÔ HÌNH │
│ - Sai số p_max < 3.5% │
│ - Sai số góc xuất hiện p_max < 1.0 deg CA │
└───────────────────────┬───────────────────────┘
│
▼
┌───────────────────────────────────────────────┐
│ MỞ RỘNG MÔ PHỎNG CHUYÊN SÂU │
│ - Khảo sát cấu hình: PI vs DI, GPI-EDI │
│ - Tối ưu hóa thời điểm phun, nhiệt độ nạp │
│ - Thiết lập bản đồ tỷ lệ ethanol tối ưu │
└───────────────────────────────────────────────┘
Quy trình nghiên cứu thực nghiệm
Hệ thống thực nghiệm được chuẩn hóa tại phòng thí nghiệm động cơ đốt trong với các thiết bị:
- Động cơ thử nghiệm: Daewoo A16DMS, 4 kỳ, 4 xilanh thẳng hàng, trục cam đôi DOHC 16 van, phun xăng đa điểm trên đường nạp, hệ thống điều khiển điện tử ECU nguyên bản.
- Băng thử công suất: Băng thử xoay chiều APA204/E/0943 cho phép điều khiển chính xác tốc độ và mô-men xoắn theo quy trình bước tĩnh (Stationary Step).
- Hệ thống đo phát thải: Thiết bị phân tích khí xả tự động KEG-500 đo nồng độ $\text{CO}$ ($0\text{--}9{,}99%$, độ phân giải $0{,}01%$), $\text{CO}_2$ ($0\text{--}20{,}0%$, độ phân giải $0{,}1%$), $\text{HC}$ ($0\text{--}9999\text{ ppm}$, độ phân giải $1\text{ ppm}$), $\text{NO}_x$ ($0\text{--}5000\text{ ppm}$, độ phân giải $1\text{ ppm}$) và $\text{O}_2$ ($0\text{--}25{,}0%$).
- Nhiên liệu khảo sát: Xăng cơ sở RON92 (E0) và các hỗn hợp pha chế ethanol nguyên chất (E10, E15, E20, E25, E30, E40, E50, E100).
- Quy trình đo lường: Mỗi điểm đo được thực hiện lặp lại 3 lần sau khi động cơ đạt trạng thái ổn định nhiệt (nhiệt độ nước làm mát $85 \pm 2^\circ\text{C}$, nhiệt độ dầu bôi trơn $90 \pm 3^\circ\text{C}$); dữ liệu được thu thập qua card ADC với tần số lấy mẫu cao đồng bộ theo tín hiệu xung góc quay trục khuỷu ($0{,}1^\circ\text{CA}$).
Thiết lập mô phỏng số 3D-CFD
Mô hình hình học 3D chi tiết của buồng cháy, nắp máy, xúp-páp nạp/thải và đường nạp (cả cấu hình 1 phía và 2 phía đối xứng) được tạo lập và chia lưới dynamic mesh trong Ansys Fluent.
- Thiết lập điều kiện ban đầu: Áp suất nạp $p_0 = 1\text{ bar}$, nhiệt độ nạp ban đầu $T_0 = 300\text{--}345\text{ K}$, thành phần lưu chất ban đầu xác định theo tỷ lệ hòa khí quy định.
- Hiệu chuẩn mô hình (Model Calibration): Kết quả so sánh diễn biến áp suất buồng cháy giữa mô phỏng và thực nghiệm ở chế độ $3250\text{ rpm} - 50%\text{THA}$ cho thấy sự trùng khớp rất cao:
- Với E0: $p_{\max\text{ (CFD)}} = 30{,}48\text{ bar}$ tại $14{,}6^\circ\text{CA}$ sau ĐCT so với $p_{\max\text{ (EXP)}} = 30{,}78\text{ bar}$ tại $14{,}2^\circ\text{CA}$ (sai số $0{,}97%$).
- Với E10: $p_{\max\text{ (CFD)}} = 30{,}52\text{ bar}$ so với $p_{\max\text{ (EXP)}} = 31{,}15\text{ bar}$ (sai số $2{,}02%$).
- Với E20: $p_{\max\text{ (CFD)}} = 31{,}74\text{ bar}$ so với $p_{\max\text{ (EXP)}} = 31{,}20\text{ bar}$ (sai số $1{,}73%$).
- Với E40: $p_{\max\text{ (CFD)}} = 28{,}74\text{ bar}$ so với $p_{\max\text{ (EXP)}} = 27{,}81\text{ bar}$ (sai số $3{,}34%$).
Tất cả sai số đều nằm dưới ngưỡng $3{,}5%$, khẳng định độ tin cậy tuyệt đối của mô hình mô phỏng để triển khai các phân tích mở rộng.
Phát hiện đột phá và implications
Những phát hiện then chốt
- Ngưỡng tối ưu của hỗn hợp hòa trộn sẵn trên động cơ nguyên bản: Kết quả thực nghiệm khẳng định tỷ lệ ethanol tối ưu về công suất ($N_e$) và mô-men xoắn ($M_e$) trên động cơ Daewoo A16DMS là từ E10 đến E20. Tại $50%\text{THA}$ và $3250\text{ rpm}$, mô-men có ích tăng từ $98{,}4\text{ Nm}$ (E0) lên $102{,}8\text{ Nm}$ (E15, tăng $4{,}47%$) và $101{,}6\text{ Nm}$ (E20). Khi tỷ lệ vượt quá E20, do ECU không thể tự động kéo dài thời gian phun vượt mức giới hạn, hệ số tương đương giảm mạnh ($f_{bl} = 0{,}85$ ở E40), dẫn đến hiện tượng cháy nghèo quá mức, làm $M_e$ giảm sâu xuống còn $89{,}2\text{ Nm}$ ở E40 (giảm $9{,}35%$ so với E0).
- Quy luật biến thiên hiệu suất nhiệt và tiêu hao nhiên liệu: Do nhiệt trị thấp của ethanol, suất tiêu hao nhiên liệu có ích $g_e$ tăng đều theo hàm lượng ethanol (ở E20 tăng khoảng $5\text{--}8%$). Tuy nhiên, hiệu suất có ích $\eta_e$ lại tăng tuyến tính từ $26{,}8%$ (E0) lên đỉnh điểm $29{,}4%$ (E20) ở chế độ tải trung bình. Điều này chứng minh quá trình cháy được cải thiện nhờ oxy tự do trong ethanol và tốc độ cháy nhanh hơn.
- Quy luật phát thải khí ô nhiễm: Nồng độ phát thải $\text{CO}$ và $\text{HC}$ giảm mạnh khi tăng tỷ lệ ethanol từ E0 lên E20 trên mọi chế độ tải (ở $50%\text{THA}$, $\text{CO}$ giảm từ $1{,}85%$ xuống $0{,}62%$, giảm $66{,}5%$; $\text{HC}$ giảm từ $245\text{ ppm}$ xuống $140\text{ ppm}$, giảm $42{,}8%$). Ngược lại, phát thải $\text{CO}_2$ tăng nhẹ do quá trình cháy kiệt hydrocarbon. Phát thải $\text{NO}_x$ có xu hướng tăng nhẹ từ E0 đến E20 ở tải cao do tốc độ cháy nhanh làm nhiệt độ đỉnh buồng cháy tăng, nhưng giảm rõ rệt ở các tỷ lệ pha cao ($> \text{E30}$) nhờ nhiệt ẩn hóa hơi lớn của ethanol làm mát môi chất đầu quá trình nén.
- Nghịch lý bay hơi khi phun trực tiếp (DI) ethanol thuần túy: Kết quả 3D-CFD chỉ ra rằng khi phun trực tiếp E100 ở nhiệt độ khí nạp thấp ($T_{kn} < 345\text{ K}$), tốc độ bay hơi của ethanol thấp hơn đáng kể so với xăng do nhiệt độ sôi cao và nhiệt ẩn hóa hơi cực lớn ($840\text{ kJ/kg}$), dẫn đến tạo màng chất lỏng dày trên đỉnh piston và gây cháy không kiệt cục bộ.
- Ưu thế đột phá của giải pháp phun riêng rẽ (GPI-EDI): Cấu hình kết hợp phun xăng trên đường nạp và phun ethanol trực tiếp trong buồng cháy (GPI-EDI) giải quyết triệt để nghịch lý trên. Xăng phun trên đường nạp bay hơi hoàn toàn tạo hòa khí mồi đồng nhất, trong khi ethanol phun trực tiếp trong quá trình nén tận dụng nhiệt nén để bay hơi nhanh, làm mát khối khí nạp, tăng khối lượng nạp và nâng cao áp suất cháy cực đại $p_{\max}$ lên $36{,}8\text{ bar}$ (so với $32{,}1\text{ bar}$ của cấu hình phun hỗn hợp trên đường nạp ở cùng tỷ lệ E50).
Áp suất buồng cháy p (bar)
^
40 │ ── GPI-EDI (Phun riêng rẽ)
│ ─── ··· PI_Blend (Trộn sẵn đường nạp)
30 │ ──── ····· --- E0 (Xăng thương mại)
│ ─── ····
20 │ ─── ···
│ ─── ···
10 │ ─── ···
│ ─── ··
0 └───┴────┴────┴────┴────┴────┴────┴────┴────> Góc trục khuỷu (°CA)
-60 -40 -20 0 20 40 60 80
ĐCT
Implications đa chiều
- Về mặt học thuật: Cung cấp bộ dữ liệu thực nghiệm và mô phỏng 3D-CFD chi tiết, làm sáng tỏ cơ chế nhiệt động học vi mô của quá trình bay hơi và cháy phân tầng nhiên liệu cồn sinh học.
- Về mặt kỹ thuật công nghệ: Cung cấp giải pháp kết cấu đường nạp 2 phía và cấu hình điều khiển vòi phun kép, cho phép các kỹ sư ô tô chuyển đổi động cơ truyền thống sang động cơ nhiên liệu linh hoạt (Flexible Fuel Vehicle - FFV) với chi phí thấp nhất.
- Về mặt chính sách và môi trường: Khẳng định tính khả thi kỹ thuật vững chắc cho việc mở rộng lưu hành xăng E10 và tiến tới E20 tại Việt Nam, đóng góp trực tiếp vào mục tiêu cắt giảm khí nhà kính quốc gia và đảm bảo an ninh năng lượng.
Limitations và Future Research
Mặc dù đạt được những kết quả đột phá, luận án ghi nhận một số giới hạn nghiên cứu:
- Nghiên cứu thực nghiệm chỉ tiến hành trên một dòng động cơ nạp khí tự nhiên đại diện (Daewoo A16DMS) với tỷ số nén cố định $\varepsilon = 9{,}5$, chưa tiến hành thực nghiệm trực tiếp trên động cơ có tăng áp (Turbocharger) hoặc hệ thống điều khiển van biến thiên (VVT).
- Chưa đánh giá sâu tác động dài hạn của ethanol tỷ lệ cao đến độ bền mài mòn hóa học của kim phun, xéc-măng và đường ống dẫn nhiên liệu trong điều kiện vận hành thực tế tại Việt Nam (độ ẩm không khí cao).
- Mô hình động học hóa học $\text{NO}_x$ sử dụng cơ chế Zeldovich mở rộng kết hợp ngọn lửa rối Zimont dù đạt độ chính xác cao nhưng đã bỏ qua một số phản ứng trung gian phức tạp của gốc $\text{N}_2\text{O}$ ở dải áp suất cực cao.
Các hướng nghiên cứu tiếp theo được đề xuất:
- Mở rộng mô hình 3D-CFD và thực nghiệm trên động cơ phun xăng trực tiếp tăng áp cỡ nhỏ (Downsizing Turbocharged GDI) sử dụng hệ thống Octane-on-Demand.
- Nghiên cứu thuật toán điều khiển thích nghi thời gian thực (Adaptive Real-Time ECU) cho hệ thống phun riêng rẽ dựa trên cảm biến nhận diện nồng độ cồn quang học hoặc cảm biến áp suất trong buồng cháy.
- Đánh giá chu trình phát thải thực tế khi lái xe (Real Driving Emissions - RDE) và phân tích vòng đời phát thải (Life Cycle Assessment - LCA) toàn diện từ khâu sản xuất sinh khối sắn/mía đến khâu đốt cháy trong động cơ.
Tác động và ảnh hưởng
Luận án tạo ra tác động lan tỏa mạnh mẽ trên nhiều phương diện:
- Đóng góp học thuật: Mở ra hướng tiếp cận mới trong mô phỏng dòng chảy rối phản ứng cháy đa nhiên liệu; dự kiến tạo nguồn trích dẫn quan trọng cho các nghiên cứu về nhiên liệu sinh học và công nghệ động cơ đốt trong sạch.
- Chuyển giao công nghiệp: Cung cấp giải pháp kỹ thuật trực tiếp cho các nhà sản xuất, lắp ráp ô tô và các trung tâm R&D chuyển đổi động cơ xăng hiện hữu sang sử dụng nhiên liệu sinh học linh hoạt mà không làm giảm tuổi thọ động cơ.
- Tác động kinh tế - xã hội: Hỗ trợ thúc đẩy chuỗi giá trị nông nghiệp bền vững thông qua việc tiêu thụ sắn, mía phế phẩm tại các nhà máy ethanol Bình Phước, Dung Quất, Phú Thọ; giảm phụ thuộc nhập khẩu xăng khoáng, giảm phát thải khí độc hại bảo vệ sức khỏe cộng đồng đô thị.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận nghiên cứu kết hợp 3D-CFD độ trung thực cao và thực nghiệm động cơ chuyên sâu, cùng bộ thông số biên nhiệt động học chuẩn xác.
- Kỹ sư R&D ngành Động lực và Ô tô: Sở hữu cơ sở dữ liệu định lượng và hướng dẫn thiết kế cấu hình vòi phun, định thời góc phun và phân tầng hòa khí tối ưu cho các dòng động cơ thế hệ mới.
- Các nhà hoạch định chính sách: Có đầy đủ căn cứ khoa học thực chứng để tự tin ban hành lộ trình bắt buộc áp dụng xăng E10, E20 trên quy mô toàn quốc, hiện thực hóa các cam kết giảm phát thải quốc tế.
Câu hỏi chuyên sâu
-
Đóng góp lý thuyết độc đáo nhất của luận án là gì?
Đó là việc tích hợp và hiệu chỉnh thành công mô hình ngọn lửa rối Zimont kết hợp mô hình phân rã giọt TAB và bay hơi Dukowicz trong môi trường 3D-CFD cho hòa khí xăng - ethanol không đồng nhất. Nghiên cứu đã làm sáng tỏ hiện tượng phân tầng bay hơi cục bộ và cơ chế làm mát buồng cháy của ethanol, từ đó mở rộng lý thuyết quá trình cháy hai nhiên liệu (Dual-Fuel Combustion Theory) trên động cơ đánh lửa cưỡng bức.
-
Điểm đổi mới phương pháp luận so với các nghiên cứu quốc tế trước đây?
So với Di Iorio và cộng sự (2014) hay Costagliola và cộng sự (2013) vốn chỉ thực nghiệm đơn thuần hoặc mô phỏng 1D động học, luận án đã xây dựng mô hình không gian 3D-CFD hoàn chỉnh mô tả từ quá trình chuyển động tia phun, va chạm màng chất lỏng với vách (Film-Wall), bay hơi đến lan truyền màng lửa rối và phát sinh phát thải. Mô hình đạt độ chính xác cao với sai số áp suất cực đại $p_{\max} < 3{,}5%$ trên toàn dải vận hành.
-
Phát hiện bất ngờ nhất từ dữ liệu nghiên cứu là gì?
Đó là hiện tượng sụt giảm áp suất cháy cực đại và tăng phát thải $\text{CO/HC}$ khi phun trực tiếp (DI) ethanol nguyên chất ở nhiệt độ khí nạp thấp ($< 345\text{ K}$). Nghiên cứu chứng minh rằng mặc dù ethanol có trị số octan vượt trội, nhiệt ẩn hóa hơi quá lớn của nó trở thành rào cản bay hơi nếu không được kết hợp với một lượng xăng mồi hòa trộn trước trên đường nạp (cấu hình GPI-EDI).
-
Luận án có cung cấp quy trình tái lập (Replication Protocol) không?
Có. Luận án mô tả chi tiết toàn bộ điều kiện biên hình học động cơ Daewoo A16DMS, thông số vật lý của các cấu tử nhiên liệu, các bước thiết lập lưới động trong Ansys Fluent, quy trình hiệu chỉnh trên băng thử APA204 và các bước phân tích khí thải bằng KEG-500, cho phép các nhóm nghiên cứu khác tái lập hoàn toàn quy trình.
-
Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
Tập trung phát triển động cơ thông minh ứng dụng công nghệ RCCI (Reactivity Controlled Compression Ignition) và HCCI (Homogeneous Charge Compression Ignition) sử dụng cặp nhiên liệu Diesel - Ethanol/Xăng; tích hợp trí tuệ nhân tạo (AI/Machine Learning) trong việc dự đoán và điều khiển tối ưu quá trình phân tầng nhiên liệu theo thời gian thực.
Kết luận
Luận án "Nghiên cứu tính năng động cơ đánh lửa cưỡng bức sử dụng nhiên liệu xăng - ethanol" của NCS. Nguyễn Quang Trung đã hoàn thành xuất sắc các mục tiêu khoa học và thực tiễn đề ra với 6 kết luận then chốt:
- Xác định giới hạn tối ưu của nhiên liệu hỗn hợp xăng - ethanol hòa trộn sẵn trên động cơ thương mại không thay đổi phần cứng là E10 đến E20, giúp tăng mô-men xoắn lên đến $4{,}47%$, tăng hiệu suất nhiệt có ích lên $29{,}4%$ và giảm mạnh phát thải $\text{CO}$ (tới $66{,}5%$), $\text{HC}$ (tới $42{,}8%$).
- Chứng minh nguyên nhân cốt lõi của sự suy giảm công suất khi tỷ lệ cồn $> \text{E20}$ trên động cơ nguyên bản là do giới hạn thời gian mở vòi phun của ECU gốc dẫn đến hòa khí bị nhạt hóa cục bộ ($f < 1$).
- Xây dựng và kiểm chứng thành công mô hình mô phỏng 3D-CFD đa pha trên nền tảng Ansys Fluent, đạt độ tương thích cao với thực nghiệm (sai số $p_{\max} < 3{,}5%$).
- Khám phá quy luật bay hơi và phân bố hòa khí của ethanol trong xilanh, chỉ ra sự cần thiết của việc gia nhiệt khí nạp ($T_{kn} \ge 345\text{ K}$) hoặc kéo dài thời gian bay hơi khi sử dụng ethanol tỷ lệ cao.
- Đề xuất và chứng minh tính ưu việt vượt trội của giải pháp phun riêng rẽ xăng/ethanol (đặc biệt là cấu hình GPI-EDI), cho phép nâng tỷ lệ ethanol sử dụng lên tới E100, triệt tiêu nguy cơ kích nổ và tối ưu hóa toàn diện quá trình cháy.
- Cung cấp luận cứ khoa học vững chắc khẳng định tính khả thi của việc nâng tỷ lệ phối trộn nhiên liệu sinh học lên E10 và E20 tại Việt Nam, đóng góp trực tiếp vào chiến lược an ninh năng lượng quốc gia và cam kết giảm phát thải toàn cầu.
