Tổng quan nghiên cứu

Thị trường ô tô sử dụng động cơ Diesel trên toàn cầu đã chứng kiến sự tăng trưởng vượt bậc trong nhiều thập kỷ qua. Theo số liệu từ Hiệp hội các nhà sản xuất Ô tô Châu Âu (ACEA), tỷ lệ ô tô con trang bị động cơ Diesel tại khu vực Tây Âu đã tăng từ mức 13,8% vào năm 1990 lên đến 53,3% vào năm 2007, thậm chí đạt mức trên 77% tại các quốc gia như Bỉ và Luxembourg. Tại khu vực Châu Á và các thị trường mới nổi, nhu cầu tiêu thụ động cơ Diesel duy trì tốc độ tăng trưởng bình quân khoảng 3% mỗi năm. Ưu thế vượt trội về hiệu suất nhiệt và tính kinh tế nhiên liệu đã giúp dòng động cơ này giữ vai trò then chốt trong vận tải hàng hóa và hành khách.

Tuy nhiên, lượng khí thải độc hại chiếm khoảng 6,5% trong tổng thành phần khí xả của động cơ Diesel, đặc biệt là bồ hóng (soot/particulate matter) và oxít nitơ (NOx), đang là rào cản lớn nhất đối với các quy chuẩn môi trường ngày càng khắt khe. Vấn đề cốt lõi nằm ở mối quan hệ mâu thuẫn đánh đổi kinh điển: các biện pháp làm giảm nhiệt độ để hạ nồng độ NOx thường dẫn đến hiện tượng cháy không hoàn toàn, làm gia tăng muội than và ngược lại.

Nghiên cứu này được thực hiện nhằm giải quyết mục tiêu tối ưu hóa đồng thời công suất động cơ và giảm thiểu phát thải ô nhiễm môi trường thông qua việc điều chỉnh các thông số phun nhiên liệu, cụ thể là góc phun sớm (Start of Injection - SOI) và thời gian phun (Duration of Injection - DOI). Nghiên cứu được triển khai trong phạm vi 12 tháng từ tháng 7 năm 2010 đến tháng 7 năm 2011, ứng dụng trên đối tượng động cơ thử nghiệm 1 xy lanh của hãng AVL tại Phòng Thí nghiệm Động cơ đốt trong thuộc Trường Đại học Bách khoa - Đại học Quốc gia TP. Hồ Chí Minh. Việc ứng dụng giải pháp mô phỏng số 3D giúp rút ngắn 60% thời gian thử nghiệm thực tế và giảm đáng kể chi phí gia công mẫu vật lý trong quá trình phát triển động cơ.

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

Khung lý thuyết áp dụng

Nền tảng lý thuyết của đề tài tập trung vào động học dòng chảy rối, quá trình tán sương nhiên liệu lỏng và các cơ chế phản ứng nhiệt hóa học trong buồng cháy động cơ Diesel phun trực tiếp. Các mô hình lý thuyết chủ đạo bao gồm:

  • Mô hình dòng chảy rối Renormalized Group k-epsilon (RNG k-e): Được sử dụng để mô tả cấu trúc dòng khí chuyển động xoáy lốc và tương tác nhiệt động lực học trong không gian buồng đốt kín với độ chính xác cao hơn so với mô hình k-epsilon tiêu chuẩn.
  • Mô hình phân rã tia phun lai Kelvin-Helmholtz và Rayleigh-Taylor (KH-RT): Mô tả chính xác hai giai đoạn phân rã của tia nhiên liệu lỏng áp suất cao thoát ra từ vòi phun. Giai đoạn đầu do lực xé khí động học bề mặt tạo ra các hạt sơ cấp, giai đoạn sau do gia tốc làm vỡ hạt thành các giọt nhiên liệu cực nhỏ với đường kính trung bình Sauter (DSM) đạt mức từ 10 đến 50 micromet.
  • Mô hình hóa học tạo bồ hóng 8 bước của Foster và cơ chế Zeldovich mở rộng: Mô hình hóa quá trình hình thành hạt mầm carbon, sự phát triển bề mặt, đông tụ và quá trình oxy hóa bồ hóng; kết hợp cùng cơ chế Zeldovich để dự báo nồng độ NOx hình thành chủ yếu ở nhiệt độ cục bộ vượt ngưỡng 1371 độ C (tương đương 2500 độ F).
  • Các khái niệm nhiệt động cốt lõi: Bao gồm thời kỳ cháy trễ (Ignition Delay), tốc độ tỏa nhiệt thực (Apparent Heat Release Rate), giai đoạn cháy hỗn hợp hòa trộn trước và giai đoạn cháy khuếch tán.

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

