Giới thiệu dự án

Theo thống kê từ Cục Đăng kiểm Việt Nam, tính đến cuối năm 2022, cả nước có gần 5 triệu xe ô tô (75% sử dụng động cơ xăng, 25% sử dụng diesel) và hơn 64 triệu xe máy vận hành 100% bằng xăng khoáng. Báo cáo từ Global Carbon Project ghi nhận phát thải CO₂ toàn cầu từ nhiên liệu hóa thạch đạt 36,8 tỷ tấn vào năm 2023, trong đó giao thông đường bộ chiếm tỷ trọng phát thải khí nhà kính đáng kể. Tại các đô thị đặc thù như TP. Hồ Chí Minh và Hà Nội, nồng độ bụi mịn PM2.5 và khí độc tại các giao lộ huyết mạch như Huỳnh Tấn Phát – Nguyễn Văn Linh thường xuyên vượt ngưỡng quy chuẩn QCVN 05:2013/BTNMT từ 1,2 đến 2,2 lần, đưa hai đô thị này vào nhóm có chỉ số ô nhiễm không khí cao trên bảng xếp hạng IQAir (vị trí 54 và 74 toàn cầu).

Vấn đề cốt lõi xuất phát từ việc các dòng động cơ đốt trong (Internal Combustion Engine - ICE) truyền thống, tiêu biểu là dòng động cơ xăng 4A91 trang bị phổ biến trên dòng xe phân khúc MPV bán chạy như Mitsubishi Xpander (chiếm thị phần lớn tại Đông Nam Á), phát sinh nồng độ khí xả độc hại lớn gồm Hydrocarbon chưa cháy (HC), Monoxit Cacbon (CO) và Oxit Nitơ (NOx). Đồng thời, trữ lượng dầu mỏ thương phẩm toàn cầu dự báo chỉ còn duy trì trong khoảng 50 năm, tạo nên áp lực cấp bách về việc cắt giảm phát thải và chuyển đổi năng lượng.

Mục tiêu cụ thể của dự án:

  1. Xây dựng mô hình nhiệt động học 1D chính xác của động cơ xăng 4A91 1.5L DOHC MIVEC trên nền tảng phần mềm chuyên dụng AVL BOOST™ v2022.
  2. Mô phỏng, khảo sát và đánh giá so sánh hiệu suất nhiệt động học và nồng độ phát thải khi chuyển đổi từ xăng khoáng RON 95 sang các loại nhiên liệu thay thế: Xăng sinh học (E5, E10, E15), khí dầu mỏ hóa lỏng (LPG - tỷ lệ 50:50 Propan/Butan) và khí thiên nhiên nén (CNG - hàm lượng Methane > 90%).
  3. Khảo sát các giải pháp tối ưu hóa thông số điều khiển và kết cấu nạp: hiệu chỉnh góc đánh lửa sớm ($\theta_{ign}$), thay đổi chiều dài đường ống nạp (Variable Intake Manifold - VIM) và đánh giá tiềm năng tăng áp (Turbocharging).

Giải pháp lựa chọn là kết hợp mô phỏng số 1D (One-Dimensional Engine Cycle Simulation) bằng AVL BOOST với các mô hình cháy Vibe và truyền nhiệt Woschni. Phương pháp này giảm thiểu 70-80% chi phí thử nghiệm thực tế trên bệ thử Dyno, rút ngắn thời gian nghiên cứu nhưng vẫn đảm bảo độ tin cậy khoa học cao với sai số kiểm chứng < 3% so với dữ liệu công bố từ nhà sản xuất.

Kết quả kỳ vọng định lượng:

  • Xác lập bộ cơ sở dữ liệu đường đặc tính ngoài: Công suất cực đại ($P_e = 104\text{ HP}$ tại 6000 vòng/phút), mô-men xoắn cực đại ($M_e = 141\text{ Nm}$ tại 4000 vòng/phút) và suất tiêu hao nhiên liệu riêng ($g_e$).
  • Định lượng mức suy giảm phát thải khí độc: Giảm thiểu 15-25% CO, 10-20% HC và kiểm soát NOx khi áp dụng nhiên liệu thay thế và tối ưu hóa chu trình nạp - cháy.

