Tổng quan nghiên cứu

Khí thải từ các phương tiện giao thông và máy móc nông nghiệp sử dụng động cơ Diesel chiếm khoảng 6,5% các thành phần khí độc hại trực tiếp gây ô nhiễm môi trường và hiệu ứng nhà kính. Nhằm kiểm soát nguồn phát thải này, Thủ tướng Chính phủ đã ban hành Quyết định số 49/2011/QĐ-TTg ngày 01/09/2011 quy định lộ trình áp dụng tiêu chuẩn khí thải mức 4 (Euro 4) từ ngày 01/01/2017 và mức 5 (Euro 5) từ ngày 01/01/2022. Tuy nhiên, phần lớn các dòng động cơ Diesel một xy lanh đang lưu hành tại Việt Nam đều có mức phát thải bồ hóng và oxit nitơ (NOx) vượt quá ngưỡng cho phép, trong khi việc thử nghiệm chế tạo cơ khí trực tiếp đòi hỏi chi phí lớn và thời gian kéo dài.

Đề tài tập trung giải quyết vấn đề tối ưu hóa kết cấu buồng đốt nhằm đồng thời nâng cao công suất và giảm thiểu ô nhiễm môi trường cho động cơ Diesel cỡ nhỏ. Mục tiêu cụ thể là nghiên cứu sự ảnh hưởng của các thông số hình học buồng cháy khoét sâu trên đỉnh piston bao gồm độ sâu, đường kính đáy và đường kính miệng buồng cháy đến công suất cực đại 12,5 kW tại tốc độ 2400 vòng/phút cùng đặc tính phát thải của động cơ Vikyno RV125-2 thông qua phương pháp mô phỏng nhiệt động 3 chiều.

Nghiên cứu được triển khai thực hiện tại Phòng Thí nghiệm Động cơ đốt trong thuộc Trường Đại học Bách khoa – ĐHQG TP.HCM trong khoảng thời gian từ tháng 02/2014 đến tháng 06/2014. Kết quả nghiên cứu mang lại ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp giải pháp thiết kế biên dạng buồng đốt tối ưu giúp tăng hiệu suất chu trình nhiệt, giảm từ 8,5% đến 14,2% lượng phát thải ô nhiễm NOx và tiết kiệm hơn 60% kinh phí chế tạo thử nghiệm khuôn mẫu ban đầu cho các nhà sản xuất máy móc nông nghiệp.

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 lý thuyết hòa khí và động học quá trình cháy trong động cơ Diesel phun trực tiếp (DI Diesel). Quá trình cháy diễn ra qua hai giai đoạn liên tiếp gồm cháy đồng nhất tự phát sau thời gian cháy trễ và cháy khuếch tán. Trong không gian buồng đốt, các dòng xoáy lốc tiếp tuyến và xoáy hướng kính đóng vai trò quyết định đến việc xé nhỏ chùm tia nhiên liệu, bốc hơi màng chất lỏng và hòa trộn không khí.

Mô hình dòng chảy rối Renormalized Group k-epsilon (RNG k-epsilon) được áp dụng nhằm mô tả cấu trúc ngọn lửa quy mô lớn và tính toán chính xác trường nhiệt độ trong không gian xy lanh. Đi kèm với đó là mô hình truyền nhiệt vách (Temperature Wall Function) giúp mô phỏng biến thiên nhiệt động lớp biên giữa khối khí và vách kim loại.

Quá trình tán sương và biến thiên kích thước hạt nhiên liệu được mô tả thông qua mô hình phân rã tia phun lai Kelvin-Helmholtz và Rayleigh-Taylor (KH-RT) với hằng số mô hình B0 xác định bằng 0,6 và chuẩn số Weber giới hạn bằng 6. Sự bay hơi của giọt nhiên liệu tuân theo quy luật bốc hơi Ranz-Marshall với chuẩn số Sherwood xác định tốc độ khuếch tán khối lượng.

Về cơ chế hình thành chất gây ô nhiễm, nồng độ khí thải NOx (chứa 97% đến 98% khí NO) được tính toán theo cơ chế nhiệt Zeldovich mở rộng, vốn phụ thuộc rất mạnh vào nồng độ oxy và nhiệt độ buồng đốt vượt ngưỡng 1800 K (tương đương 1371 độ C). Quá trình sinh ra và oxy hóa hạt bồ hóng được tính toán bằng mô hình động học 8 bước của Foster, mô tả phản ứng nhiệt phân n-tetradecan (C14H30) thành hạt cơ sở C50 và acetylen (C2H2) trước khi trải qua chu trình tích tụ và oxy hóa bề mặt.

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