Nghiên cứu kết hợp giữa phương pháp thực nghiệm thu thập dữ liệu và phương pháp mô phỏng số động lực học dòng chảy (CFD) bằng phần mềm KIVA-3V ERC Release 2.

Cỡ mẫu nghiên cứu được thiết lập dựa trên lưới tính toán 3 chiều mô phỏng hình học buồng đốt piston AVL dạng chỏm cầu với cung góc quay 72 độ. Việc lựa chọn góc 72 độ xuất phát từ phương pháp lấy mẫu đối xứng chu kỳ cho vòi phun 5 lỗ (360 độ chia 5 bằng 72 độ), giúp tối ưu hóa tài nguyên phần cứng và giảm 80% thời gian tính toán mà vẫn đảm bảo tính hội tụ vật lý của dòng môi chất.

Phương pháp phân tích dựa trên sự đồng nhất hóa dữ liệu thực nghiệm thu thập từ động cơ AVL 5402 Single Cylinder Research Engine (áp suất có ích trung bình 620 kPa, dải tốc độ làm việc 1000 vòng/phút, tỷ số nén 16:1). Quá trình hiệu chỉnh bao gồm việc khớp nối đường đặc tính áp suất xy lanh khi không phun nhiên liệu và khi có phun nhiên liệu ở chế độ 80% tải nhằm đảm bảo sai số giữa mô phỏng và thực tế nằm dưới 5%. Timeline nghiên cứu hoàn thiện toàn bộ chuỗi mô phỏng trong vòng 365 ngày.

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 trên phần mềm KIVA-3V đã mang lại những số liệu định lượng chi tiết về sự thay đổi của áp suất, nhiệt độ, công suất và phát thải theo các thông số phun:

  • Tác động của góc phun sớm (SOI): Khảo sát dải thời điểm phun từ 10 đến 22 độ trước điểm chết trên (degBTDC). Khi tăng góc phun sớm từ 12 lên 22 degBTDC, thời điểm đạt nhiệt độ lớn nhất trong buồng cháy dịch chuyển từ 8 độ sau điểm chết trên (degATDC) về đúng điểm chết trên (0 degTDC). Áp suất cực đại trong xy lanh tăng hơn 25%, kéo theo lượng NOx phát thải tăng vọt hơn 45% do nhiệt độ đỉnh buồng đốt vượt quá 2200 K. Ngược lại, lượng muội than bồ hóng giảm từ 35% đến 40% do thời gian hòa trộn nhiên liệu với không khí được kéo dài.
  • Tác động của thời gian phun (DOI): Khảo sát dải thời gian phun từ 4 CAD đến 18 CAD (độ quay trục khuỷu). Kéo dài thời gian phun (đồng nghĩa với việc giảm áp suất phun khi giữ nguyên lượng nhiên liệu chu trình) làm giảm tốc độ tỏa nhiệt cực đại, khiến nhiệt độ trung bình buồng cháy giảm và lượng phát thải NOx giảm mạnh tới hơn 50%.
  • Phát hiện về cơ chế hai giai đoạn của bồ hóng: Khi tăng DOI từ 4 CAD lên 12 CAD, lượng bồ hóng hình thành tăng mạnh do tia nhiên liệu có áp suất thấp, đường kính hạt Sauter lớn làm giảm tốc độ bay hơi. Tuy nhiên, khi tiếp tục tăng DOI từ 12 CAD lên 18 CAD, lượng bồ hóng phát thải lại có xu hướng giảm nhẹ do lượng nhiên liệu cháy muộn trong kỳ giãn nở được oxy hóa bổ sung dưới tác động của dòng xoáy lốc tồn dư.

Thảo luận kết quả

Cơ chế tạo phát thải trong động cơ Diesel phun trực tiếp phụ thuộc chặt chẽ vào nhiệt độ ngọn lửa cục bộ và hệ số dư lượng không khí. Khi phun nhiên liệu quá sớm (20 đến 22 degBTDC), phần lớn nhiên liệu bốc hơi và hòa trộn trước khi tự bén lửa, tạo ra đỉnh tỏa nhiệt ban đầu rất dốc. Điều này giải thích tại sao nồng độ NOx tăng đột biến theo cơ chế nhiệt Zeldovich. Ngược lại, góc phun muộn làm dịch chuyển pha cháy vào kỳ giãn nở, áp suất và nhiệt độ giảm giúp ức chế NOx nhưng làm tăng lượng bồ hóng do thiếu nhiệt độ để oxy hóa muộn các hạt cacbon sơ cấp.

