Tổng quan nghiên cứu

Động cơ diesel đóng vai trò then chốt trong việc cung cấp nguồn động lực cơ giới hóa nông nghiệp tại Việt Nam, trong đó dòng động cơ Vikyno RV125-2 chiếm hơn 40% tổng sản lượng tiêu thụ của Tổng công ty Máy động lực và Máy nông nghiệp Việt Nam (VEAM). Dù sở hữu độ bền bỉ cao nhờ chuyển giao công nghệ từ Kubota, mẫu động cơ 01 xylanh phun trực tiếp này đang đối mặt với thách thức lớn về khí thải khi kết quả thử nghiệm theo chu trình tiêu chuẩn ISO 8178-C1 tại Phòng Thí nghiệm Trọng điểm Động cơ đốt trong (Đại học Quốc gia TP.HCM) cho thấy hàm lượng phát thải vượt xa ngưỡng quy chuẩn US Tier 2.

Vấn đề cốt lõi đặt ra là sự mâu thuẫn nhiệt động học kinh điển (trade-off) giữa hai chất ô nhiễm chính trong khí xả động cơ diesel: oxit nitơ (NOx) và bồ hóng (soot). Quá trình cháy khuếch tán truyền thống với kim phun đơn tạo ra các vùng cục bộ có nhiệt độ vượt ngưỡng 2000 K làm gia tăng mạnh mẽ NOx, đồng thời các vùng giàu nhiên liệu thiếu oxy lại sản sinh lượng lớn hạt muội rắn. Mục tiêu cụ thể của nghiên cứu là đánh giá định lượng ảnh hưởng của chiến lược phun nhiên liệu nhiều giai đoạn (phun tách hai lần trong một chu trình công tác) nhằm giảm phát thải ô nhiễm mà vẫn duy trì các chỉ tiêu kinh tế - kỹ thuật.

Phạm vi nghiên cứu tập trung vào động cơ Vikyno RV125-2 vận hành ở chế độ 80% phụ tải toàn phần tại tốc độ định mức 2400 vòng/phút, được thực hiện trong mốc thời gian từ tháng 09/2013 đến tháng 06/2014. Ý nghĩa của công trình thể hiện qua việc ứng dụng công cụ mô phỏng 3D hiện đại để tìm ra quy luật cấp nhiên liệu tối ưu, giúp cắt giảm tới 62,27% phát thải NOx, đóng góp cơ sở dữ liệu học thuật quan trọng cho lộ trình nâng chuẩn khí thải các dòng máy nông nghiệp nội địa.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng hệ thống nền tảng lý thuyết nhiệt động học và động học hóa học chuyên sâu cho quá trình cháy hỗn hợp không đồng nhất trong buồng đốt động cơ diesel:

  • Cơ chế hình thành NOx nhiệt Zeldovich mở rộng: Dựa trên chuỗi phản ứng dây chuyền oxy hóa phân tử nitơ trong không khí nạp bởi các nguyên tử oxy hoạt tính ở điều kiện nhiệt độ cao trên 2000 K. Nồng độ NOx phát sinh phụ thuộc hàm mũ vào nhiệt độ cháy cực đại và hệ số dư lượng không khí cục bộ.
  • Mô hình động học muội than 8 bước của Foster và Kazakov: Mô tả sự chuyển hóa phân tử nhiên liệu n-tetradecan (C14H30) qua quá trình nhiệt phân thành phân tử tiền chất Acetylen (C2H2) và các hạt cơ sở (C50), trải qua các giai đoạn hình thành bề mặt, tích tụ hạt và phản ứng oxy hóa cạnh tranh dưới tác động của các gốc tự do O2 và OH.
  • Mô hình dòng chảy rối Renormalized Group k-epsilon (RNG k-ε): Do Han và Reitz phát triển, tối ưu hóa khả năng dự báo cấu trúc ngọn lửa quy mô lớn và truyền nhiệt qua vách buồng đốt chính xác hơn mô hình k-ε tiêu chuẩn.
  • Mô hình phân rã tia phun lai Kelvin-Helmholtz và Rayleigh-Taylor (KH-RT): Giải thích sự phát triển sóng bề mặt không ổn định do chênh lệch vận tốc pha khí - lỏng và lực khí động học, kiểm soát đường kính trung bình Sauter (SMD) và độ xuyên thâu của chùm tia phun.

