Tổng quan nghiên cứu

Khí thải từ các phương tiện giao thông vận tải sử dụng động cơ đốt trong đang là một trong những nguyên nhân hàng đầu gây suy thoái chất lượng không khí toàn cầu. Khí xả động cơ Diesel thông thường chứa từ 50 đến 1000 ppm oxit nitơ (NOx), 10 đến 500 ppm monoxit carbon (CO) và đặc biệt là mật độ bồ hóng cao với các hạt sơ cấp kích thước siêu mịn từ 10 đến 80 nm. Vấn đề cốt lõi trong kỹ thuật động cơ là sự tồn tại của mối quan hệ nghịch biến kinh điển: các giải pháp làm giảm phát thải bồ hóng thường kéo theo sự gia tăng nhiệt độ cháy cục bộ, làm bùng phát nồng độ NOx độc hại và ngược lại.

Mục tiêu trọng tâm của đề tài là giải quyết nghịch lý trên bằng cách nghiên cứu định lượng ảnh hưởng của tỷ lệ xoáy lốc dòng khí nạp và kỹ thuật phun nhiên liệu nhiều lần (phun tách) đến hiệu suất công tác cùng đặc tính phát thải. Phạm vi nghiên cứu tập trung vào không gian buồng đốt của động cơ nghiên cứu chuyên dụng 1 xy lanh AVL 5402 có công suất tối đa 8 kW, tích hợp hệ thống phun nhiên liệu tích áp Common Rail hiện đại. Công trình được thực hiện thông qua mô phỏng nhiệt động và động lực học chất lưu 3D đa chiều trên nền tảng phần mềm chuyên dụng KIVA-3V. Kết quả đạt được mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở dữ liệu xác thực giúp các kỹ sư tối ưu hóa quá trình tạo hòa khí, cắt giảm từ 20% đến 35% lượng phát thải bồ hóng, đồng thời khống chế sự gia tăng của NOx mà vẫn bảo toàn mức công suất định mức 8 kW của động cơ.

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

Khung lý thuyết áp dụng

Quá trình hình thành hỗn hợp cháy và phát sinh ô nhiễm trong xy lanh động cơ Diesel được phân tích dựa trên sự kết hợp chặt chẽ của các mô hình lý thuyết nhiệt - thủy khí động lực học tiên tiến:

  • Mô hình dòng chảy rối Renormalized Group k-epsilon (RNG k-e): Mô tả chính xác chuyển động rối của dòng khí xoáy lốc trong không gian buồng cháy phức tạp.
  • Mô hình phân rã tia phun lai Kelvin-Helmholtz và Rayleigh-Taylor (KH-RT): Mô phỏng cơ chế phân rã tia nhiên liệu lỏng ở áp suất cao từ 200 đến 1700 atm thành các chùm hạt cực nhỏ dưới tương tác khí động học.
  • Mô hình bốc hơi hạt và phân bố kích thước hạt Sauter Mean Diameter (DSM) theo công thức bán thực nghiệm Hiroyasu - Kadota.
  • Mô hình hình thành và oxy hóa bồ hóng 8 bước của Foster: Mô tả chi tiết các giai đoạn từ hình thành hạt cơ sở (precursor), phát triển bề mặt, đông tụ thành chuỗi đến quá trình oxy hóa hạt rắn.
  • Cơ chế Zeldovich mở rộng: Xác định tốc độ tạo thành oxit nitơ nhiệt (thermal NOx) trong các vùng cục bộ có nhiệt độ vượt ngưỡng 2000 K và giàu oxy.

Các khái niệm then chốt bao gồm: tỷ lệ xoáy lốc (Swirl Ratio - tỷ số giữa vận tốc góc của dòng khí nạp và vận tốc góc của trục khuỷu), độ xuyên thâu của tia phun, quá trình cháy khuếch tán và thời gian cháy trễ.

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

Nghiên cứu ứng dụng phương pháp mô phỏng số động lực học dòng chảy đa chiều (CFD) bằng mã nguồn KIVA-3V của Phòng thí nghiệm Quốc gia Los Alamos. Nguồn dữ liệu hình học và điều kiện vận hành được lấy chuẩn hóa từ mẫu động cơ đơn xy lanh AVL 5402.

