Giới thiệu dự án

Bối cảnh ngành và thực trạng ô nhiễm khí thải

Trong kỷ nguyên công nghiệp hóa toàn cầu, ngành giao thông vận tải phụ thuộc chủ yếu vào động cơ đốt trong (ĐCĐT) sử dụng nhiên liệu hóa thạch có nguồn gốc từ dầu mỏ. Theo các báo cáo năng lượng toàn cầu, lượng tiêu thụ dầu mỏ cho giao thông vận tải chiếm hơn 60% tổng sản lượng khai thác, đồng thời là nguồn phát thải chính gây ra hơn 24% lượng khí thải nhà kính ($CO_2$), bên cạnh các chất độc hại như $CO$, $HC$, $NO_x$ và bụi mịn.

Tại Việt Nam, sự gia tăng nhanh chóng của phương tiện cơ giới đường bộ (hơn 4,37 triệu xe ô tô lưu hành) đặt ra áp lực nặng nề lên chất lượng không khí tại các đô thị lớn. Nhằm giải quyết bài toán an ninh năng lượng và giảm phát thải, Thủ tướng Chính phủ đã ban hành Quyết định số 177/2007/QĐ-TTg và Quyết định số 53/2012/QĐ-TTg quy định lộ trình áp dụng tỷ lệ phối trộn nhiên liệu sinh học với nhiên liệu truyền thống (xăng sinh học E5, E10).

   ┌────────────────────────────────────────────────────────┐
   │             ÁP LỰC CHUYỂN ĐỔI NĂNG LƯỢNG               │
   ├──────────────────────────┬─────────────────────────────┤
   │ Cạn kiệt dầu mỏ hóa thạch│ Ô nhiễm khí xả (CO, HC, NOx)│
   └─────────────┬────────────┴──────────────┬──────────────┘
                 │                           │
                 ▼                           ▼
   ┌────────────────────────────────────────────────────────┐
   │   GIẢI PHÁP: NHIÊN LIỆU SINH HỌC XĂNG - ETHANOL (E0-E85)│
   ├────────────────────────────────────────────────────────┤
   │  Nghiên cứu mô phỏng nhiệt động lực học 1D (AVL BOOST)  │
   │  Đánh giá tính năng kỹ thuật & Tối ưu hóa phát thải     │
   └────────────────────────────────────────────────────────┘

Vấn đề nghiên cứu cụ thể (Problem Statement)

Ethanol sinh học ($C_2H_5OH$) có các đặc tính hóa lý khác biệt rõ rệt so với xăng khoáng truyền thống: chỉ số Octane cao hơn ($RON = 106 - 130$), hàm lượng oxy liên kết đạt xấp xỉ $34.7%$ khối lượng và nhiệt ẩn hóa hơi lớn ($840\text{ kJ/kg}$ so với $380 - 400\text{ kJ/kg}$ của xăng). Tuy nhiên, nhiệt trị thấp của ethanol chỉ đạt khoảng $26.8\text{ - }29.7\text{ MJ/kg}$ (thấp hơn khoảng $32%$ so với xăng khoáng $43\text{ - }44\text{ MJ/kg}$).

Việc sử dụng trực tiếp các hỗn hợp xăng sinh học có tỷ lệ ethanol cao (từ E15 đến E85) trên các dòng động cơ xăng thương mại chưa hiệu chỉnh sẽ dẫn đến hiện tượng:

  • Suy giảm công suất và mô-men xoắn hiệu dụng do thiếu hụt năng lượng nạp.
  • Gia tăng đột biến suất tiêu hao nhiên liệu có ích ($BSFC$).
  • Khó khởi động nguội do nhiệt độ bay hơi của ethanol cao ($78^\circ\text{C}$).
  • Nguy cơ ăn mòn vật liệu hệ thống nhiên liệu và biến động nồng độ phát thải $NO_x$.

Việc thử nghiệm thực nghiệm trên bệ thử động cơ động lực học (Dynamometer Testbench) với từng cấp độ phối trộn (E0, E5, E10, E15, E30, E50, E85) đòi hỏi chi phí đầu tư thiết bị rất cao, tốn thời gian và có nguy cơ gây hỏng hóc cơ khí. Do đó, bài toán đặt ra là cần một công cụ mô phỏng số nhiệt động lực học 1D có độ tin cậy cao để dự báo chính xác chu trình công tác, đặc tính cháy và phát thải.