Nguồn dữ liệu thực nghiệm được thu thập trực tiếp trên băng thử động cơ một xy lanh Vikyno RV125-2 với dung tích 624 cm3, đường kính hành trình piston 90 mm, trang bị hệ thống phun nhiên liệu Bosch CP1 áp suất 1350 bar và vòi phun 4 lỗ tia diện tích 0,2867e-4 cm2. Hệ thống thu thập tín hiệu bao gồm cảm biến áp suất buồng cháy kiểu áp điện AVL-EU21C, cảm biến tốc độ trục khuỷu dạng quang kết hợp encoder chia độ 3600 xung/vòng, và cân nhiên liệu điện tử TE6001 có độ chính xác 0,05 g.

Phương pháp phân tích số trị sử dụng bộ phần mềm KIVA3V-Release kết hợp công cụ tiền xử lý tạo lưới K3PREP và công cụ hậu xử lý Tecplot và Origin. Lưới tính toán 3D được thiết lập đối xứng cho góc 90 độ (tương ứng 1/4 buồng đốt và 1 tia phun) với các vùng logic được phân chia chi tiết theo bản vẽ kỹ thuật piston.

Quy mô khảo sát bao gồm chế độ tốc độ danh định 2400 vòng/phút ở mức 80% tải với cỡ mẫu 10 chu trình nén - cháy liên tiếp để lấy giá trị trung bình ổn định. Phương pháp mô phỏng số 3D CFD được lựa chọn vì cho phép khảo sát độc lập từng biến số hình học buồng đốt trong điều kiện chuẩn hóa nghiêm ngặt tỉ số nén 16,5 và hiệu chỉnh khe hở squish 0,366 cm, điều mà các phương pháp thực nghiệm cơ khí không thể thực hiện nhanh chóng. Quá trình phân tích và đồng nhất dữ liệu mô phỏng với thực nghiệm được hoàn tất trong mốc thời gian từ tháng 03/2014 đến tháng 06/2014.

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

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

Quá trình mô phỏng các phương án buồng cháy khoét sâu trên đỉnh piston với việc thay đổi thông số hình học (độ sâu buồng cháy, đường kính đáy và đường kính miệng) đã mang lại các phát hiện khoa học quan trọng:

  • Tăng đường kính buồng cháy (bao gồm mở rộng đường kính đáy và đường kính miệng) giúp công suất động cơ tăng vượt mức danh định 12,5 kW tại số vòng quay 2400 vòng/phút, đồng thời làm tăng áp suất trung bình chỉ thị trong xy lanh từ 4,2% đến 6,8% so với buồng cháy nguyên bản.
  • Lượng phát thải oxit nitơ (NOx) giảm rõ rệt từ 8,5% đến 14,2% ở dạng buồng cháy mở rộng đường kính miệng và đáy. Khối lượng phát thải NOx đo được tại góc quay trục khuỷu 15 độ sau điểm chết trên (15 deg ATDC) ghi nhận sự suy giảm mật độ các vùng quá nhiệt cục bộ trên 1800 K.
  • Quá trình phát thải hạt bồ hóng diễn ra theo hai giai đoạn đặc trưng: ở giai đoạn cháy ban đầu, lượng muội than hình thành tăng khoảng 12% do tốc độ phân hủy nhiệt gia tăng, nhưng sau đó bước vào giai đoạn cháy kiệt, quá trình oxy hóa bồ hóng diễn ra mạnh mẽ làm lượng bồ hóng tồn đọng giảm sâu, mức phát thải cuối cùng chỉ cao hơn dạng buồng cháy ban đầu khoảng 3% đến 4,8%.
  • Tốc độ tỏa nhiệt cực đại đạt đỉnh cao hơn và dịch chuyển gần về điểm chết trên, giúp suất tiêu hao nhiên liệu riêng ở chế độ tải 80% giảm được khoảng 15 g/kWh so với kết cấu ban đầu của nhà sản xuất.

Thảo luận kết quả

Nguyên nhân cốt lõi của việc cải thiện công suất và giảm phát thải khi tăng đường kính buồng cháy nằm ở cơ chế tương tác giữa tia phun nhiên liệu và dòng khí nạp. Khi mở rộng miệng khoét lõm, tỉ số đường kính xy lanh trên đường kính miệng buồng cháy tạo nên cường độ dòng xoáy hướng kính mạnh mẽ, tỉ lệ thuận với bình phương tỉ số kích thước. Dòng xoáy này đưa oxy từ vùng ngoại vi vào trung tâm ngọn lửa nhanh hơn, giúp quá trình oxy hóa muộn các hạt bồ hóng diễn ra triệt để.

Bên cạnh đó, việc nới rộng không gian đáy buồng cháy kéo dài quãng đường chuyển động tự do của tia nhiên liệu phun ra từ vòi phun 1350 bar, hạn chế sự va đập trực tiếp của nhiên liệu lỏng lên thành piston và ngăn ngừa hiện tượng tích tụ các điểm nhiệt độ cục bộ quá cao. Nhờ đó, nhiệt độ đỉnh ngọn lửa không bị tập trung quá mức, ức chế các phản ứng oxy hóa Nitơ theo chuỗi động học Zeldovich.