Do tính chất đối xứng tâm của vòi phun đa lỗ, không gian buồng cháy đỉnh piston được phân chia thành mô hình lưới quét cung tròn 72 độ với hàng chục nghìn phần tử khối. Việc lựa chọn phương pháp sai phân và thể tích hữu hạn trong KIVA-3V là tối ưu vì cho phép giải đồng thời hệ phương trình Navier-Stokes phi tuyến, phương trình bảo toàn năng lượng và truyền chất nhiều pha mà phương pháp thực nghiệm khó đo đạc chi tiết từng mili-giây.

Ma trận khảo sát bao gồm 5 trường hợp phân bố lượng phun:

  • Phun đơn truyền thống (100% nhiên liệu).
  • 4 trường hợp phun tách 2 lần với các tỷ lệ khối lượng: 25-3-75, 50-3-50, 75-3-25 và 90-3-10 (trong đó số 3 biểu thị khoảng cách ngắt quãng 3 độ góc quay trục khuỷu giữa hai lần phun).
  • Các dải tỷ lệ xoáy lốc thay đổi từ giá trị cơ sở 1.0 lên các giá trị cao hơn.

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

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

Quá trình trích xuất và xử lý dữ liệu từ mô phỏng CFD 3D đã làm sáng tỏ 4 phát hiện quan trọng:

  • Thứ nhất, tăng tỷ lệ xoáy lốc dòng nạp từ mức 1.0 lên các giá trị cao hơn giúp tăng cường động năng xé tơi tia nhiên liệu, gia tăng nhanh chóng chiều rộng góc mở tia và thúc đẩy quá trình oxy hóa muội than, giúp tổng khối lượng bồ hóng phát thải giảm từ 20% đến 35%.
  • Thứ hai, khi chuyển từ chế độ phun đơn sang chế độ phun tách với tỷ lệ 25-3-75, lượng phát thải NOx giảm rõ rệt khoảng 15% đến 20% ở giai đoạn cháy ban đầu so với trường hợp phun 1 lần, do giảm bớt lượng nhiên liệu tham gia vào giai đoạn cháy nhanh tự bốc cháy.
  • Thứ ba, ở các tỷ lệ phun tách tiếp theo gồm 50-3-50, 75-3-25 và 90-3-10, lượng bồ hóng tiếp tục giữ xu hướng giảm nhưng nồng độ NOx lại có chiều hướng gia tăng do lượng nhiên liệu phun đợt 2 bốc cháy trong điều kiện nền nhiệt độ xy lanh đã bị nung nóng ở mức cao.
  • Thứ tư, đường đặc tính công suất đầu ra và áp suất chỉ thị trung bình duy trì ổn định xung quanh ngưỡng 7.5 đến 8.0 kW, suất tiêu hao nhiên liệu không bị suy giảm tiêu cực, khẳng định tính khả thi của các phương án hòa trộn cải tiến.

Thảo luận kết quả

Nguyên nhân cốt lõi của sự suy giảm bồ hóng khi tăng xoáy lốc là do luồng không khí chuyển động xoắn cuốn theo các hạt nhiên liệu bay hơi, mở rộng diện tích tiếp xúc với oxy tự do trong buồng đốt. Nhờ đó, các hạt bồ hóng sơ cấp có đường kính 15 đến 30 nm liên tục bị oxy hóa trước khi kịp đông tụ thành chuỗi hạt rắn lớn có kích thước 100 đến 150 nm.

Tuy nhiên, xoáy lốc quá mạnh cũng làm tăng cường độ truyền nhiệt qua vách xy lanh và đẩy nhiệt độ cục bộ tại vùng màng lửa lên trên 2000 K, kích hoạt mạnh mẽ phản ứng giữa N2 và O2 tạo thành NOx. Kết quả này hoàn toàn tương đồng với các công bố khoa học của Fuchs và Rutland hay nghiên cứu của nhóm tác giả Y. Zha và N. Ladommatos về sự tồn tại của một ngưỡng tỷ lệ xoáy lốc tối ưu.