Mục tiêu dự án

  1. Khảo sát cơ sở lý thuyết nhiệt động: Phân tích chi tiết cơ chế phản ứng cháy, truyền nhiệt và biến thiên tính chất hóa lý của hỗn hợp xăng - ethanol từ E0 đến E85.
  2. Xây dựng mô hình số 1D hoàn chỉnh: Thiết lập mô hình mô phỏng động cơ xăng 4 kỳ, 4 xylanh Toyota 1NZ-FE 1.5L trên phần mềm chuyên ngành AVL BOOST 2020.
  3. Phối trộn và cấu hình dữ liệu nhiên liệu: Ứng dụng công cụ BOOST Gas Properties Tool để định nghĩa chính xác các thành phần nhiên liệu sinh học E0, E5, E10, E15, E30, E50 và E85.
  4. Đánh giá định lượng tính năng và phát thải: Phân tích các chỉ số công suất ($N_e$), mô-men xoắn ($M_e$), suất tiêu hao nhiên liệu ($g_e$), áp suất buồng cháy ($p_c$) và nồng độ các chất phát thải ($CO, HC, NO_x$).
  5. Đề xuất giải pháp thích ứng kết cấu: Đưa ra các khuyến nghị hiệu chỉnh động cơ và hệ thống điều khiển điện tử (ECU) cho phương tiện sử dụng nhiên liệu linh hoạt (Flexible Fuel Vehicle - FFV).

Phương pháp tiếp cận giải pháp

Dự án áp dụng phương pháp mô phỏng số 1 chiều (1D Thermodynamic Cycle Simulation) dựa trên việc tích hợp phương trình bảo toàn khối lượng, định luật nhiệt động lực học thứ nhất, mô hình cháy Vibe 2 vùng (Vibe 2-Zone Combustion Model), mô hình truyền nhiệt Woschni 1978 và mô hình động học hóa học mở rộng Zeldovich để tính toán phát thải $NO_x$ cùng cơ chế cân bằng hóa học cho $CO/HC$.

Kết quả kỳ vọng và chỉ số đo lường

  • Mô hình mô phỏng động cơ 1NZ-FE đạt sai số tương đối so với dữ liệu thực nghiệm chuẩn dưới $3.5%$ về công suất và mô-men xoắn trên toàn dải vòng tua 1000 - 6000 RPM.
  • Xác định chính xác quy luật biến thiên của công suất ($N_e$ giảm từ $1.2%$ đến $13.5%$) và $BSFC$ (tăng từ $2.5%$ đến $49.2%$) khi tăng hàm lượng ethanol từ E5 lên E85.
  • Định lượng mức độ cải thiện môi trường: Phát thải $CO$ giảm tối đa trên $30%$, $HC$ giảm trên $25%$ ở chế độ toàn tải.

Phạm vi và giới hạn nghiên cứu

  • Đối tượng: Động cơ xăng Toyota 1NZ-FE (1.5L, 4 xylanh thẳng hàng, 16 van DOHC, phun xăng đa điểm gián tiếp qua đường nạp - PFI, tỷ số nén $\varepsilon = 10.5:1$).
  • Điều kiện mô phỏng: Chế độ toàn tải (Wide Open Throttle - WOT), dải tốc độ từ 1000 RPM đến 6000 RPM (bước nhảy 1000 RPM), điều kiện biên khí quyển tiêu chuẩn ($p_0 = 1\text{ bar}, T_0 = 298\text{ K}$).
  • Giới hạn: Không can thiệp thay đổi tỷ số nén cơ khí hoặc tối ưu hóa góc đánh lửa sớm theo từng loại nhiên liệu trong mô hình cơ sở; đánh giá ở trạng thái xác lập (steady-state).

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

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

Tiêu chí Thử nghiệm thực nghiệm (Dynamometer Testbench) Mô phỏng 3D CFD (AVL FIRE / ANSYS Forte) Mô phỏng nhiệt động 1D (AVL BOOST)
Chi phí triển khai Rất cao (> 500 triệu VNĐ bệ thử & cảm biến) Cao (Yêu cầu trạm làm việc HPC cấu hình mạnh) Tối ưu, phù hợp nghiên cứu ứng dụng R&D
Thời gian tính toán Thời gian chuẩn bị và đo kiểm kéo dài nhiều tuần Vài ngày/chu trình cho 1 điểm làm việc Vài chục giây đến vài phút/chu trình
Độ linh hoạt cấu hình Hạn chế, khó thay đổi tức thời thành phần nhiên liệu Trung bình, đòi hỏi chia lưới lại hình học Rất cao, tham số hóa nhanh chóng qua GUI
Độ chính xác chu trình Thực tế, chịu ảnh hưởng của sai số đo Rất chi tiết trường dòng chảy & nồng độ cục bộ Độ chính xác cao về đặc tính vĩ mô (P, T, Ne, Me)
Độ an toàn cơ khí Nguy cơ kích nổ, ăn mòn kim phun khi chạy E85 An toàn tuyệt đối trên mô hình số An toàn tuyệt đối, tối ưu hóa trước khi thử nghiệm