Kết quả này hoàn toàn tương thích và bổ trợ cho các nghiên cứu quốc tế của Arturo de Risi tại Đại học Lecce và các mô hình tối ưu hóa của Rolf Reitz tại Trung tâm Nghiên cứu Động cơ thuộc Đại học Wisconsin – Madison. Trong các báo cáo khoa học, toàn bộ diễn biến nhiệt động học được trực quan hóa sinh động thông qua đồ thị biến thiên áp suất xy lanh theo góc quay trục khuỷu (đồ thị P-theta), biểu đồ tốc độ tỏa nhiệt (ROHR), bảng so sánh suất tiêu hao nhiên liệu và các hình ảnh cắt lớp 3D biểu diễn trường nhiệt độ cũng như nồng độ NOx tại thời điểm 20 độ góc quay trục khuỷu sau khi phun.

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

Nhằm chuyển hóa các kết quả nghiên cứu vào thực tế sản xuất và vận hành động cơ Diesel, các giải pháp kỹ thuật và quản lý được đề xuất cụ thể như sau:

  • Cải tiến biên dạng gia công đỉnh piston: Doanh nghiệp sản xuất máy động lực (tiêu biểu như Tổng Công ty Máy động lực và Máy nông nghiệp Việt Nam - VEAM và Công ty Vikyno) cần tiến hành điều chỉnh khuôn đúc và quy trình tiện CNC đỉnh piston động cơ RV125-2 theo hướng tăng đường kính miệng buồng cháy thêm 8% đến 12% và nới rộng đường kính đáy, đặt mục tiêu nâng công suất lên 13 kW và giảm 10% lượng NOx phát thải trong kế hoạch sản phẩm giai đoạn 2025-2027.
  • Tinh chỉnh thông số hệ thống phun nhiên liệu: Bộ phận R&D tại các nhà máy chế tạo động cơ Diesel cần phối hợp việc mở rộng đường kính buồng cháy với việc nâng áp suất kim phun từ 1350 bar lên mức 1500 bar và tối ưu hóa góc tia phun nhằm triệt tiêu hoàn toàn lượng bồ hóng dư thừa trong thời hạn 12 tháng tới.
  • Chuẩn hóa phần mềm mô phỏng trong nghiên cứu động lực học: Các trường đại học kỹ thuật và viện nghiên cứu chuyên ngành cơ khí ô tô máy kéo cần tích hợp quy trình mô phỏng KIVA3V-Release kết hợp thiết bị đo áp điện AVL-EU21C vào chương trình đào tạo sau đại học từ quý 3 năm nay, giúp rút ngắn thời gian thử nghiệm thực tế khoảng 50% cho các dự án nghiên cứu động cơ đốt trong.
  • Xây dựng quy chuẩn kiểm soát khí thải máy nông nghiệp: Cục Đăng kiểm Việt Nam và Bộ Giao thông Vận tải cần ban hành quy định bắt buộc áp dụng tiêu chuẩn tương đương Euro 4 đối với động cơ Diesel máy nông lâm ngư nghiệp một xy lanh, thiết lập lộ trình kiểm định phát thải định kỳ từ giai đoạn 2026-2028.

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

Công trình nghiên cứu cung cấp khối lượng kiến thức chuyên sâu và dữ liệu số hóa có giá trị cao cho nhiều nhóm đối tượng:

  • Kỹ sư thiết kế và chuyên gia R&D ngành động cơ đốt trong: Ứng dụng trực tiếp thông số kích thước hình học tối ưu và quy trình tạo lưới K3PREP để nâng cấp các dòng động cơ Diesel một xy lanh hiện hữu, cắt giảm kinh phí gia công thử nghiệm.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Cơ khí Động lực: Sử dụng toàn bộ khung lý thuyết về mô hình dòng rối RNG k-epsilon, mô hình tán sương KH-RT và cơ chế tạo muội than 8 bước làm tài liệu giảng dạy và nền tảng phát triển đề tài mô phỏng nhiệt động 3D.
  • Doanh nghiệp sản xuất và lắp ráp máy móc nông nghiệp: Khai thác giải pháp cải tiến buồng đốt để gia tăng sức cạnh tranh cho các sản phẩm máy cày, máy gặt đập liên hợp, máy phát điện và ghe thuyền cỡ nhỏ sử dụng động cơ Vikyno.
  • Cán bộ kỹ thuật và chuyên viên quản lý môi trường: Tiếp cận phương pháp luận phân tích phát thải bồ hóng và NOx để xây dựng các giải pháp kỹ thuật giảm thiểu ô nhiễm không khí tại các vùng nông thôn và đô thị.

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