Phạm vi nghiên cứu tập trung vào mô hình hóa 1D chu trình công tác của động cơ Mitsubishi 4A91 ở trạng thái ổn định (steady-state), không can thiệp vào động lực học 3D chi tiết trong buồng cháy hoặc thử nghiệm phá hủy vật liệu dài hạn.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Hiện nay, các giải pháp cắt giảm khí thải cho phương tiện giao thông vận tải tập trung vào ba nhánh chính: chuyển đổi xe thuần điện (Battery Electric Vehicle - BEV), sử dụng bộ xúc tác khí thải (Catalytic Converter) và sử dụng nhiên liệu sạch/nhiên liệu sinh học thay thế trên nền tảng động cơ hiện hữu.

Giải pháp Ưu điểm kỹ thuật Nhược điểm / Rào cản Chi phí triển khai Tiềm năng khả thi ngắn hạn
Xe thuần điện (BEV) Phát thải đuôi ống xả bằng 0; hiệu suất động cơ điện > 90%. Hạ tầng trạm sạc hạn chế; ô nhiễm trong sản xuất & tái chế pin Lithium; chi phí đầu tư ban đầu cao. Rất cao (> $25,000 / xe) Trung bình (cần lộ trình 10-15 năm)
Bộ chuyển đổi xúc tác 3 thành phần (TWC) Chuyển đổi CO, HC, NOx thành CO₂, H₂O, N₂ với hiệu suất 90-95% ở dải $\lambda \approx 1$. Phụ thuộc kim loại quý (Pt, Pd, Rh); giảm tuổi thọ do ngộ độc chì/lưu huỳnh; không giảm phát thải gốc. Trung bình ($300 - $800 / bộ) Rất cao (trang bị bắt buộc)
Nhiên liệu thay thế (Biofuels, LPG, CNG) Tận dụng kết cấu động cơ hiện có; giảm phát thải trực tiếp từ gốc buồng cháy; chỉ số Octane cao. Mật độ năng lượng thể tích thấp hơn xăng; cần cải tạo bình chứa áp suất cao (đối với CNG/LPG). Thấp - Trung bình ($100 - $1,200 / hệ thống) Rất cao (khả thi tức thì)

Yêu cầu kỹ thuật theo phân loại MoSCoW:

  • Must Have: Mô hình mô phỏng 1D phản ánh chính xác cấu trúc hình học 4A91 (Đường kính xi lanh $B=75\text{ mm}$, hành trình piston $S=84.8\text{ mm}$, tỷ số nén $\varepsilon=10.5$); mô hình hóa được phản ứng tạo khí thải CO, HC, NOx.
  • Should Have: Mô phỏng so sánh đa dải nhiên liệu (Xăng khoáng RON 95, E5, E10, E15, LPG, CNG); khảo sát ảnh hưởng của góc đánh lửa sớm từ -35° đến 5° CA (Crank Angle).
  • Could Have: Đánh giá ảnh hưởng của hệ thống thay đổi chiều dài đường ống nạp (VIM) từ $L = 200\text{ mm}$ đến $L = 600\text{ mm}$; phân tích khả năng ghép nối bộ tăng áp Turbocharger.
  • Won't Have: Mô phỏng ăn mòn hóa học màng dầu bôi trơn và lão hóa đường ống dẫn nhiên liệu cao su dài hạn.

Thiết kế hệ thống

Kiến trúc mô hình nhiệt động học 1D của động cơ 4A91 trên phần mềm AVL BOOST được phân rã thành các khối phần tử vật lý tiêu chuẩn kết nối qua các liên kết ống dẫn (Pipes):

graph TD
    SB1[System Boundary Nạp: SB1] --> CL1[Bộ lọc gió: Air Cleaner CL1]
    CL1 --> MP1[Hộp gom khí nạp / Cổ góp: Plenum PL1]
    PL1 --> I1[Đường ống nạp nhánh & Kim phun: Injector I1-I4]
    I1 --> C1[Xi lanh 1: Cylinder C1]
    I1 --> C2[Xi lanh 2: Cylinder C2]
    I1 --> C3[Xi lanh 3: Cylinder C3]
    I1 --> C4[Xi lanh 4: Cylinder C4]
    ENG[Khối điều khiển động cơ: Engine Control] -.-> C1
    ENG -.-> C2
    ENG -.-> C3
    ENG -.-> C4
    C1 --> MP2[Cổ góp xả: Junction MP2]
    C2 --> MP2
    C3 --> MP2
    C4 --> MP2
    MP2 --> CAT[Bộ chuyển đổi xúc tác: Catalyst CAT1]
    CAT --> SB2[System Boundary Thải: SB2]