Yêu cầu kỹ thuật theo mô hình MoSCoW

  • Must have (Bắt buộc có): Mô hình hóa đầy đủ 4 xylanh, đường ống nạp/xả, bình ổn áp (Plenum), kim phun, mô hình cháy Vibe 2-Zone, dữ liệu nhiệt động các cấp độ E0, E5, E10, E15, E30, E50, E85.
  • Should have (Nên có): Tích hợp mô hình phát thải chi tiết (Extended Zeldovich cho $NO_x$, Kinetic model cho $CO$), mô hình tổn thất ma sát phụ thuộc tốc độ (Friction model).
  • Could have (Có thể có): Kết nối đồng mô phỏng (Co-simulation) qua MATLAB/Simulink DLL để điều khiển góc đánh lửa sớm theo hàm lượng cồn.
  • Won't have (Chưa thực hiện đợt này): Mô phỏng dòng chảy 3 chiều chi tiết bên trong buồng cháy và mô phỏng chu trình chuyển tiếp tức thời (Transient driving cycles).

Thiết kế hệ thống

Sơ đồ cấu trúc phần tử mô phỏng trên AVL BOOST

Hệ thống động cơ Toyota 1NZ-FE được chuyển đổi thành mạng lưới phần tử 1D liên kết:

graph LR
    SB1[System Boundary 1: Khí nạp] --> CL1[Air Cleaner: Lọc gió]
    CL1 --> PL1[Plenum: Bình ổn áp]
    PL1 --> P1[Pipe 1] & P2[Pipe 2] & P3[Pipe 3] & P4[Pipe 4]
    
    subgraph NhiênLiệu_Va_Nap [Hệ thống nạp & Kim phun]
        I1[Injector 1] -.-> P1
        I2[Injector 2] -.-> P2
        I3[Injector 3] -.-> P3
        I4[Injector 4] -.-> P4
        P1 --> MP1[Intake Port 1]
        P2 --> MP2[Intake Port 2]
        P3 --> MP3[Intake Port 3]
        P4 --> MP4[Intake Port 4]
    end

    subgraph KhoiDongCo [Khối sinh công 4 Xylanh]
        E1[Engine: Điều khiển trung tâm] --- C1[Cylinder 1]
        E1 --- C2[Cylinder 2]
        E1 --- C3[Cylinder 3]
        E1 --- C4[Cylinder 4]
        MP1 --> C1
        MP2 --> C2
        MP3 --> C3
        MP4 --> C4
    end

    subgraph HeThongThai [Đường xả & Đo kiểm]
        C1 --> EP1[Exhaust Port 1]
        C2 --> EP2[Exhaust Port 2]
        C3 --> EP3[Exhaust Port 3]
        C4 --> EP4[Exhaust Port 4]
        EP1 & EP2 & EP3 & EP4 --> CAT1[Catalytic Converter: Bộ xử lý khí thải]
        CAT1 --> SB2[System Boundary 2: Môi trường xả]
    end

Thông số kỹ thuật của động cơ đối tượng (Toyota 1NZ-FE)

  • Kiểu động cơ: 4 kỳ, 4 xylanh thẳng hàng, làm mát bằng dung dịch.
  • Đường kính xylanh (Bore - $D$): $75.0\text{ mm}$.
  • Hành trình piston (Stroke - $S$): $84.7\text{ mm}$.
  • Dung tích công tác ($V_h$): $1497\text{ cm}^3$ ($1.5\text{L}$).
  • Tỷ số nén ($\varepsilon$): $10.5:1$.
  • Chiều dài thanh truyền ($L_{rod}$): $140.5\text{ mm}$.
  • Công suất cực đại thiết kế: $80\text{ kW}$ tại $6000\text{ RPM}$.
  • Mô-men xoắn cực đại thiết kế: $141\text{ N}\cdot\text{m}$ tại $4200\text{ RPM}$.
  • Hệ thống phân phối khí: DOHC 16 van, điều khiển xú-páp biến thiên VVT-i.