Các khái niệm then chốt bao gồm: góc phun sớm (SOI), thời gian cháy trễ (ID), góc quay trục khuỷu (°CA), và khoảng nghỉ giữa các lần phun (dwell delay).

Phương pháp nghiên cứu

  • Nguồn dữ liệu và thiết lập thực nghiệm: Nghiên cứu thu thập dữ liệu áp suất chỉ thị xylanh thông qua hệ thống đo đạc chuyên dụng trên băng thử công suất tại Đại học Bách Khoa TP.HCM ở chế độ vận hành 2400 vòng/phút, tải 80%. Dữ liệu thực nghiệm này đóng vai trò chuẩn hóa mô hình tính toán.
  • Quy mô mẫu và trường hợp khảo sát: Phương pháp chọn mẫu tập trung vào 5 phương án phân chia tỷ lệ khối lượng phun lần đầu (90%, 80%, 75%, 50%, 35% tổng lượng cấp nhiên liệu) và 6 mức độ trễ thời gian nghỉ (2°CA, 4°CA, 6°CA, 8°CA, 10°CA, 12°CA). Việc phân bố mẫu này đảm bảo bao quát toàn bộ vùng làm việc chuyển tiếp giữa chế độ phun mồi và phun chính.
  • Lý do lựa chọn phương pháp phân tích: Phần mềm mô phỏng khí động lực học 3 chiều KIVA-3V-ERC được lựa chọn vì khả năng giải đồng thời hệ phương trình bảo toàn khối lượng, động lượng, năng lượng cùng các mô hình hóa học phản ứng chi tiết trong không gian buồng đốt chuyển động thực. Phương pháp này tiết kiệm hơn 75% chi phí chế tạo mẫu thử cơ khí và đo đạc thực nghiệm phức tạp.
  • Quy trình chuẩn hóa: Mô hình hình học buồng đốt piston được chia lưới cấu trúc di động, sau đó tiến hành hiệu chỉnh đồng nhất thông số áp suất nén khi không phun và áp suất cháy khi có phun với sai lệch thực nghiệm dưới 5%.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

  • Cắt giảm vượt bậc phát thải oxit nitơ (NOx): Khi áp dụng chiến lược phun tách theo tỷ lệ 80% cho lần phun đầu và 20% cho lần phun sau với khoảng nghỉ 4°CA, khối lượng NOx phát thải giảm mạnh 62,27% so với chế độ phun đơn tiêu chuẩn.
  • Gia tăng phát thải muội than (bồ hóng): Ở phương án tối ưu NOx (tỷ lệ 80/20, khoảng nghỉ 4°CA), lượng phát thải bồ hóng tăng 53,15% so với phun đơn. Nguyên nhân do lượng nhiên liệu phun lần hai cháy trễ trong giai đoạn giãn nở, nơi nhiệt độ và nồng độ oxy khả dụng đã suy giảm.
  • Ảnh hưởng của tỷ lệ phân chia lượng phun: Khi hạ thấp tỷ lệ phun ban đầu xuống mức 50% hoặc 35%, áp suất cháy cực đại trong xylanh giảm rõ rệt từ 64,8 bar xuống dưới 54,2 bar, kéo theo sự sụt giảm nhiệt độ trung bình cực đại hơn 220 K, khiến công suất chỉ thị của động cơ giảm khoảng 3,8% đến 5,2%.
  • Ảnh hưởng của thời gian nghỉ giữa hai lần phun: Khi kéo dài khoảng nghỉ từ 2°CA lên 12°CA với cùng tỷ lệ 80/20, phát thải NOx giảm dần từ mức giảm 48,1% xuống mức giảm sâu 68,4%, tuy nhiên lượng muội than có xu hướng tăng nhanh khi khoảng nghỉ vượt quá 8°CA do không đủ thời gian để quá trình oxy hóa C2H2 và hạt C50 hoàn tất.

Thảo luận kết quả

Cơ chế vật lý đằng sau sự sụt giảm NOx nằm ở việc kiểm soát lượng nhiệt giải phóng trong giai đoạn cháy sớm. Ở chế độ phun đơn, lượng lớn hòa khí tự bốc cháy đồng loạt tạo ra xung áp suất và đỉnh nhiệt độ vượt 2200 K. Ngược lại, việc chia tách 80% nhiên liệu ở pha đầu giúp khống chế đỉnh nhiệt độ cháy cục bộ tại thời điểm 20°ATDC xuống dưới 1950 K, ức chế trực tiếp tốc độ phản ứng oxy hóa nitơ theo cơ chế Zeldovich.