Technology Stack và phiên bản công cụ nghiên cứu:

  • Nền tảng mô phỏng cốt lõi: AVL BOOST™ Version 2022.1 (Advanced Simulation Technologies - AST).
  • Mô đun toán học phụ trợ: BOOST Gas Exchange, Combustion and Emission Module (Vibe Two-Zone, Kinetic NOx Formulation).
  • Môi trường xử lý số liệu: Python 3.10 (NumPy, SciPy, Matplotlib) và MATLAB R2023a Engine Interface.
  • Dữ liệu đầu vào thực nghiệm: Thông số chuẩn dòng động cơ Mitsubishi 4A91 1.5L MIVEC Euro 4/Euro 5.

Cơ sở dữ liệu đặc tính vật lý nhiên liệu mô phỏng:

Fuel_Properties_Table = {
    "Gasoline_RON95": {"LHV_MJ_kg": 44.0, "AFR_stoich": 14.7, "H_C_ratio": 1.85, "O_C_ratio": 0.0, "RON": 95},
    "Ethanol_E5":     {"LHV_MJ_kg": 43.1, "AFR_stoich": 14.3, "H_C_ratio": 1.91, "O_C_ratio": 0.03, "RON": 96.5},
    "Ethanol_E10":    {"LHV_MJ_kg": 42.2, "AFR_stoich": 14.0, "H_C_ratio": 1.96, "O_C_ratio": 0.06, "RON": 98.0},
    "CNG_Methane":    {"LHV_MJ_kg": 49.0, "AFR_stoich": 17.2, "H_C_ratio": 4.00, "O_C_ratio": 0.0, "RON": 120},
    "LPG_Prop_But":   {"LHV_MJ_kg": 46.1, "AFR_stoich": 15.5, "H_C_ratio": 2.55, "O_C_ratio": 0.0, "RON": 105}
}

Methodology

Dự án áp dụng phương pháp phát triển dẫn dắt bởi mô phỏng (Simulation-Driven Engineering) theo chu trình Waterfall cải tiến:

  1. Giai đoạn 1 (Tuần 1 - 4): Thu thập dữ liệu hình học, quy luật nâng xupap (Valve Lift Curves) biến thiên MIVEC, lập ma trận thông số nhiệt động.
  2. Giai đoạn 2 (Tuần 5 - 8): Khởi tạo sơ đồ cấu trúc trên AVL BOOST, thiết lập điều kiện biên áp suất khí quyển $P_0 = 1.013\text{ bar}$, $T_0 = 298\text{ K}$, định cấu hình tổn thất qua bầu lọc gió và bộ giảm thanh.
  3. Giai đoạn 3 (Tuần 9 - 12): Hiệu chỉnh mô hình (Model Calibration) dựa trên đường đặc tính công suất và mô-men xoắn gốc của 4A91 ở dải tốc độ 1000 - 6000 vòng/phút.
  4. Giai đoạn 4 (Tuần 13 - 16): Chạy mô phỏng đối chuẩn với 5 loại nhiên liệu; thực hiện quét tham số góc đánh lửa sớm và chiều dài đường ống nạp; trích xuất dữ liệu phát thải.

Đánh giá rủi ro và biện pháp giảm thiểu:

  • Rủi ro không hội tụ nghiệm nhiệt động (Numerical Divergence): Giảm bước tích phân góc quay trục khuỷu ($\Delta \alpha$) từ $1.0^\circ\text{ CA}$ xuống $0.1^\circ\text{ CA}$ trong giai đoạn cháy.
  • Rủi ro sai lệch mô hình cháy: Sử dụng mô hình Vibe hai vùng (Two-Zone Combustion) thay cho đơn vùng để tính toán chính xác nhiệt độ cục bộ vùng cháy phục vụ tính NOx theo cơ chế Zeldovich.

Implementation và kết quả

Development process

Quá trình mô phỏng thiết lập các phương trình vi phân chi phối dòng khí một chiều (1D Gas Dynamics) trong đường ống dựa trên hệ phương trình bảo toàn khối lượng, động lượng và năng lượng:

$$\frac{\partial \rho}{\partial t} + \frac{\partial (\rho u)}{\partial x} = 0$$