Phương pháp luận nghiên cứu (Methodology)

   ┌────────────────────────────────────────────────────────┐
   │ BƯỚC 1: KHẢO SÁT HÌNH HỌC & ĐẶC TÍNH NHIỆT ĐỘNG 1NZ-FE  │
   └──────────────────────────┬─────────────────────────────┘
                              │
                              ▼
   ┌────────────────────────────────────────────────────────┐
   │ BƯỚC 2: THIẾT LẬP MÔ HÌNH TRÊN AVL WORKSPACE GUI       │
   │  - Khai báo Engine, Cylinders, Pipes, Restrictors...   │
   │  - Khai báo profile đóng/mở xú-páp (Intake/Exhaust)    │
   └──────────────────────────┬─────────────────────────────┘
                              │
                              ▼
   ┌────────────────────────────────────────────────────────┐
   │ BƯỚC 3: CẤU HÌNH NHIÊN LIỆU (BOOST GAS PROPERTIES)      │
   │  - E0 (RON 92) -> E5 -> E10 -> E15 -> E30 -> E50 -> E85│
   │  - Khai báo LHV, A/F stoichiometric, C-H-O ratio       │
   └──────────────────────────┬─────────────────────────────┘
                              │
                              ▼
   ┌────────────────────────────────────────────────────────┐
   │ BƯỚC 4: GIẢI HỆ PHƯƠNG TRÌNH VI PHÂN (RUNGE-KUTTA 4)   │
   │  - Bảo toàn năng lượng, Vibe 2-Zone, Woschni Heat Loss │
   │  - Tính toán động học phát thải NOx, CO, HC            │
   └──────────────────────────┬─────────────────────────────┘
                              │
                              ▼
   ┌────────────────────────────────────────────────────────┐
   │ BƯỚC 5: HẬU XỬ LÝ (IMPRESS CHART) & KIỂM CHỨNG SAI SỐ  │
   │  - So sánh đồ thị Ne, Me với số liệu Catalogue chuẩn   │
   │  - Xuất bảng dữ liệu đánh giá tối ưu hóa phát thải     │
   └────────────────────────────────────────────────────────┘
  • Kế hoạch thực hiện: Dự án được tiến hành trong 21 tuần (từ ngày 17/03/2022 đến ngày 09/08/2022), phân chia thành các mốc đánh giá mô hình định kỳ hàng tháng cùng giảng viên hướng dẫn.
  • Đảm bảo chất lượng (QA): Kiểm chứng tính hội tụ của chu trình nhiệt động qua các vòng lặp tính toán liên tiếp, đảm bảo độ chênh lệch áp suất và nhiệt độ giữa chu trình $k$ và $k+1$ đạt dưới $0.01%$.

Triển khai thực nghiệm và kết quả

Quy trình thiết lập mô hình và nền tảng thuật toán

1. Phương trình bảo toàn năng lượng trong buồng cháy (First Law of Thermodynamics)

Trạng thái nhiệt động của môi chất trong xylanh tại mỗi góc quay trục khuỷu ($\alpha$) được xác định qua phương trình vi phân:

$$\frac{d(m_c \cdot u)}{d\alpha} = -p_c \cdot \frac{dV}{d\alpha} + \frac{dQ_F}{d\alpha} - \sum \frac{dQ_W}{d\alpha} - h_{BB} \cdot \frac{dm_{BB}}{d\alpha} + \sum h_i \cdot \frac{dm_i}{d\alpha} - \sum h_e \cdot \frac{dm_e}{d\alpha}$$

Trong đó:

  • $m_c, u, p_c, V$: Khối lượng, nội năng riêng, áp suất và thể tích tức thời của môi chất trong xylanh.
  • $\frac{dQ_F}{d\alpha}$: Tốc độ tỏa nhiệt từ quá trình cháy của nhiên liệu.
  • $\frac{dQ_W}{d\alpha}$: Tốc độ truyền nhiệt tổn thất qua vách xylanh, nắp máy và đỉnh piston.
  • $h_{BB} \cdot \frac{dm_{BB}}{d\alpha}$: Tổn thất enthalpy do hiện tượng lọt khí qua khe hở xéc-măng (Blow-by).
  • $h_i \frac{dm_i}{d\alpha}, h_e \frac{dm_e}{d\alpha}$: Enthalpy mang vào/ra qua xú-páp nạp và xú-páp thải.

2. Mô hình quy luật cháy Vibe 2 vùng (Vibe 2-Zone Function)

Tốc độ cháy và phần trăm khối lượng nhiên liệu đã cháy ($x$) theo góc quay trục khuỷu ($\alpha$) được mô tả bởi phương trình:

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

  • $\alpha_0$: Thời điểm bắt đầu cháy (Start of Combustion - SOC).
  • $\Delta \alpha_c$: Thời gian cháy kéo dài (Combustion Duration), được hiệu chỉnh theo nồng độ ethanol ($\approx 40^\circ - 55^\circ$ góc quay trục khuỷu).
  • $m$: Tham số hình dáng đặc trưng Vibe (Vibe shape parameter, dao động từ $1.2$ đến $2.0$).

3. Mô hình động học phát thải $NO_x$ mở rộng (Extended Zeldovich Mechanism)

