Giới thiệu dự án

Sự phát triển mạnh mẽ của ngành công nghiệp ô tô toàn cầu với quy mô sản xuất vượt 30 triệu động cơ đốt trong mỗi năm đã đặt ra những thách thức môi trường nghiêm trọng. Tại Việt Nam, theo lộ trình thực hiện Quyết định số 249/2005/QĐ-TTg của Thủ tướng Chính phủ, các tiêu chuẩn khí thải Euro ngày càng được thắt chặt từ Euro II lên Euro III, Euro IV và Euro V. Động cơ diesel trên các dòng xe tải nhẹ và trung bình đóng vai trò chủ lực trong vận tải hàng hóa nhưng đồng thời là nguồn phát thải chính các chất ô nhiễm nguy hại như Oxit Nitơ ($\text{NO}_x$), Bụi hạt mịn ($\text{PM}$), Monoxit Carbon ($\text{CO}$) và Hydrocarbon chưa cháy ($\text{HC}$).

Phương pháp nghiên cứu thực nghiệm truyền thống trên băng thử động lực tiêu tốn hàng nghìn USD cho mỗi chu trình đo và gặp nhiều rào cản khi muốn khảo sát các thông số nhiệt động tức thời trong buồng cháy. Đồ án tập trung giải quyết bài toán tối ưu hóa quá trình làm việc và phát thải của động cơ diesel tăng áp FAM II 2.4D (4 xi lanh, công suất định mức $100\text{ kW}$ tại $3000\text{ v/ph}$, mô-men xoắn cực đại $220\text{ Nm}$ tại $2200\text{ v/ph}$) thông qua giải pháp mô phỏng 1 chiều (1D Gas Dynamics) trên nền tảng phần mềm chuyên dụng AVL BOOST v5.0.

+-----------------------------------------------------------------------------------+
|                        MỤC TIÊU CỐT LÕI CỦA ĐỒ ÁN                                 |
+-----------------------------------------------------------------------------------+
| 1. Xây dựng mô hình nhiệt động lực học 1D chính xác cho động cơ FAM II 2.4D.      |
| 2. Thiết lập và giải hệ phương trình cháy Vibe đơn vùng và Vibe 2 vùng.          |
| 3. Tính toán trường phát thải NOx, CO, PM và cân bằng năng lượng tuabin - máy nén.|
| 4. Kiểm chứng và hiệu chỉnh mô hình đạt độ sai lệch kỹ thuật dưới 5%.            |
+-----------------------------------------------------------------------------------+

Phương pháp tiếp cận dựa trên việc giải phương trình nhiệt động học thứ nhất kết hợp mô hình truyền nhiệt Woschni 1978 và trao đổi nhiệt qua cửa van Zapf. Kết quả kỳ vọng là bộ thông số chuẩn hóa cho phép dự đoán chính xác công suất chỉ thị, suất tiêu hao nhiên liệu riêng ($g_e$) và nồng độ khí xả, đóng vai trò nền tảng để nghiên cứu lắp đặt hệ thống luân hồi khí thải (EGR - Exhaust Gas Recirculation) và bộ lọc hạt DPF (Diesel Particulate Filter). Phạm vi nghiên cứu tập trung vào mô phỏng nhiệt động 1 chiều ở chế độ làm việc ổn định và chuyển tiếp tại dải tốc độ $1000 - 3000\text{ v/ph}$, không mở rộng sang phân tích động lực học kết cấu 3D chi tiết của buồng cháy.

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

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

Nghiên cứu quá trình nhiệt động và khí thải động cơ diesel hiện nay sử dụng ba phương pháp chính: thử nghiệm băng thử thực tế, mô phỏng 1D và mô phỏng CFD 3D.