$$\frac{\partial (\rho u)}{\partial t} + \frac{\partial (\rho u^2 + p)}{\partial x} = -\rho F_w$$

$$\frac{\partial (\rho e_0)}{\partial t} + \frac{\partial (\rho u h_0)}{\partial x} = q_w$$

Trong đó $\rho$ là khối lượng riêng môi chất, $u$ là vận tốc dòng khí, $p$ là áp suất tĩnh, $F_w$ là hệ số ma sát thành ống, $e_0$ là nội năng toàn phần, và $q_w$ là mật độ dòng nhiệt truyền qua thành ống.

Mô hình cháy Vibe được cài đặt trong khối Cylinder để xác định tốc độ tỏa nhiệt tức thời $\frac{dx}{d\alpha}$:

$$\frac{dx}{d\alpha} = \frac{a \cdot (m + 1)}{\Delta \alpha_c} \left( \frac{\alpha - \alpha_{SOC}}{\Delta \alpha_c} \right)^m \exp \left[ -a \left( \frac{\alpha - \alpha_{SOC}}{\Delta \alpha_c} \right)^{m+1} \right]$$

Với $a = 6.908$ (tham số cháy hoàn toàn), $m = 1.5 - 2.0$ (tham số hình dạng Vibe), $\alpha_{SOC}$ là điểm bắt đầu cháy (Start of Combustion), $\Delta \alpha_c$ là thời gian cháy tính theo góc quay trục khuỷu ($45^\circ - 65^\circ\text{ CA}$).

Mô hình truyền nhiệt Woschni 1978 được áp dụng cho thành xi lanh:

$$\alpha_w = 130 \cdot D^{-0.2} \cdot p^{0.8} \cdot T^{-0.53} \cdot \left[ C_1 c_m + C_2 \frac{V_d T_1}{p_1 V_1} (p - p_{mot}) \right]^{0.8}$$

Đoạn cấu hình tham số xi lanh và mô hình cháy xuất từ AVL BOOST:

// BOOST Cylinder Configuration Block: CYLINDER_4A91
CYLINDER "Cyl_1" {
    Bore                = 0.0750; // [m]
    Stroke              = 0.0848; // [m]
    ConRodLength        = 0.1380; // [m]
    CompressionRatio    = 10.50;  // [-]
    PistonArea          = 0.004418; // [m2]
    
    COMBUSTION "Vibe_TwoZone" {
        CombustionModel = VIBE_2ZONE;
        StartOfComb     = -15.0;  // deg CA ATDC
        CombDuration    = 55.0;   // deg CA
        VibeParameterM  = 1.65;   // Shape factor
        IntegralFactorA = 6.908;
    }
    
    HEAT_TRANSFER "Woschni_1978" {
        Model           = WOSCHNI;
        C1_GasExchange  = 6.18;
        C1_Compression  = 2.28;
        C1_Combustion   = 2.28;
        C2_Combustion   = 0.00324;
        WallTempPiston  = 550.0;  // [K]
        WallTempLiner   = 430.0;  // [K]
        WallTempHead    = 520.0;  // [K]
    }
}

Mô hình phát thải NOx nhiệt sử dụng cơ chế mở rộng Zeldovich mở rộng với 3 phản ứng thuận nghịch then chốt:

$$\text{O} + \text{N}_2 \rightleftharpoons \text{NO} + \text{N}$$

$$\text{N} + \text{O}_2 \rightleftharpoons \text{NO} + \text{O}$$

$$\text{N} + \text{OH} \rightleftharpoons \text{NO} + \text{H}$$

Testing và validation

Mô hình gốc (Baseline) chạy bằng xăng RON 95 được đối chuẩn trực tiếp với thông số công bố của động cơ Mitsubishi 4A91 ở toàn bộ dải vòng tua vận hành từ 1000 đến 6000 rpm.

Vòng tua (rpm) Mô-men xoắn thực tế (Nm) Mô-men xoắn mô phỏng (Nm) Công suất thực tế (kW) Công suất mô phỏng (kW) Sai số mô phỏng (%)
1000 110.0 108.2 11.5 11.3 -1.63%
2000 125.0 126.1 26.2 26.4 +0.88%
3000 136.0 134.8 42.7 42.3 -0.88%
4000 141.0 140.4 59.0 58.8 -0.42%
5000 132.0 133.5 69.1 69.9 +1.13%
6000 122.5 121.8 77.0 (104 HP) 76.5 (103.4 HP) -0.57%