Cơ chế hình thành Nitric Oxide ($NO$) nhiệt độ cao được tính toán qua 3 phản ứng thuận nghịch:

$$N_2 + O \rightleftharpoons NO + N \quad (k_{f1}, k_{r1})$$ $$N + O_2 \rightleftharpoons NO + O \quad (k_{f2}, k_{r2})$$ $$N + OH \rightleftharpoons NO + H \quad (k_{f3}, k_{r3})$$

Tốc độ biến thiên nồng độ $NO$ phụ thuộc vào nhiệt độ vùng cháy tức thời ($T_{burned}$) và nồng độ oxy nguyên tử $[O]$.

// Cấu hình tham số môi chất & biên nạp trong AVL BOOST (Trích đoạn cấu hình phần tử Cylinder)
CYLINDER_PROPERTIES {
    Bore                = 0.075;      // [m]
    Stroke              = 0.0847;     // [m]
    ConRodLength        = 0.1405;     // [m]
    CompressionRatio    = 10.5;       // [-]
    CombustionModel     = VIBE_2_ZONE;
    HeatTransferModel   = WOSCHNI_1978;
    PistonTemp          = 550.0;      // [K]
    HeadTemp            = 520.0;      // [K]
    LinerTemp           = 450.0;      // [K]
}

FUEL_PROPERTIES_E15 {
    LowerHeatingValue   = 40.82e6;    // [J/kg]
    StoichAirFuelRatio  = 13.85;      // [kg_air / kg_fuel]
    Density             = 754.2;      // [kg/m^3]
    OxygenContentMass   = 0.052;      // 5.2% O2 by weight
}

Kết quả mô phỏng và kiểm chứng mô hình

Kiểm chứng độ chính xác mô hình cơ sở (Xăng khoáng RON 92 - E0)

Mô hình mô phỏng sau khi hiệu chỉnh tham số dòng chảy tại miệng xú-páp ($\mu\sigma$) và tổn thất ma sát được đối chiếu trực tiếp với đường đặc tính ngoài thực nghiệm của động cơ 1NZ-FE:

Tốc độ động cơ (RPM) Công suất mô phỏng (kW) Công suất thực nghiệm (kW) Sai số công suất (%) Mô-men mô phỏng (N·m) Mô-men thực nghiệm (N·m) Sai số mô-men (%)
1000 12.85 13.10 -1.91% 122.7 125.1 -1.92%
2000 27.90 28.50 -2.10% 133.2 136.0 -2.06%
3000 43.15 44.00 -1.93% 137.4 140.0 -1.86%
4000 58.60 59.50 -1.51% 139.9 141.0 -0.78%
5000 72.40 73.50 -1.50% 138.3 140.4 -1.50%
6000 79.10 80.00 -1.12% 125.9 127.3 -1.10%

Đánh giá: Sai số tối đa giữa mô phỏng và thực nghiệm chỉ đạt $2.1%$, hoàn toàn thỏa mãn tiêu chuẩn kiểm chứng mô hình nhiệt động lực học ĐCĐT (cho phép sai số dưới $5%$).

Ảnh hưởng của tỷ lệ pha trộn Ethanol đến tính năng kỹ thuật và phát thải (Tại tốc độ toàn tải 4000 RPM)

Loại nhiên liệu Công suất $N_e$ (kW) Mô-men xoắn $M_e$ (N·m) $BSFC$ ($g_e$) (g/kW·h) Phát thải $CO$ (g/kW·h) Phát thải $HC$ (g/kW·h) Phát thải $NO_x$ (g/kW·h)
E0 (RON 92) 58.60 139.90 248.5 42.10 1.85 8.20
E5 58.25 (-0.60%) 139.06 (-0.60%) 254.8 (+2.53%) 40.10 (-4.75%) 1.82 (-1.62%) 8.25 (+0.61%)
E10 57.80 (-1.36%) 137.99 (-1.36%) 262.1 (+5.47%) 38.20 (-9.26%) 1.77 (-4.32%) 8.32 (+1.46%)
E15 57.30 (-2.22%) 136.79 (-2.22%) 270.4 (+8.81%) 36.40 (-13.54%) 1.71 (-7.57%) 8.38 (+2.20%)
E30 55.40 (-5.46%) 132.25 (-5.46%) 295.2 (+18.79%) 34.10 (-19.00%) 1.59 (-14.05%) 8.45 (+3.05%)
E50 53.10 (-9.38%) 126.76 (-9.38%) 328.6 (+32.23%) 31.80 (-24.46%) 1.48 (-20.00%) 8.51 (+3.78%)
E85 50.68 (-13.51%) 120.99 (-13.51%) 370.8 (+49.21%) 29.11 (-30.86%) 1.38 (-25.41%) 8.58 (+4.63%) 
    TÍNH NĂNG ĐỘNG CƠ (4000 RPM)       │       PHÁT THẢI Ô NHIỄM (4000 RPM)
  100% ┌──────────────────────┐        │     100% ┌──────────────────────┐
       │████████ Công suất Ne │        │          │▒▒▒▒▒▒▒▒ Khí thải CO  │
   80% │██████░░ (-13.5%)     │        │      80% │▒▒▒▒▒▒░░ (-30.9%)     │
       │                      │        │          │                      │
   60% │                      │        │      60% │                      │
       │                      │        │          │                      │
   40% │                      │        │      40% │                      │
       │                      │        │          │                      │
    0% └──────────────────────┘        │       0% └──────────────────────┘
          E0   ───>   E85              │             E0   ───>   E85     
      [Suất tiêu hao ge tăng +49.2%]   │         [HC giảm -25.4%, NOx biến thiên nhẹ]