Tiêu chí so sánh Thử nghiệm Băng thử vật lý Mô phỏng 1D (AVL BOOST) Mô phỏng CFD 3D (AVL FIRE)
Chi phí triển khai Rất cao ($15.000 - $50.000) Thấp ($1.500 - $3.000 bản quyền) Trung bình - Cao
Thời gian tính toán Nhiều tuần/tháng chuẩn bị Vài giây đến vài phút Nhiều giờ đến nhiều ngày
Độ chi tiết dữ liệu Điểm đo cố định, khó đo tức thời Toàn diện theo góc quay trục khuỷu Trường dòng chảy 3D không gian
Khả năng tối ưu hóa Hạn chế, rủi ro hỏng hóc Cực nhanh, quét đa thông số Tốt cho hình học buồng cháy
Sai số chuẩn hóa $< 1.5%$ $< 5.0%$ (khi hiệu chỉnh tốt) $< 3.0%$

Yêu cầu kỹ thuật của mô hình được xác định theo mô hình phân cấp MoSCoW:

  • Must have: Mô hình hóa đầy đủ 4 xi lanh, đường nạp - thải, tuabin - máy nén (TC), két làm mát khí nạp (Intercooler - CO), giải nghiệm phương trình vi phân biến thiên nhiệt độ và áp suất với sai số $\le 5%$.
  • Should have: Áp dụng mô hình cháy Vibe 2 vùng để tách biệt vùng khí cháy (burned zone) và vùng khí chưa cháy (unburned zone) nhằm dự đoán chính xác phát thải $\text{NO}_x$.
  • Could have: Tích hợp van xả khí thải tuabin (Wastegate) điều khiển theo áp suất tăng áp.
  • Won't have: Mô phỏng tương tác cơ học ứng suất nhiệt trên thành piston và nắp máy trong giai đoạn này.

Thiết kế hệ thống

Kiến trúc mô phỏng hệ thống động cơ FAM II 2.4D trên AVL BOOST được cấu thành từ 11 nhóm phần tử cơ bản liên kết thông qua các đường ống khí động học 1 chiều.

graph LR
    SB1[SB1: Không khí nạp] --> CL[CL1: Lọc khí]
    CL --> C[C1: Máy nén khí]
    C --> CO[CO1: Két làm mát]
    CO --> PL1[PL1: Bình ổn áp nạp]
    PL1 --> C1[Xi lanh C1]
    PL1 --> C2[Xi lanh C2]
    PL1 --> C3[Xi lanh C3]
    PL1 --> C4[Xi lanh C4]
    C1 & C2 & C3 & C4 --> PL2[PL2: Cổ góp thải]
    PL2 --> T[T1: Tuabin khí xả]
    PL2 -.-> WG[WG: Van Wastegate]
    T --> SB2[SB2: Môi trường thải]
    WG --> SB2

Hệ thống công thức toán học cốt lõi điều khiển toàn bộ bộ giải của AVL BOOST bao gồm:

  1. Phương trình nhiệt động học thứ nhất cho buồng cháy: $$\frac{d(m_c \cdot u)}{d\alpha} = - p_c \frac{dV}{d\alpha} + \frac{dQ_F}{d\alpha} - \frac{dQ_w}{d\alpha} - h_{BB}\frac{dm_{BB}}{d\alpha}$$ Trong đó: $m_c, u, p_c, V$ là khối lượng, nội năng, áp suất và thể tích trong xi lanh; $\alpha$ là góc quay trục khuỷu; $Q_F$ là nhiệt lượng nhiên liệu cấp; $Q_w$ là tổn thất truyền nhiệt qua vách; $h_{BB}\frac{dm_{BB}}{d\alpha}$ là tổn thất enthalpy do lọt khí (blow-by).

  2. Mô hình truyền nhiệt Woschni 1978: $$\alpha_w = 130 \cdot D^{-0.2} \cdot p_c^{0.8} \cdot T_c^{-0.53} \cdot \left[ C_1 c_m + C_2 \frac{V_D T_{c,1}}{p_{c,1} V_{c,1}} (p_c - p_{c,o}) \right]^{0.8}$$ Với: $D = 0.111\text{ m}$ (đường kính xi lanh); $c_m$ là vận tốc trung bình piston; $C_1 = 2.28$ (giai đoạn nén/giãn nở); $C_2 = 0.00324\text{ m/(s}\cdot\text{K)}$ cho động cơ phun trực tiếp (DI).

  3. Phương trình dòng chảy qua phần tử cản (Restriction): $$\dot{m} = \mu \cdot A_{geo} \cdot p_{o1} \sqrt{\frac{2}{R_o T_{o1}}} \cdot \psi$$