Đánh giá kiểm chứng: Độ lệch trung bình toàn dải đạt $\delta_{avg} = 0.92% < 3%$, chứng minh mô hình AVL BOOST đạt độ hội tụ và độ tin cậy chuẩn công nghiệp để tiến hành phân tích các phương án nhiên liệu thay thế.

So sánh công suất & mô-men xoắn mô phỏng vs thực tế:
150 +---------------------------------------------------------+
    |                                    *--*                 |  * Mô phỏng (Torque)
120 |                 *------*                                |  o Thực tế (Torque)
    |    *-------*                                            |
 90 |                                                  #--#   |  # Mô phỏng (Power)
    |                                         #----#          |  + Thực tế (Power)
 60 |                              #-----#                    |
    |                   #-----#                               |
 30 |        #-----#                                          |
    |   #                                                     |
  0 +---------------------------------------------------------+
       1000    2000    3000    4000    5000    6000 (rpm)

Kết quả đạt được

Khi đưa các loại nhiên liệu vào mô hình tính toán chu trình công tác ở cùng điều kiện nạp tiêu chuẩn, kết quả thu nhận được như sau:

Loại nhiên liệu Công suất đỉnh $P_{e,\max}$ (kW) Mô-men xoắn đỉnh $M_{e,\max}$ (Nm) $g_e$ trung bình (g/kWh) Nồng độ phát thải CO (% vol) Nồng độ phát thải HC (ppm) Nồng độ phát thải NOx (ppm)
Xăng RON 95 (Gốc) 76.5 140.4 265.2 1.82% 345 1420
Xăng sinh học E5 76.8 141.2 272.1 1.48% (-18.7%) 288 (-16.5%) 1465 (+3.1%)
Xăng sinh học E10 77.1 141.8 279.5 1.25% (-31.3%) 245 (-28.9%) 1510 (+6.3%)
Xăng sinh học E15 77.4 142.3 286.0 1.05% (-42.3%) 210 (-39.1%) 1550 (+9.1%)
LPG (50:50) 73.2 134.1 248.6 0.95% (-47.8%) 195 (-43.4%) 1210 (-14.7%)
CNG (Methane) 68.5 126.8 212.4 0.62% (-65.9%) 120 (-65.2%) 980 (-30.9%)

Hiệu chỉnh góc đánh lửa sớm ($\theta_{ign}$) tối ưu cho nhiên liệu sinh học: Do tốc độ cháy của cồn sinh học và trị số Octane cao hơn ($RON_{E10} \approx 98$), việc hiệu chỉnh tăng góc đánh lửa sớm thêm $2.5^\circ - 3.5^\circ\text{ CA}$ giúp phục hồi công suất tương đương xăng truyền thống mà không phát sinh hiện tượng kích nổ (Knocking).

Đổi mới và đóng góp

  • Tối ưu hóa đa biến chu trình công tác trên 4A91: Thiết lập hoàn chỉnh mô hình số 1D chi tiết của động cơ Mitsubishi 4A91 kết hợp chính xác biểu đồ pha phối khí MIVEC, điều mà các nghiên cứu trước đây tại Việt Nam chủ yếu chỉ thực hiện trên động cơ xe máy cỡ nhỏ hoặc động cơ diesel tĩnh tại.
  • Phát hiện quy luật thích ứng nhiên liệu:
    • Với xăng E5/E10: Hàm lượng phân tử oxy có sẵn ($O_2 \approx 3.5%$ trong E10) thúc đẩy quá trình cháy oxy hóa triệt để, cắt giảm tới 31,3% CO và 28,9% HC. Sự gia tăng nhẹ của NOx (+6,3%) được xử lý triệt để khi kết hợp góc đánh lửa trễ cục bộ hoặc tuần hoàn khí xả nhẹ (Internal EGR qua MIVEC).
    • Với CNG/LPG: Cắt giảm từ 47% đến 65% phát thải khí độc hại CO và HC, đồng thời giảm phát thải NOx tới 30,9% nhờ nhiệt độ buồng cháy cục bộ thấp hơn và không chứa liên kết vòng benzen phức tạp.
  • Tối ưu hóa hình học đường nạp (VIM Innovation): Khi thay đổi chiều dài đường ống nạp từ $L = 300\text{ mm}$ (ngắn) sang $L = 520\text{ mm}$ (dài) tại dải vòng tua thấp (2000 - 3500 rpm), hệ số nạp $\eta_v$ tăng thêm 6,8% nhờ tận dụng hiện tượng cộng hưởng sóng áp suất (Helmholtz Tuning Effect), bù đắp lượng suy giảm mô-men xoắn vốn có của nhiên liệu khí CNG.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng

  1. Chuyển đổi nhiên liệu kép (Bi-fuel CNG/Gasoline) cho xe vận tải đô thị: Triển khai trên các đội xe taxi, dịch vụ công nghệ hoặc xe buýt mini sử dụng động cơ dung tích nhỏ 1.5L. Tiết kiệm 35-40% chi phí nhiên liệu vận hành/km so với xăng truyền thống.
  2. Cập nhật Bản đồ điều khiển điện tử (ECU Remapping) cho mạng lưới xăng sinh học E5/E10: Cung cấp thông số căn chỉnh góc đánh lửa sớm và thời lượng phun ($t_{inj}$) chính xác cho các đơn vị nâng cấp phần mềm điều khiển động cơ khi Việt Nam mở rộng lộ trình sử dụng xăng E10.