Đổi mới và đóng góp khoa học

  1. Ứng dụng công cụ mô phỏng 1D tiên tiến vào nghiên cứu năng lượng sạch: Dự án đã chứng minh khả năng số hóa hoàn toàn quá trình đánh giá nhiên liệu sinh học trên động cơ thương mại phổ biến tại Việt Nam, rút ngắn hơn $80%$ thời gian khảo sát so với phương pháp thử nghiệm cơ học.
  2. Khảo sát toàn diện dải tỷ lệ pha trộn rộng (E0 đến E85): Cung cấp bộ cơ sở dữ liệu chi tiết về đặc tính áp suất buồng cháy ($p_c$), chỉ số công suất, mô-men và hàm lượng phát thải theo từng cấp phối trộn, làm tiền đề cho việc xây dựng bản đồ nạp (Fuel Map) cho ECU xe linh hoạt nhiên liệu (FFV).
  3. Làm rõ cơ chế đối nghịch giữa hiệu suất năng lượng và giảm phát thải:
    • Năng lượng: Sự sụt giảm công suất ($13.5%$) và tăng suất tiêu hao ($49.2%$) ở E85 được giải thích trực tiếp qua mật độ năng lượng riêng và tỷ lệ nạp nhiên liệu cần thiết để duy trì hệ số dư lượng không khí tối ưu ($\lambda = 1$).
    • Môi trường: Hàm lượng oxy liên kết nội tại trong phân tử ethanol giúp tăng cường quá trình oxy hóa ở các vùng hỗn hợp giàu, triệt tiêu các mầm tạo muội than và giảm lượng $CO$ tới $30.86%$, $HC$ tới $25.41%$.
    • $NO_x$: Mức tăng phát thải $NO_x$ được kiềm chế ở mức rất thấp ($+4.63%$) do hiệu ứng làm mát khí nạp từ nhiệt ẩn hóa hơi cực lớn của ethanol ($840\text{ kJ/kg}$) đã kìm hãm nhiệt độ cháy đỉnh trong buồng đốt.

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

   ┌────────────────────────────────────────────────────────┐
   │         LỘ TRÌNH ỨNG DỤNG THỰC TIỄN CHO ĐỘNG CƠ         │
   ├──────────────────────────┬─────────────────────────────┤
   │ XĂNG SINH HỌC THẤP (E5 - E15)│ XĂNG SINH HỌC CAO (E30 - E85)│
   ├──────────────────────────┼─────────────────────────────┤
   │ - Tương thích 100% cơ khí│ - Thay thế ống dẫn, gioăng cao su│
   │ - Giữ nguyên bản đồ ECU  │ - Cảm biến cồn (Fuel Composition)│
   │ - Áp dụng đại trà tức thì│ - Remap góc đánh lửa & kim phun│
   └──────────────────────────┴─────────────────────────────┘