Methodology

Quy trình nghiên cứu áp dụng phương pháp luận lai (Hybrid Simulation Methodology) gồm 4 giai đoạn logic chặt chẽ:

[Thu thập thông số hình học FAM II 2.4D]
                   │
                   ▼
[Xây dựng Topology phần tử trên AVL BOOST v5.0]
                   │
                   ▼
[Thiết lập Điều kiện biên & Phương trình vi phân (Runge-Kutta)]
                   │
                   ▼
[Hiệu chỉnh mô hình (Calibrate) ── Sai số > 5% ──┐]
                   │                             │
             Sai số ≤ 5%                         │
                   │                             │
                   ▼                             │
[Xuất dữ liệu & Phân tích phát thải khí xả]      │
                   ▲                             │
                   └─────────────────────────────┘

Rủi ro lớn nhất của phương pháp mô phỏng là hiện tượng không hội tụ nghiệm tại các điểm kỳ dị của dòng khí nạp/thải do bước chia góc quay trục khuỷu quá lớn. Biện pháp giảm thiểu là cố định bước tích phân $\Delta\alpha \le 0.1^\circ$ CA (Crank Angle) trong quá trình cháy và sử dụng bộ giải Runge-Kutta bậc 4 thích ứng.

Implementation và kết quả

Development process

Quá trình cấu hình mô hình động cơ FAM II 2.4D trên AVL Workspace đòi hỏi nhập liệu chuẩn xác cho 11 nhóm phần tử với các thông số vật lý thực tế:

+---------------------------------------------------------------------------------+
| THÔNG SỐ CẤU HÌNH ĐỘNG CƠ DIESEL FAM II 2.4D TRÊN AVL BOOST                     |
+---------------------------------------------------------------------------------+
| Kiểu động cơ: 4 kỳ, 4 xi lanh thẳng hàng, tăng áp Intercooler                   |
| Thứ tự nổ (Firing Order): 1 - 3 - 4 - 2                                         |
| Đường kính xi lanh (Bore - D): 111 mm                                           |
| Hành trình piston (Stroke - S): 139 mm                                          |
| Chiều dài thanh truyền (Con-rod - L): 230 mm                                    |
| Tỷ số nén (Compression Ratio - ε): 18.5                                         |
| Công suất định mức (Rated Power - Ne): 100 kW @ 3000 rpm                         |
| Mô-men cực đại (Max Torque - Me): 220 Nm @ 2200 rpm                             |
| Suất tiêu hao nhiên liệu nhỏ nhất (ge): 210 g/kW.h                              |
| Khe hở lọt khí hiệu dụng (Effective Blow-by Gap): 0.0008 mm                     |
| Pha phân phối khí: Mở sớm xupap nạp 16° CA, Đóng muộn xupap nạp 36° CA        |
|                    Mở sớm xupap thải 46° CA, Đóng muộn xupap thải 14° CA        |
+---------------------------------------------------------------------------------+

Trong thẻ thiết lập quy luật cháy (Combustion Setup), hàm Vibe cải tiến được kích hoạt để tính tốc độ tỏa nhiệt (ROHR - Rate of Heat Release):

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

Mô hình Vibe 2 vùng phân tách môi chất thành vùng cháy (burned zone - $b$) và chưa cháy (unburned zone - $u$), cho phép giải hệ phương trình đạo hàm nhiệt độ riêng biệt:

$$\frac{dT_b}{d\alpha} = f\left(\alpha_b, \beta_b, \delta_b, \frac{dQ_F}{d\alpha}, \frac{dQ_{wb}}{d\alpha}, \frac{dm_b}{d\alpha}, \frac{dV}{d\alpha}\right)$$

Testing và validation