Toàn bộ các diễn biến nhiệt - động học trên được phản ánh trực quan qua đồ thị biểu diễn áp suất trung bình và nhiệt độ cực đại theo góc quay trục khuỷu (từ trước điểm chết trên BTDC đến sau điểm chết trên ATDC), kết hợp bảng đối chiếu nồng độ phát thải bồ hóng và NOx tương ứng với từng tỷ lệ cấp nhiên liệu.

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

Nhằm hiện thực 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, 4 giải pháp trọng tâm được đề xuất như sau:

  • Thiết lập bản đồ điều khiển phun phân tách đa tầng thích ứng (Adaptive Multi-Injection Mapping) trên bộ điều khiển động cơ điện tử (ECU), cài đặt tỷ lệ phun mồi 20% đến 25% ở dải tải trung bình nhằm giảm 15% đến 20% nồng độ NOx và hạ 25% lượng bồ hóng; thời gian hoàn thiện thuật toán trong vòng 6 đến 9 tháng bởi các kỹ sư điều khiển hệ thống động lực.
  • Tái thiết kế biên dạng cổng nạp khí xoắn ốc (Helical Port) để cố định tỷ lệ xoáy lốc tối ưu trong khoảng 1.2 đến 1.5, giúp tăng tốc độ hòa trộn không khí - nhiên liệu thêm 12%; lộ trình thực hiện 12 tháng do các phòng R&D cơ khí chế tạo động cơ đảm nhiệm.
  • Tích hợp đồng bộ kỹ thuật phun tách 25-3-75 với hệ thống hồi lưu khí xả làm mát (Cooled EGR), hướng tới mục tiêu kiểm soát nồng độ NOx trong khí xả dưới ngưỡng 200 ppm mà không làm tăng quá 10% lượng bồ hóng; triển khai thử nghiệm thực tế trong 1 đến 2 năm tại các trung tâm thử nghiệm khí thải.
  • Ứng dụng quy trình mô phỏng số KIVA-3V vào giai đoạn thiết kế tiền sản xuất tại các doanh nghiệp ô tô, giúp cắt giảm 40% chi phí gia công mẫu thử vật lý và rút ngắn thời gian thử nghiệm thực nghiệm từ 18 tháng xuống còn 6 tháng.

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

Công trình nghiên cứu mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chuyên môn:

  • Kỹ sư nghiên cứu và phát triển (R&D) ô tô: Khai thác các thông số tối ưu về áp suất phun Common Rail (200 đến 1700 atm), góc mở tia phun và dạng buồng cháy chỏm cầu để hoàn thiện thiết kế động cơ Diesel thế hệ mới đạt tiêu chuẩn khí thải cao.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Động lực: Sử dụng toàn bộ phương pháp luận thiết lập lưới buồng đốt 72 độ, điều kiện biên và hệ phương trình vi phân mô phỏng CFD 3D (RNG k-e, KH-RT, Foster 8 bước) làm khung tham chiếu chuẩn cho các đề tài chuyên sâu.
  • Cán bộ quản lý kỹ thuật và hoạch định tiêu chuẩn môi trường: Tham khảo các số liệu định lượng về đặc tính phát thải CO, HC, SO2, bồ hóng và NOx để xây dựng quy chuẩn kỹ thuật quốc gia về kiểm soát khí thải phương tiện giao thông cơ giới.
  • Giảng viên các trường đại học kỹ thuật: Ứng dụng mô hình tính toán và hình ảnh trường nhiệt độ, trường nồng độ hạt trong buồng đốt làm học liệu trực quan cho các môn học Động cơ đốt trong, Lý thuyết quá trình cháy và Ô nhiễm môi trường.

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