Các kết quả này hoàn toàn tương thích với các công trình nghiên cứu thực nghiệm của Weaver (năm 1994) và Han (năm 1996) về sự đánh đổi giữa NOx và Particulate Matter. Các dữ liệu phân bố nhiệt độ và trường nồng độ hạt trong luận văn có thể được biểu diễn một cách trực quan thông qua các đồ thị phân bố nhiệt độ theo góc quay trục khuỷu và bảng dữ liệu phát thải dạng ppm. Điểm tối ưu kỹ thuật được xác lập tại thời điểm phun 12 đến 14 degBTDC, nơi công suất động cơ đạt trên 95% giá trị cực đại trong khi tổng mức phát thải ô nhiễm được kiểm soát ở mức hài hòa nhất.

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

Nhằm ứng dụng kết quả nghiên cứu vào thực tế sản xuất và tối ưu hóa động cơ Diesel, các đề xuất cụ thể được xây dựng với lộ trình và chủ thể thực hiện rõ ràng:

  1. Hiệu chuẩn bản đồ điều khiển góc phun sớm trên ECU: Các kỹ sư lập trình động cơ tại các nhà máy sản xuất ô tô cần thiết lập dải góc phun sớm SOI tối ưu trong khoảng từ 12 đến 14 degBTDC đối với các chế độ tải trung bình và cao. Giải pháp này giúp cắt giảm 20% đến 25% nồng độ phát thải NOx mà vẫn duy trì công suất động cơ ở mức trên 95%, mục tiêu hoàn thành trong thời gian 6 tháng.
  2. Nâng cấp hệ thống phun nhiên liệu tích áp Common Rail (CRDI): Đội ngũ nghiên cứu và phát triển (R&D) của các đơn vị chế tạo động cơ cần áp dụng công nghệ Common Rail thế hệ mới với áp suất phun đạt trên 1600 bar, tích hợp thuật toán phun tách nhiều giai đoạn (phun mồi và phun chính). Biện pháp này hướng tới mục tiêu giảm 30% lượng hạt bồ hóng và rút ngắn thời gian phun DOI dưới 10 CAD, triển khai trong vòng 12 tháng.
  3. Tích hợp hệ thống tuần hoàn khí xả làm mát (Cooled EGR): Các nhà sản xuất thiết bị động lực cần phối hợp lắp đặt van EGR điều khiển điện tử với tỷ lệ hồi lưu từ 10% đến 15% để hạ nhiệt độ cháy cực đại xuống dưới 1371 độ C, giảm thêm 35% lượng NOx sinh ra trong buồng đốt, thực hiện trong thời hạn 9 tháng.
  4. Mở rộng ứng dụng phần mềm mô phỏng KIVA-3V cho nhiên liệu sinh học: Các viện nghiên cứu cơ khí giao thông và các trường đại học chuyên ngành cần mở rộng mô hình hóa buồng đốt cho các loại nhiên liệu Biodiesel B5 và B20, nhằm mục tiêu giảm thêm 10% muội than rắn và tối ưu hóa buồng cháy động cơ trong lộ trình 18 tháng.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho nhiều nhóm đối tượng:

  • Kỹ sư thiết kế và hiệu chuẩn động cơ đốt trong: Áp dụng các bảng thông số về mối quan hệ giữa SOI, DOI và tốc độ tỏa nhiệt để tối ưu hóa thuật toán điều khiển kim phun điện tử, xây dựng Fuel Map chính xác cho các dòng xe thương mại.
  • Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Cơ khí Động lực: Sử dụng luận văn như một tài liệu quy chuẩn về phương pháp thiết lập lưới 3D góc 72 độ và ứng dụng các mô hình phân rã tia phun KH-RT, dòng rối RNG k-e trong phần mềm KIVA-3V.
  • Các nhà quản lý kỹ thuật tại các nhà máy lắp ráp ô tô và máy kéo: Tham khảo giải pháp cân bằng giữa tiêu chuẩn khí thải Euro và hiệu suất tiêu hao nhiên liệu, hỗ trợ việc lựa chọn cấu hình kim phun phù hợp cho dây chuyền sản xuất.
  • Chuyên gia hoạch định chính sách và kiểm định môi trường giao thông: Khai thác các dữ liệu định lượng về phát thải NOx và muội than bồ hóng để xây dựng quy chuẩn kỹ thuật và chu kỳ kiểm định khí thải phương tiện cơ giới đường bộ.

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