Đối với bồ hóng, dữ liệu nghiên cứu cho thấy sự cạnh tranh khốc liệt giữa tốc độ tạo tiền chất Acetylen và tốc độ oxy hóa bề mặt hạt rắn. Khi nhiên liệu được cấp vào giai đoạn muộn (phun lần hai), thể tích xylanh giãn nở làm loãng mật độ oxy, cản trở phản ứng oxy hóa muội than.

Các dữ liệu nghiên cứu này được trực quan hóa sinh động qua hệ thống đồ thị diễn biến áp suất chỉ thị theo góc quay trục khuỷu từ -20°CA đến +60°CA, đồ thị trường phân bố nhiệt độ cục bộ 3D tại 20°ATDC, và bảng đối chiếu đa biến giữa nồng độ phát thải (tính bằng g/chu trình) với các mức thời gian nghỉ khác nhau. Kết quả này hoàn toàn tương thích với các nghiên cứu thực nghiệm quốc tế về phun phân đoạn trên động cơ diesel cỡ nhỏ.

Đề xuất và khuyến nghị

  • Nâng cấp hệ thống phun nhiên liệu điện tử Common Rail: Triển khai thay thế bơm cao áp cơ khí truyền thống bằng hệ thống kim phun điện tử Common Rail (CRDI) áp suất phun trên 1200 bar, cho phép kiểm soát linh hoạt thời điểm phun và định lượng chính xác tỷ lệ phun 80/20 cùng khoảng cách 4°CA. Thời gian thực hiện mục tiêu trong vòng 12 đến 18 tháng, chủ thể thực hiện là đội ngũ R&D Tổng công ty VEAM phối hợp cùng các viện nghiên cứu.
  • Tối ưu hóa thiết kế buồng cháy đỉnh piston: Cải tiến biên dạng hốc buồng đốt piston từ dạng trụ tròn mở sang biên dạng lượn sâu dạng re-entrant bowl nhằm tăng cường tỷ số xoáy lốc (swirl ratio) lên mức trên 2,5. Giải pháp này giúp gia tăng cường độ hòa trộn khí - nhiên liệu cho lần phun thứ hai, mục tiêu nâng hiệu suất oxy hóa muội than thêm 25-30% trong giai đoạn 2025-2026 do các kỹ sư thiết kế cơ khí động lực đảm nhiệm.
  • Tích hợp hệ thống tuần hoàn khí xả làm mát (Cooled EGR): Ứng dụng tỷ lệ luân hồi khí xả từ 10% đến 15% kết hợp cùng chiến lược phun tách 80/20 nhằm kéo giảm thêm 15-20% nồng độ NOx mà không làm suy giảm hiệu suất nhiệt của động cơ.
  • Hoàn thiện chính sách quy chuẩn khí thải máy nông nghiệp: Kiến nghị Bộ Giao thông Vận tải phối hợp cùng Bộ Nông nghiệp và Phát triển Nông thôn ban hành lộ trình áp dụng bắt buộc tiêu chuẩn khí thải mức Euro 4 hoặc Tier 4 cho toàn bộ dòng động cơ tĩnh tại và máy cơ giới nông nghiệp từ năm 2026 đến năm 2030.

Đối tượng nên tham khảo luận văn

  • Kỹ sư R&D tại các nhà máy chế tạo động cơ đốt trong: Các chuyên gia kỹ thuật tại SVEAM, Vinappro và VEAM có thể trực tiếp ứng dụng bộ thông số tối ưu (tỷ lệ 80/20, khoảng nghỉ 4°CA) để cải tiến hệ thống cung cấp nhiên liệu trên các dòng động cơ diesel 1 xylanh hiện hành.
  • Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Động lực: Luận văn là tài liệu tham khảo mẫu mực về quy trình thiết lập lưới tính toán động 3D, cài đặt các điều kiện biên và tích hợp các mô hình hóa học phản ứng phức tạp trong phần mềm KIVA-3V.
  • Các nhà hoạch định chính sách môi trường và cơ quan kiểm định: Cung cấp cơ sở dữ liệu khoa học định lượng về đặc tính phát thải thực tế của dòng động cơ nông nghiệp nông thôn, phục vụ việc xây dựng lộ trình quy chuẩn kỹ thuật quốc gia.
  • Doanh nghiệp phát triển thiết bị xử lý khí thải sau buồng đốt: Cung cấp thông số phát thải đầu vào chính xác (nồng độ NOx và kích thước hạt bồ hóng) để định cỡ và thiết kế tối ưu các bộ lọc muội DPF hoặc bộ xúc tác oxy hóa DOC.