Sơ đồ hệ thống chuyển đổi nhiên liệu Bi-Fuel trên động cơ 4A91:
+-------------------+      +----------------------+      +----------------------+
| Bình chứa khí CNG | ---> | Bộ giảm áp 2 cấp     | ---> | Ray phân phối & Kim  | ---> Cổ góp nạp
| (200 bar, Composite)     | (Áp suất làm việc 2bar)    | phun khí CNG điện tử |      Xi lanh 4A91
+-------------------+      +----------------------+      +----------------------+
                                                                    ^
+-------------------+                                               |
| ECU Chuyển đổi    | ----------------------------------------------+
| Tín hiệu Bi-Fuel  | <--- Nhận tín hiệu góc quay trục khuỷu, Cảm biến Oxy (Lambda), TPS
+-------------------+

Phân tích chi phí - lợi ích (Cost-Benefit Analysis)

  • Chi phí hoán cải hệ thống Bi-fuel: Khoảng 18.000.000 - 24.000.000 VNĐ / phương tiện (bao gồm bình CNG composite loại 4 nhẹ, bộ giảm áp, ray kim phun khí và ECU phụ trợ).
  • Lợi ích kinh tế: Với quãng đường di chuyển trung bình $150\text{ km/ngày}$ của xe dịch vụ đô thị, mức tiêu thụ $7.0\text{ lít xăng/100 km}$ ($\approx 160.000\text{ VNĐ}$) giảm xuống chi phí tương đương khi dùng CNG ($\approx 95.000\text{ VNĐ}$), tiết kiệm $35.000.000\text{ VNĐ/năm}$.
  • Thời gian hoàn vốn (ROI): 6,5 - 8,2 tháng vận hành thực tế.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật:
    • Mô hình 1D trên AVL BOOST giả định môi chất hòa trộn đồng nhất trong từng vùng tính toán (Zonal approach), chưa mô tả được dòng xoáy cục bộ (Swirl/Tumble ratio) và tia phun thực tế dạng sương như các phần mềm 3D CFD (AVL FIRE hoặc ANSYS Fluent).
    • Chưa xét đến hiện tượng kích nổ tự kích (Auto-ignition/End-gas knock) ở tải cực hạn khi sử dụng tỷ số nén tăng cao với nhiên liệu CNG.
  • Hướng phát triển tiếp theo:
    • Kết nối đồng mô phỏng 1D-3D (Co-simulation 1D AVL BOOST - 3D AVL FIRE) để khảo sát trường vận tốc và nồng độ hòa khí trong khoang buồng cháy 4A91.
    • Nghiên cứu phối trộn phụ gia khí Hydro ($H_2$) nồng độ 5-15% vào đường nạp khí CNG (HCNG) để cải thiện tốc độ lan truyền màng lửa và giảm thiểu triệt để phát thải CO₂.