Tại sao trong động cơ Diesel luôn xuất hiện sự mâu thuẫn giữa phát thải NOx và bồ hóng?
Bồ hóng sinh ra tại các vùng hòa khí đậm đặc thiếu oxy, còn NOx hình thành ở vùng hòa khí có nhiệt độ cao vượt 2000 K có thừa oxy. Khi tăng xoáy lốc hoặc tăng góc phun sớm để đốt cháy triệt để bồ hóng, nhiệt độ ngọn lửa tăng cao lại kích hoạt phản ứng nhiệt Zeldovich khiến nồng độ NOx bùng phát.

Chiến lược phun tách 25-3-75 mang lại lợi ích gì vượt trội so với phun đơn?
Chiến lược 25-3-75 phun trước 25% lượng nhiên liệu đóng vai trò mồi lửa, sau đó ngừng 3 độ góc quay trục khuỷu rồi mới phun 75% còn lại. Cơ chế này kiểm soát tốc độ tăng áp suất, hạn chế nhiệt độ cực đại trong buồng đốt, từ đó cắt giảm 15% đến 20% lượng NOx ban đầu mà vẫn giữ phát thải bồ hóng ở mức thấp.

Có nên tăng tỷ lệ xoáy lốc lên mức tối đa để triệt tiêu hoàn toàn muội than không?
Không nên tăng xoáy lốc quá mức giới hạn. Khi tỷ lệ xoáy lốc quá cao, độ xuyên thâu của tia nhiên liệu bị suy giảm khiến hòa khí không phân bố đều ra sát thành xy lanh. Đồng thời, tổn thất nhiệt truyền qua vách tăng mạnh làm giảm hiệu suất nhiệt và gây cháy không hoàn toàn ở một số vùng.

Tại sao phần mềm mô phỏng KIVA-3V lại được ưa chuộng trong nghiên cứu quá trình cháy động cơ?
KIVA-3V giải quyết xuất sắc các bài toán động lực học dòng chảy rối đa chiều, quá trình phân rã tia phun, bốc hơi và động học hóa học phức tạp. Sử dụng KIVA-3V giúp tiết kiệm từ 40% đến 60% chi phí chế tạo mẫu thử thực nghiệm và cho phép quan sát chi tiết diễn biến bên trong xy lanh.

Đặc tính kỹ thuật cơ bản của động cơ nghiên cứu AVL 5402 là gì?
AVL 5402 là động cơ Diesel nghiên cứu tiêu chuẩn 1 xy lanh, sử dụng hệ thống phun nhiên liệu tích áp Common Rail phun trực tiếp với công suất phát tối đa 8 kW. Kết cấu buồng cháy dạng chỏm cầu trên piston rất thuận lợi cho việc mô hình hóa lưới CFD sector 72 độ đối xứng.

Kết luận

  • Đã thiết lập thành công mô hình mô phỏng CFD 3D hoàn chỉnh trên phần mềm KIVA-3V để khảo sát động lực học quá trình cháy trên động cơ AVL 5402.
  • Chứng minh việc tăng tỷ lệ xoáy lốc khí nạp giúp gia tăng hòa trộn oxy, cắt giảm hiệu quả từ 20% đến 35% lượng phát thải hạt bồ hóng.
  • Xác định chiến lược phun tách 25-3-75 là giải pháp tối ưu giúp hạ thấp 15% đến 20% nồng độ phát thải NOx ở giai đoạn đầu mà không làm suy giảm công suất danh định 8 kW.
  • Cung cấp luận cứ khoa học giải quyết mối quan hệ đánh đổi phức tạp giữa bồ hóng và NOx trong động cơ Diesel phun trực tiếp.
  • Đặt nền móng vững chắc cho việc ứng dụng công nghệ điều khiển phun điện tử Common Rail đa giai đoạn kết hợp hệ thống hồi lưu khí xả EGR trong thực tế.

Trong giai đoạn 12 tháng tiếp theo, nhóm nghiên cứu khuyến nghị tiếp tục mở rộng mô phỏng kết hợp hệ thống EGR làm mát và tiến hành thực nghiệm đối chứng trên băng thử động cơ chuyên dụng. Hãy tham khảo chi tiết toàn văn luận văn để ứng dụng ngay các giải pháp kỹ thuật buồng đốt tiên tiến này vào công trình của bạn!