Tại sao nghiên cứu lại chọn đối tượng thực nghiệm là động cơ Vikyno RV125-2?

Động cơ Diesel một xy lanh Vikyno RV125-2 với dung tích 624 cm3 và công suất cực đại 12,5 kW là mẫu động cơ máy nông nghiệp được sử dụng vô cùng rộng rãi tại Việt Nam. Việc nghiên cứu nâng cao tính năng hoạt động và hạ mức khí thải độc hại trên dòng động cơ này mang lại giá trị kinh tế trực tiếp và đáp ứng đúng lộ trình Quyết định số 49/2011/QĐ-TTg của Chính phủ.

Việc thay đổi hình dạng buồng cháy có làm thay đổi tỉ số nén của động cơ không?

Trong suốt quá trình mô phỏng, tỉ số nén của động cơ luôn được giữ cố định ở giá trị nguyên bản là 16,5. Chương trình tính toán tự động bù trừ thể tích khoét lõm với khe hở squish (hiệu chỉnh về mức 0,366 cm) để đảm bảo các thay đổi về công suất và khí thải chỉ phát sinh duy nhất từ yếu tố biên dạng hình học.

Mô hình phân rã tia phun lai KH-RT đóng vai trò gì trong việc nâng cao độ chính xác mô phỏng?

Mô hình kết hợp sóng bề mặt Kelvin-Helmholtz và Rayleigh-Taylor với chuẩn số Weber giới hạn bằng 6 giúp mô tả chính xác quá trình biến đổi từ tia nhiên liệu lỏng áp suất 1350 bar thành các chùm hạt sương kích thước siêu nhỏ, tái hiện trung thực hiện tượng bay hơi và phân bố hòa khí trước khi tự bốc cháy.

Vì sao tăng đường kính buồng cháy lại giúp giảm lượng phát thải khí NOx?

Tăng đường kính đáy và miệng buồng cháy giúp phân tán không gian hòa trộn nhiên liệu đồng đều hơn, tránh hiện tượng cháy quá đậm đặc tại tâm buồng đốt. Điều này làm hạ nhiệt độ ngọn lửa cục bộ xuống dưới ngưỡng kích hoạt nhiệt phân 1800 K, qua đó ức chế phản ứng oxy hóa Nitơ theo cơ chế Zeldovich.

Làm thế nào để bảo đảm kết quả mô phỏng số phản ánh đúng thực tế của động cơ?

Nghiên cứu đã tiến hành đồng nhất nghiêm ngặt đường đặc tính áp suất xy lanh giữa mô phỏng KIVA3V và thực nghiệm băng thử có lắp cảm biến áp điện AVL-EU21C ở tốc độ 2400 vòng/phút và 80% tải. Sai số áp suất cực đại giữa mô phỏng và thực tế được khống chế dưới 3%, bảo đảm độ tin cậy tuyệt đối cho các kết luận.

Kết luận

  • Xây dựng thành công mô hình mô phỏng nhiệt động 3 chiều trên phần mềm KIVA3V-Release có độ tương thích cao với dữ liệu thực nghiệm trên động cơ Vikyno RV125-2 tại số vòng quay 2400 vòng/phút.
  • Chứng minh phương án mở rộng đường kính đáy và miệng buồng cháy giúp nâng cao công suất chỉ thị, giảm suất tiêu hao nhiên liệu khoảng 15 g/kWh và cải thiện hiệu suất nhiệt tổng thể của chu trình.
  • Xác định rõ cơ chế giảm thiểu phát thải NOx từ 8,5% đến 14,2% nhờ việc kiểm soát trường nhiệt độ đỉnh buồng đốt dưới ngưỡng 1800 K.
  • Làm sáng tỏ quy luật hai giai đoạn của quá trình tạo bồ hóng, khẳng định dòng xoáy lốc hướng kính tăng cường đóng vai trò quyết định trong việc oxy hóa muộn muội than.
  • Cung cấp giải pháp thiết kế số hóa chuẩn xác giúp các cơ sở nghiên cứu và sản xuất cơ khí tiết kiệm hơn 60% thời gian và chi phí chế tạo mẫu thử.

Kế hoạch phát triển tiếp theo trong giai đoạn 2025-2026 sẽ tập trung tích hợp hệ thống phun dầu điện tử Common Rail đa điểm áp suất trên 1800 bar và thử nghiệm với nhiên liệu sinh học Biodiesel B20. Hãy tải toàn văn luận văn thạc sĩ kỹ thuật ô tô máy kéo để tham khảo toàn bộ bộ thông số lưới và quy trình mô phỏng buồng đốt chuyên sâu cho các dự án cơ khí động lực của bạn.