Kịch bản triển khai trong thực tế

  • Giai đoạn sử dụng trực tiếp (E5 - E15): Các phương tiện sử dụng động cơ xăng truyền thống (như Toyota Vios, Yaris) có thể vận hành ổn định ngay lập tức với xăng E5 và E10 mà không cần sửa đổi kết cấu hay hiệu chỉnh phần mềm điều khiển. Độ sụt giảm công suất là không đáng kể ($< 2.5%$), nằm trong ngưỡng dung sai vận hành của người lái.
  • Giai đoạn chuyển đổi xe linh hoạt nhiên liệu (FFV - E30 đến E85):
    • Nâng cấp phần cứng: Thay thế các đường ống nhiên liệu cao su tự nhiên bằng cao su tổng hợp chịu cồn (Fluorocarbon elastomers - FKM), nâng cấp bơm xăng và tăng lưu lượng kim phun thêm $30 - 40%$ để đáp ứng nhu cầu cấp nhiên liệu lớn hơn.
    • Tích hợp cảm biến thành phần nhiên liệu (Ethanol Content Sensor): Giúp ECU nhận biết tỷ lệ % ethanol tức thời, từ đó tự động hiệu chỉnh thời gian mở kim phun ($Pulse Width$) và tăng góc đánh lửa sớm để tận dụng triệt để chỉ số Octane cao ($RON > 108$) của ethanol, bù đắp hoàn toàn sự suy giảm công suất.

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

  • Kinh tế người dùng: Mặc dù suất tiêu hao nhiên liệu theo thể tích tăng khi sử dụng ethanol nồng độ cao, nhưng việc tận dụng nguồn cồn sinh học sản xuất từ phụ phẩm nông nghiệp (sắn lát, rỉ đường mía, bã ngô) với giá thành rẻ sẽ giúp cân bằng tổng chi phí vận hành trên mỗi $100\text{ km}$.
  • Lợi ích xã hội & môi trường: Giảm hàng triệu tấn khí độc hại $CO, HC$ phát tán vào môi trường đô thị, thúc đẩy phát triển chuỗi giá trị nông nghiệp bền vững và bảo đảm an ninh năng lượng quốc gia.

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

Hạn chế kỹ thuật của đề tài

  1. Mô hình 1 chiều lý tưởng hóa: Chưa thể hiện chi tiết sự phân lớp hòa khí cục bộ và quá trình phun sương dính ướt thành buồng cháy (Wall-wetting effect) như các phần mềm CFD 3D.
  2. Chưa khảo sát hiện tượng kích nổ thực tế: Chưa kích hoạt mô hình tự kích nổ (Auto-ignition kinetics) để tìm góc đánh lửa sớm giới hạn (Knock Limit Spark Advance - KLSA) nhằm tối ưu hóa công suất cho E85.
  3. Đánh giá ở trạng thái tĩnh: Chưa thực hiện mô phỏng chu trình lái xe thực tế (NEDC, WLTP) để đánh giá tiêu hao nhiên liệu và phát thải ở chế độ khởi động nguội và tăng tốc đột ngột.

Định hướng phát triển tiếp theo

  • Tích hợp mô phỏng đa chiều 1D-3D Co-simulation: Kết nối AVL BOOST với AVL FIRE (thông qua FIRE Link) để mô phỏng chi tiết dòng khí xoáy lốc và quá trình bắt cháy của tia lửa điện.
  • Nghiên cứu động cơ tăng áp phun trực tiếp (GDI Turbo): Kết hợp nhiên liệu E85 với động cơ tăng áp tỷ số nén cao để khai thác tối đa khả năng chống kích nổ, nâng cao hiệu suất nhiệt lên trên $40%$.
  • Xây dựng thuật toán ECU thích ứng: Thiết lập mô hình điều khiển vòng kín Lambda sử dụng mạng nơ-ron nhúng trong MATLAB DLL tích hợp vào AVL BOOST.

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

1. Sinh viên và Học viên ngành Kỹ thuật Ô tô / Động lực

  • Tiếp cận quy trình nghiên cứu mô phỏng chuẩn công nghiệp với phần mềm chuyên ngành hàng đầu thế giới (AVL BOOST).
  • Nắm vững mối liên hệ giữa nhiệt động lực học lý thuyết, hóa học nhiên liệu và vận hành thực tế của động cơ đốt trong.

2. Kỹ sư R&D và Lập trình viên ECU (Automotive Calibration Engineers)

  • Sử dụng bộ dữ liệu mô phỏng để tham chiếu khi lập bản đồ nạp (Fuel Mapping) và bản đồ đánh lửa (Ignition Mapping) cho các dòng động cơ chuyển đổi sang sử dụng nhiên liệu sinh học hoặc cồn E85.

3. Doanh nghiệp Vận tải và Nhà hoạch định chính sách Năng lượng

  • Có cơ sở khoa học định lượng vững chắc để đánh giá mức độ suy giảm tiêu hao nhiên liệu và tiềm năng cắt giảm phát thải độc hại khi ban hành các quy định bắt buộc về xăng sinh học.

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

1. Yêu cầu hệ thống phần cứng và phần mềm để triển khai mô phỏng AVL BOOST là gì?

Hệ thống cần máy tính chạy hệ điều hành Windows 10/11 64-bit hoặc Linux RedHat/Ubuntu LTS, cấu hình tối thiểu CPU 4 nhân (khuyến nghị Intel Core i7 hoặc AMD Ryzen 7), RAM tối thiểu 8GB (khuyến nghị 16GB), ổ cứng trống 10GB và đã cài đặt bộ phần mềm AVL BOOST (phiên bản 2019 đến 2022) có bản quyền license hợp lệ.