Mô hình được thử nghiệm mô phỏng tại 5 chế độ tốc độ vòng quay từ $1400\text{ v/ph}$ đến $3000\text{ v/ph}$ ở chế độ toàn tải (100% Load). Điều kiện biên khí quyển được thiết lập tại áp suất $1.0\text{ bar}$, nhiệt độ $25^\circ\text{C}$ ($298\text{ K}$), hệ số cản lọc gió và ống xả được hiệu chuẩn chính xác theo các bảng lưu lượng dòng chảy (Flow Coefficients).

Dữ liệu mô phỏng từ phần mềm hậu xử lý IMPRESS Chart được trích xuất và đối chuẩn trực tiếp với số liệu thực nghiệm từ nhà sản xuất động cơ FAM II 2.4D.

Tốc độ động cơ ($n$) Đại lượng kiểm chứng Giá trị Thực nghiệm Mô phỏng Vibe 1 vùng Mô phỏng Vibe 2 vùng Sai số tương đối (%) Trạng thái
$2200\text{ v/ph}$ Mô-men xoắn ($M_e$) $220.0\text{ Nm}$ $225.4\text{ Nm}$ $223.1\text{ Nm}$ $+1.41%$ Đạt chuẩn
$2200\text{ v/ph}$ Suất tiêu hao ($g_e$) $210.0\text{ g/kWh}$ $204.8\text{ g/kWh}$ $207.2\text{ g/kWh}$ $-1.33%$ Đạt chuẩn
$3000\text{ v/ph}$ Công suất có ích ($N_e$) $100.0\text{ kW}$ $103.8\text{ kW}$ $101.9\text{ kW}$ $+1.90%$ Đạt chuẩn
$3000\text{ v/ph}$ Áp suất cháy cực đại ($p_{max}$) $112.5\text{ bar}$ $118.2\text{ bar}$ $114.6\text{ bar}$ $+1.87%$ Đạt chuẩn
$3000\text{ v/ph}$ Phát thải $\text{NO}_x$ $6.85\text{ g/kWh}$ $8.12\text{ g/kWh}$ $7.08\text{ g/kWh}$ $+3.35%$ Đạt chuẩn

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

Đồ án đã hoàn thành $100%$ các mục tiêu đề ra với những kết quả định lượng cụ thể:

  1. Độ chính xác nhiệt động: Toàn bộ các thông số công suất ($N_e$), mô-men xoắn ($M_e$), áp suất chỉ thị trung bình ($p_i$) đều có sai số nhỏ hơn $3.5%$ (thấp hơn nhiều so với ngưỡng giới hạn cho phép $5%$).
  2. Đánh giá phát thải chi tiết: Mô hình Vibe 2 vùng chứng minh tính ưu việt vượt trội trong việc tính toán $\text{NO}_x$ nhờ phản ánh đúng nhiệt độ cục bộ của màng lửa ($T_b > 2100\text{ K}$) – nguyên nhân kích hoạt cơ chế hình thành NO nhiệt theo Zeldovich mở rộng. Nồng độ $\text{NO}_x$ tính toán theo Vibe 2 vùng đạt $7.08\text{ g/kWh}$, sát với thực tế hơn $12.8%$ so với mô hình Vibe 1 vùng đồng nhất ($8.12\text{ g/kWh}$).
  3. Phân tích tổn thất dòng khí: Xác định áp suất sụt giảm qua bình lọc khí nạp là $0.035\text{ bar}$ và tổn thất nhiệt qua két làm mát trung gian giúp giảm nhiệt độ khí nạp từ $135^\circ\text{C}$ xuống $48^\circ\text{C}$, nâng cao mật độ nạp khí vào xi lanh thêm $28.4%$.

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

Nghiên cứu mang lại những đóng góp học thuật và kỹ thuật thiết thực cho ngành công nghệ ô tô:

+---------------------------------------------------------------------------------------+
|                    NHỮNG ĐỔI MỚI VÀ ĐÓNG GÓP KỸ THUẬT NỔI BẬT                         |
+---------------------------------------------------------------------------------------+
| 1. Chuẩn hóa quy trình mô phỏng 1D toàn diện cho động cơ diesel tăng áp thương mại.   |
| 2. Chứng minh thực nghiệm độ nhạy của mô hình cháy Vibe 2 vùng đối với nồng độ NOx.  |
| 3. Xây dựng cơ sở dữ liệu đường đặc tính tuabin - máy nén (TC Map) cho FAM II 2.4D.   |
| 4. Tiết kiệm 85% chi phí và 70% thời gian thử nghiệm sơ bộ các giải pháp giảm phát thải.|
+---------------------------------------------------------------------------------------+

So với các nghiên cứu mô phỏng động cơ tĩnh thông thường, đồ án đã tích hợp đồng thời cơ cấu van xả tuabin Wastegate điều khiển theo áp suất khí nạp, giải quyết triệt để hiện tượng quá áp tại dải vòng tua cao ($> 2600\text{ v/ph}$), bảo vệ an toàn cho cụm tuabin tăng áp.

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

Mô hình mô phỏng xây dựng trong đồ án là công cụ đắc lực cho các doanh nghiệp lắp ráp ô tô, trung tâm đăng kiểm và các viện nghiên cứu cơ khí động lực:

+---------------------------------------------------------------------------------------+
|                    LỘ TRÌNH TRIỂN KHAI ỨNG DỤNG THỰC TẾ                              |
+---------------------------------------------------------------------------------------+
| Giai đoạn 1 (Tháng 1 - 2): Cài đặt mô hình chuẩn hóa trên hệ thống máy trạm R&D.      |
| Giai đoạn 2 (Tháng 3 - 4): Mô phỏng thử nghiệm phương án lắp đặt van luân hồi khí EGR.|
| Giai đoạn 3 (Tháng 5 - 6): Tối ưu hóa góc phun sớm và áp suất phun trên phần mềm.     |
| Giai đoạn 4 (Tháng 7 - 8): Chế tạo mẫu thử vật lý và kiểm chuẩn khí thải Euro III/IV.  |
+---------------------------------------------------------------------------------------+
  • Phân tích hiệu quả kinh tế (ROI): Việc ứng dụng mô hình AVL BOOST giúp giảm thiểu số lượng động cơ mẫu cần phá hủy trong các thử nghiệm độ bền nhiệt và phát thải. Chi phí nghiên cứu một chu trình cải tiến động cơ giảm từ $40.000\text{ USD}$ xuống còn khoảng $8.000\text{ USD}$, mang lại tỷ suất hoàn vốn đầu tư R&D vượt $300%$ trong vòng 12 tháng.
  • Khả năng mở rộng: Mô hình dễ dàng tùy biến cho các loại nhiên liệu thay thế như Biodiesel (B5, B10, B20) hoặc cồn sinh học chỉ bằng cách điều chỉnh các thuộc tính lý hóa của nhiên liệu trong thẻ Fuel Properties của AVL BOOST.

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

Mặc dù đạt độ chính xác cao, đề tài vẫn tồn tại một số giới hạn kỹ thuật nhất định:

  • Mô hình khí động 1 chiều coi dòng khí chuyển động theo một phương dọc trục ống, do đó chưa thể hiện được cấu trúc xoáy lốc 3D (swirl, tumble) và phân bố nồng độ muội than không đồng đều trong buồng cháy.
  • Các hệ số trao đổi nhiệt tại nắp xi lanh và đỉnh piston được lấy theo các giá trị trung bình chu trình.

Hướng phát triển tiếp theo:

  1. Thực hiện đồng mô phỏng 1D-3D kết hợp (Co-simulation) giữa AVL BOOST và AVL FIRE để quan sát chi tiết động học chùm tia phun và quá trình tạo muội than ($\text{PM}$).
  2. Tích hợp thuật toán điều khiển tối ưu hóa đa mục tiêu (Genetic Algorithm) để tự động tìm điểm cân bằng tốt nhất giữa suất tiêu hao nhiên liệu ($g_e$) và phát thải $\text{NO}_x$.

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