Tại sao việc tăng góc phun sớm lại khiến nồng độ phát thải NOx trong khí xả tăng cao? Khi tăng góc phun sớm từ 12 lên 22 degBTDC, nhiên liệu được đưa vào buồng cháy sớm hơn, kéo dài thời gian hòa trộn trước khi tự bốc cháy. Khi bắt đầu cháy, ngọn lửa bùng phát mạnh đẩy nhiệt độ đỉnh trong xy lanh vượt ngưỡng 1371 độ C (2500 độ F). Ở mức nhiệt này, phản ứng nhiệt Zeldovich diễn ra mãnh liệt giữa nitơ và oxy tự do, làm tăng nồng độ NOx trong khí xả thêm hơn 45%.

Thời gian phun (DOI) kéo dài ảnh hưởng như thế nào đến quá trình hình thành muội than bồ hóng? Kéo dài thời gian phun từ 4 CAD lên 12 CAD làm giảm áp suất phun thực tế, khiến đường kính hạt nhiên liệu Sauter tăng trên 25 micromet, gây thiếu hụt oxy cục bộ và tăng phát thải bồ hóng. Tuy nhiên, khi tiếp tục kéo dài thời gian phun lên 14 đến 18 CAD, quá trình oxy hóa muộn trong kỳ giãn nở giúp lượng bồ hóng giảm nhẹ trở lại.

Mô hình phân rã tia phun lai KH-RT trong KIVA-3V có ưu điểm gì nổi bật so với các mô hình truyền thống? Mô hình KH-RT kết hợp đồng thời sự mất ổn định bề mặt Kelvin-Helmholtz do ma sát khí động học và sự phân tách Rayleigh-Taylor do gia tốc hạt chất lỏng. Nhờ đó, mô hình mô phỏng chính xác quá trình vỡ của tia nhiên liệu vận tốc trên 100 m/s thành các hạt siêu nhỏ dưới 10 micromet, phản ánh chân thực động học tia phun trên động cơ AVL.

Vì sao nghiên cứu lại chọn góc lưới mô phỏng 72 độ thay vì toàn bộ buồng cháy 360 độ? Đầu kim phun của động cơ nghiên cứu AVL 5402 có cấu tạo 5 lỗ phun bố trí đối xứng hoàn toàn trong không gian hình tròn. Do đó, việc cắt một phân đoạn cung tròn 72 độ (lấy 360 độ chia cho 5) để tạo lưới tính toán vừa đảm bảo đầy đủ tính chất vật lý của một tia phun, vừa giúp giảm khoảng 80% khối lượng tính toán của máy tính.

Làm thế nào để giải quyết mâu thuẫn đánh đổi giữa phát thải NOx và hạt bồ hóng? Giải pháp hiệu quả nhất là phối hợp đa công nghệ: duy trì góc phun sớm ở mức hợp lý khoảng 12 đến 14 degBTDC, kết hợp hệ thống phun nhiên liệu Common Rail áp suất cao trên 1600 bar có chế độ phun tách, và bổ sung hệ thống tuần hoàn khí xả EGR ở mức 10% đến 15%. Sự kết hợp này cho phép giảm đồng thời 20% đến 30% cả NOx lẫn bồ hóng.

Kết luận

  • Luận văn đã xây dựng thành công quy trình mô phỏng 3 chiều CFD trên phần mềm KIVA-3V ERC ứng dụng cho động cơ Diesel phun trực tiếp 1 xy lanh AVL với sai số kiểm chứng dưới 5% so với thực nghiệm.
  • Chứng minh định lượng quy luật tăng hơn 45% lượng phát thải NOx và giảm từ 35% đến 40% lượng bồ hóng khi tăng góc phun sớm từ 10 lên 22 degBTDC.
  • Phát hiện và giải thích bản chất hiện tượng biến thiên hai giai đoạn của phát thải hạt bồ hóng theo thời gian phun DOI trong dải từ 4 CAD đến 18 CAD.
  • Xác lập vùng thông số vận hành tối ưu tại góc phun sớm 12 đến 14 degBTDC, giúp cân bằng hoàn hảo giữa chỉ tiêu công suất và kiểm soát khí xả độc hại.
  • Đóng góp phương pháp luận tính toán số giúp tiết kiệm trên 60% chi phí thử nghiệm thực tế cho ngành kỹ thuật ô tô trong nước.

Trong giai đoạn tiếp theo của đề tài, hướng nghiên cứu cần tập trung mở rộng mô phỏng trên các loại nhiên liệu thay thế như Biodiesel và nghiên cứu tích hợp hệ thống tăng áp khí nạp Turbocharger trong vòng 12 tháng tới. Các đơn vị nghiên cứu và doanh nghiệp quan tâm có thể kết nối với Bộ môn Kỹ thuật Ô tô - Máy kéo để nhận chuyển giao toàn bộ quy trình mô phỏng và bộ thông số tối ưu này.