Câu hỏi thường gặp

  • Tại sao chiến lược phun tách 80/20 lại giảm được hơn 62% lượng NOx? Việc tách nhiên liệu thành hai giai đoạn giúp phân tán lượng nhiệt giải phóng trong xylanh. Đỉnh nhiệt độ cháy cục bộ tại thời điểm 20°ATDC được khống chế dưới 1950 K, ngăn cản các phản ứng nhiệt phân oxy và nitơ theo cơ chế Zeldovich vốn chỉ hoạt động mạnh ở nhiệt độ trên 2000 K.
  • Nguyên nhân chính khiến lượng phát thải bồ hóng gia tăng 53,15% là gì? Lượng nhiên liệu 20% phun ở giai đoạn sau đi vào buồng đốt khi thể tích xylanh đang giãn nở nhanh. Nhiệt độ môi trường cháy giảm xuống cùng với nồng độ oxy khả dụng thấp khiến các chuỗi hydrocacbon thơm không kịp oxy hóa hoàn toàn, dẫn đến sự tồn đọng muội than rắn.
  • Độ chính xác của phần mềm mô phỏng KIVA-3V được kiểm chứng như thế nào? Mô hình mô phỏng được chuẩn hóa trực tiếp với đường đặc tính áp suất xylanh thu thập từ băng thử thực nghiệm tại 2400 vòng/phút, 80% tải. Sai số tương đối về áp suất cháy cực đại giữa mô phỏng và thực tế được kiểm soát chặt chẽ ở mức dưới 5%.
  • Khoảng nghỉ giữa hai lần phun tối ưu nhất trên động cơ Vikyno RV125-2 là bao nhiêu? Khoảng nghỉ tối ưu xác định được là 4°CA. Mức độ trễ này đủ lớn để tách biệt hai đỉnh giải phóng nhiệt giúp hạ nhiệt độ cháy nhằm triệt tiêu NOx, nhưng không quá trễ để hạn chế sự tích tụ quá mức của các tiền chất muội than C2H2.
  • Phương pháp phun tách có gây suy giảm công suất danh định của động cơ không? Ở phương án tối ưu 80% lượng phun đầu và khoảng nghỉ 4°CA, công suất chỉ thị của động cơ chỉ suy giảm không đáng kể khoảng 1,2% đến 1,8% so với phun đơn, trong khi các chỉ tiêu kinh tế nhiên liệu vẫn được đảm bảo ổn định.

Kết luận

  • Luận văn đã mô hình hóa thành công và toàn diện quá trình cháy cùng cơ chế phát thải ô nhiễm trên động cơ diesel Vikyno RV125-2 bằng công cụ CFD 3D KIVA-3V.
  • Xác định phương án tối ưu hóa quá trình cấp nhiên liệu với tỷ lệ phun ban đầu 80%, lượng phun sau 20% và khoảng nghỉ 4°CA, giúp cắt giảm 62,27% phát thải NOx.
  • Làm sáng tỏ quy luật đánh đổi nhiệt động học giữa NOx và bồ hóng (tăng 53,15%) khi thay đổi linh hoạt các thông số phân chia tia phun.
  • Kiểm chứng độ tin cậy cao của phương pháp mô phỏng với sai số áp suất xylanh dưới 5% so với thực nghiệm đo đạc trên băng thử động lực.
  • Cung cấp cơ sở khoa học vững chắc phục vụ việc nâng cấp công nghệ phun điện tử cho dòng máy nông nghiệp chiếm hơn 40% thị phần Việt Nam.

Kế hoạch phát triển tiếp theo cần tập trung chế tạo mẫu kim phun điều khiển điện tử Common Rail để thử nghiệm thực tế trên băng thử công suất trong 12 tháng tới. Các viện nghiên cứu và nhà sản xuất động cơ hãy chủ động ứng dụng dữ liệu nghiên cứu này vào quy trình nâng cấp sản phẩm nhằm kiến tạo thế hệ động cơ nông nghiệp sạch, tiết kiệm nhiên liệu và thân thiện với môi trường.