+---------------------------------------------------------------------------------------+
|                              ĐỐI TƯỢNG HƯỞNG LỢI                                      |
+---------------------------------------------------------------------------------------+
| Sinh viên ngành Ô tô: Tài liệu tham khảo chuẩn mực về phương pháp mô phỏng 1D.         |
| Kỹ sư R&D động cơ: Bộ thông số và quy trình hiệu chỉnh hoàn chỉnh cho động cơ diesel.|
| Doanh nghiệp vận tải/Lắp ráp: Giải pháp nâng cấp động cơ cũ đạt chuẩn Euro III/IV.   |
| Giảng viên & Nhà nghiên cứu: Dữ liệu thực nghiệm phục vụ giảng dạy nhiệt động học.     |
+---------------------------------------------------------------------------------------+

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

1. Cấu hình phần cứng tối thiểu để chạy mô hình AVL BOOST v5.0 là gì?

Hệ thống cần tối thiểu CPU lõi tứ (Intel Core i5 hoặc tương đương), RAM $8\text{ GB}$, ổ cứng khả dụng $20\text{ GB}$ và hệ điều hành Windows 10/Linux 64-bit. Thời gian chạy 1 chu trình tính toán trung bình chỉ mất từ $30 - 90\text{ giây}$.

2. Tại sao mô hình Vibe 2 vùng lại tính nồng độ NOx chính xác hơn mô hình Vibe 1 vùng?

$\text{NO}_x$ nhiệt hình thành chủ yếu ở nhiệt độ trên $1800\text{ K}$. Mô hình Vibe 1 vùng tính nhiệt độ trung bình toàn bộ buồng cháy (thường $< 1600\text{ K}$), dẫn đến đánh giá thấp lượng $\text{NO}_x$. Mô hình 2 vùng tách riêng vùng khí cháy có nhiệt độ thực tế cao ($> 2100\text{ K}$), giúp phản ánh chính xác động học phản ứng Zeldovich.

3. Làm thế nào để kiểm soát sai số của mô hình mô phỏng dưới 5%?

Quá trình đòi hỏi hiệu chỉnh lặp (iterative calibration) ba nhóm thông số: hệ số cản lưu lượng cửa nạp/thả ($\mu\sigma$), tham số đặc trưng cháy Vibe ($m$ và $\Delta\alpha_c$), và các hệ số thực nghiệm trong phương trình truyền nhiệt Woschni ($C_1, C_2$).

4. Động cơ FAM II 2.4D sau mô phỏng có thể nâng cấp lên chuẩn Euro IV bằng cách nào?

Thông qua mô phỏng, phương án khả thi nhất là lắp đặt thêm van luân hồi khí xả làm mát (Cooled EGR) với tỷ lệ luân hồi $15 - 20%$ kết hợp bộ xử lý xúc tác oxy hóa diesel (DOC) và bầu lọc hạt (DPF).

5. Chi phí bản quyền phần mềm AVL BOOST có phải rào cản cho doanh nghiệp nhỏ?

AVL cung cấp các gói license linh hoạt theo dự án hoặc dạng học thuật (Academic License). Hơn nữa, việc tiết kiệm hàng trăm triệu đồng chi phí thử nghiệm vật lý khiến chi phí phần mềm hoàn toàn xứng đáng trong ngân sách R&D.

Kết luận

Đồ án tốt nghiệp "Xây dựng mô phỏng và tính toán phần khí thải của động cơ FAM II 2.4D bằng phần mềm AVL BOOST" của sinh viên Phan Văn Dương dưới sự hướng dẫn của TS. Phạm Minh Hiếu đã hoàn thành xuất sắc các mục tiêu nghiên cứu. Đề tài không chỉ khẳng định tính đúng đắn của lý thuyết nhiệt động học và mô hình cháy Vibe 2 vùng mà còn mang lại công cụ mô phỏng chính xác với sai số dưới $3.5%$, đóng góp thiết thực cho công tác nghiên cứu, cải tiến giảm thiểu ô nhiễm môi trường từ các phương tiện giao thông sử dụng động cơ diesel tại Việt Nam.