2. Tại sao công suất động cơ lại giảm khi sử dụng xăng sinh học có tỷ lệ ethanol cao?

Nguyên nhân chủ yếu do nhiệt trị thấp của ethanol ($26.8\text{ MJ/kg}$) thấp hơn đáng kể so với xăng khoáng ($44\text{ MJ/kg}$). Mặc dù lượng không khí lý thuyết để đốt cháy ethanol ít hơn ($9.0\text{ kg kk/kg nl}$ so với $14.7$ của xăng) và lượng nạp thể tích có thể tăng nhẹ nhờ nhiệt ẩn hóa hơi, nhưng tổng năng lượng nhiệt giải phóng trong một chu trình công tác vẫn thấp hơn so với xăng khoáng nếu không tăng tỷ số nén hoặc tăng góc đánh lửa sớm.

3. Phương tiện cũ sử dụng xăng E5 và E10 có bị ăn mòn kim loại hoặc hư hỏng gioăng cao su không?

Đối với các tỷ lệ ethanol thấp như E5 ($5%$) và E10 ($10%$), các nghiên cứu thực nghiệm và tiêu chuẩn sản xuất xe từ năm 1996 trở lại đây đều chứng minh hệ thống dẫn nhiên liệu hoàn toàn tương thích và không xảy ra hiện tượng biến tính gioăng phớt hay ăn mòn chi tiết kim loại. Hiện tượng ăn mòn chỉ xuất hiện khi tỷ lệ ethanol vượt quá $15 - 20%$ trên các vật liệu không tương thích.

4. Tại sao phát thải $CO$ và $HC$ lại giảm rõ rệt khi tăng nồng độ ethanol?

Phân tử ethanol ($C_2H_5OH$) chứa sẵn $34.7%$ oxy theo khối lượng. Lượng oxy liên kết này tham gia trực tiếp vào phản ứng cháy bên trong xylanh, đóng vai trò như một chất oxy hóa bổ trợ tại các vùng hòa khí cục bộ giàu nhiên liệu, giúp quá trình cháy diễn ra hoàn toàn hơn và chuyển hóa phần lớn $CO$ thành $CO_2$ cũng như đốt cháy triệt để các hydrocacbon chưa cháy ($HC$).

5. Làm thế nào để bù đắp sự suy giảm công suất khi sử dụng xăng E85 trong thực tế?

Để tận dụng tối đa xăng E85 và bù đắp suy giảm công suất, kỹ sư có thể áp dụng 3 giải pháp:

  1. Tăng góc đánh lửa sớm ($Advance Ignition Timing$) từ $4^\circ - 8^\circ$ góc quay trục khuỷu do ethanol có chỉ số Octane rất cao ($RON \approx 108$).
  2. Tăng áp suất phun và thời gian phun nhiên liệu để cung cấp đủ lưu lượng khối lượng.
  3. Trên các động cơ thiết kế mới, tăng tỷ số nén động cơ ($\varepsilon$ từ $10.5:1$ lên $12.5:1$ hoặc cao hơn) hoặc sử dụng hệ thống tăng áp Turbocharger để ép thêm lượng khí nạp vào xylanh.

Kết luận

Đồ án tốt nghiệp "Nghiên cứu ứng dụng AVL BOOST mô phỏng động cơ xăng - ethanol" đã xây dựng thành công mô hình nhiệt động lực học 1 chiều chính xác cho động cơ Toyota 1NZ-FE 1.5L. Kết quả nghiên cứu khẳng định tính khả thi vượt trội của công cụ mô phỏng số trong việc dự báo chu trình công tác và phát thải khí xả của động cơ đốt trong khi chuyển đổi sang sử dụng các nguồn nhiên liệu tái tạo.

Việc ứng dụng xăng sinh học mang lại giá trị môi trường to lớn với mức độ cắt giảm phát thải khí độc $CO$ lên đến $30.86%$ và $HC$ lên đến $25.41%$ khi sử dụng E85. Mặc dù công suất động cơ có sự suy giảm theo tỷ lệ thuận với hàm lượng cồn pha trộn, các giải pháp kỹ thuật về điều khiển điện tử và cải tiến hệ thống cung cấp nhiên liệu hoàn toàn có thể khắc phục triệt để nhược điểm này. Đây là cơ sở khoa học và công nghệ quan trọng đóng góp vào lộ trình chuyển dịch năng lượng xanh và phát triển phương tiện giao thông bền vững tại Việt Nam.