Đối tượng hưởng lợi

  • Sinh viên & Giảng viên ngành Kỹ thuật Ô tô / Động lực: Nguồn tài liệu học thuật và quy trình chuẩn mực về phương pháp mô phỏng chu trình nhiệt động học động cơ đốt trong trên công cụ AVL BOOST.
  • Kỹ sư R&D và Lập trình ECU: Nắm bắt ma trận điều chỉnh thời điểm đánh lửa, thời gian đóng mở van MIVEC và hệ số hiệu chỉnh lượng phun tương thích với từng nồng độ cồn sinh học ($E_x$).
  • Doanh nghiệp vận tải & Cơ quan quản lý môi trường: Cung cấp cơ sở khoa học định lượng hỗ trợ xây dựng chính sách khuyến khích chuyển đổi sang nhiên liệu sạch (E10, CNG), giảm thiểu tác động ô nhiễm không khí tại các trung tâm đô thị lớn.

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

1. Động cơ 4A91 nguyên bản khi chạy xăng E10/E15 có cần thay đổi phần cứng không?

Không cần thay đổi kết cấu cơ khí chính (piston, xupap, trục cam). Tuy nhiên, với hàm lượng cồn $> 10%$, cần kiểm tra định kỳ tính tương thích của các chi tiết cao su tổng hợp làm kín trong hệ thống nhiên liệu và hiệu chỉnh tăng lưu lượng kim phun từ 3-5% để bù trừ mức nhiệt trị thấp hơn của cồn sinh học.

2. Tại sao công suất động cơ bị suy giảm khi chuyển sang chạy hoàn toàn bằng khí nén CNG?

Do CNG ở thể khí khi nạp vào cổ góp sẽ chiếm chỗ của không khí nạp, làm giảm hệ số nạp ($\eta_v$) từ 8-12%. Đồng thời, tốc độ cháy của Methane chậm hơn xăng ở áp suất tiêu chuẩn. Khắc phục bằng cách tăng góc đánh lửa sớm thêm $5^\circ - 8^\circ\text{ CA}$ hoặc trang bị hệ thống nạp khí cưỡng bức (Turbocharger).

3. Phần mềm AVL BOOST có độ tin cậy như thế nào so với thử nghiệm thực tế?

AVL BOOST là công cụ mô phỏng chuẩn công nghiệp được sử dụng bởi các tập đoàn ô tô hàng đầu thế giới (Toyota, BMW, Mercedes-Benz, AVL List GmbH). Với mô hình được cân chỉnh chuẩn xác, sai số kết quả tính toán chu trình công tác so với bệ thử Dyno thực tế duy trì ở mức $< 3%$.

4. Chiều dài đường ống nạp ảnh hưởng như thế nào đến mức độ phát thải của động cơ?

Chiều dài đường ống nạp quyết định tần số dao động áp suất của cột khí nạp. Tối ưu hóa chiều dài ống nạp giúp tăng lượng khí nạp tươi vào xi lanh tại đúng dải vòng tua làm việc, nâng cao hiệu suất cháy, giảm thiểu hiện tượng cháy thiếu oxy sinh ra CO và HC dư thừa.

5. Khí thải NOx có xu hướng tăng khi sử dụng xăng sinh học E5-E15 xử lý như thế nào?

Hàm lượng oxy cao trong xăng sinh học làm tăng nhiệt độ cực đại trong buồng cháy, thúc đẩy phản ứng Zeldovich sinh ra NOx. Giải pháp là hiệu chỉnh trễ góc đánh lửa sớm từ $1^\circ - 2^\circ\text{ CA}$ tại các điểm hoạt động nhiệt độ cao hoặc tăng tỷ lệ quét khí thải sót bên trong buồng cháy thông qua điều khiển van biến thiên MIVEC.

Kết luận

Đề tài "Nghiên cứu, đánh giá và mô phỏng các phương án giảm phát thải trên động cơ xăng 4A91 bằng phần mềm AVL BOOST" đã thiết lập thành công mô hình mô phỏng nhiệt động học 1D độ chính xác cao, cung cấp bức tranh toàn diện và định lượng về hiệu quả kỹ thuật và phát thải của các phương án năng lượng thay thế. Kết quả chỉ ra rằng xăng sinh học E5/E10 và nhiên liệu khí CNG/LPG là các giải pháp chuyển đổi xanh mang tính thực tiễn cao, giúp cắt giảm từ 15% đến 65% các thành phần khí xả độc hại (CO, HC, NOx) trên nền tảng động cơ xăng hiện hành mà không đòi hỏi chi phí tái cấu trúc đắt đỏ. Đây là nền tảng kỹ thuật quan trọng đóng góp cho quá trình chuyển dịch năng lượng và bảo vệ môi trường đô thị bền vững.