Giới thiệu dự án
Sự gia tăng nhanh chóng của phương tiện giao thông cá nhân, đặc biệt là xe mô tô và xe gắn máy hai bánh tại các đô thị lớn ở các quốc gia đang phát triển, đặt ra thách thức nghiêm trọng về ô nhiễm không khí. Theo các báo cáo môi trường giao thông đô thị, khí thải xe máy chiếm hơn 60% tổng lượng phát thải hydrocacbon (HC) và carbon monoxide (CO). Để kiểm soát tình trạng này, các quy chuẩn kiểm định phát thải nghiêm ngặt như Euro 3, Euro 4 theo tiêu chuẩn quốc gia TCVN 9726:2013 và chu trình thử nghiệm toàn cầu WMTC (World Harmonized Motorcycle Test Cycle) đã được ban hành, đòi hỏi các phòng thí nghiệm phải đo đạc chính xác nồng độ phát thải tức thời trên băng thử con lăn (Chassis Dynamometer).
Hệ thống lấy mẫu với thể tích không đổi (Constant Volume Sampling - CVS) là phương pháp quy chuẩn quốc tế để đo nồng độ và khối lượng phát thải. Tuy nhiên, một bài toán kỹ thuật lớn phát sinh trong quá trình thử nghiệm liên tục: tín hiệu nồng độ khí thải ghi nhận tại tủ phân tích CEBII (Combustion Emission Bench II) luôn bị trễ thời gian (time delay) và san phẳng tín hiệu (signal smoothing/dispersion) so với thời điểm xả thực tế tại cửa xả động cơ. Sự sai lệch này xuất phát từ quá trình vận chuyển dọc theo đường ống dẫn, sự hòa trộn dòng chảy rối tại ống chữ T (Mixing-T), sự giãn nở dòng khí trong bình chứa trung gian và quán tính phản hồi của cảm biến. Nếu không có thuật toán hiệu chỉnh độ trễ, việc đồng bộ giữa vận tốc xe theo chu trình lái và nồng độ khí thải tức thời sẽ bị sai lệch nghiêm trọng.
Mục tiêu cụ thể của dự án nghiên cứu bao gồm:
- Xây dựng mô hình hình học 3D chi tiết của hệ thống CVS tại Phòng thí nghiệm Động cơ đốt trong (Bộ môn Hệ thống Động lực Ô tô – Trường Cơ khí, Đại học Bách khoa Hà Nội).
- Thiết lập mô hình mô phỏng động lực học chất lỏng tính toán (Computational Fluid Dynamics - CFD) bằng phần mềm ANSYS Fluent để giải bài toán dòng chảy rối nhiều thành phần (Species Transport) không ổn định (transient).
- Khảo sát phân bố trường vận tốc, nhiệt độ và nồng độ khí thải theo các hệ số pha loãng ($DF = 4, 8, 12, 16, 20$) dưới hai chế độ lưu lượng ống Venturi ($1.0\text{ m}^3/\text{phút}$ và $2.5\text{ m}^3/\text{phút}$).
- Định lượng thời gian trễ của dòng khí từ đầu ống xả đến mặt cắt lấy mẫu $xOy$ và xây dựng hàm đa thức nội suy Lagrange phục vụ hiệu chỉnh tín hiệu đo liên tục.
Phương pháp giải pháp được lựa chọn là ứng dụng mô phỏng 3D CFD trên nền tảng ANSYS Fluent sử dụng mô hình rối $k-\varepsilon$ tiêu chuẩn và thuật toán phân tách (Segregated Solver). Phương pháp này cho phép tái hiện chi tiết cơ chế động lực học dòng khí mà các phương pháp đo thực nghiệm cục bộ không thể quan sát toàn diện. Kết quả kỳ vọng đạt độ chính xác hội tụ sai số thặng dư dưới $10^{-6}$, cung cấp bảng thông số thời gian trễ chính xác tới phần mười giây (ví dụ $7.9\text{ s}$ tại $DF=4$, $Q=1\text{ m}^3/\text{phút}$).
Phạm vi nghiên cứu giới hạn trên hệ thống CVS thực tế tại phòng thí nghiệm C15 ĐHBK Hà Nội với các thông số hình học cố định, môi chất được xem là hỗn hợp khí lý tưởng và bỏ qua quá trình bức xạ nhiệt cũng như các phản ứng hóa học sau ống xả.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trong thử nghiệm phát thải động cơ, các phương pháp đo lường hiện nay tồn tại những ưu và nhược điểm rõ rệt khi so sánh với giải pháp mô phỏng hiệu chỉnh CVS:
| Phương pháp đo |
Ưu điểm |
Nhược điểm |
Độ chính xác thời gian thực |
| Đo trực tiếp khí thô (Raw Gas Analysis) |
Đáp ứng nhanh, không cần hệ thống hòa trộn pha loãng lớn. |
Dễ ngưng tụ hơi nước làm hỏng cảm biến; tắc nghẽn muội than; áp suất dao động mạnh. |
Trung bình - Kém (do sai số áp suất và độ ẩm) |
| Lấy mẫu túi khí tĩnh CVS (Bag Sampling CVS) |
Cho kết quả trung bình tổng thể toàn chu trình chuẩn xác, triệt tiêu sai số tức thời. |
Không cung cấp được dữ liệu phát thải tức thời (Modal Analysis) theo từng pha tăng/giảm tốc. |
Không áp dụng cho đo tức thời |
| Đo liên tục có hiệu chỉnh CFD (CFD-Calibrated CVS) |
Bảo vệ thiết bị đo, phân tích chính xác nồng độ tức thời theo từng giây của chu trình WMTC. |
Yêu cầu mô hình hóa dòng khí phức tạp và thuật toán khử trễ tín hiệu chuyên sâu. |
Cao (sai số thời gian $< 0.2\text{ s}$) |
Phân tích yêu cầu hệ thống theo mô hình MoSCoW:
- Must have: Mô hình 3D chính xác kích thước thực của ống mềm inox, ống pha loãng chữ T, bình chứa và ống Venturi; thiết lập điều kiện biên vận tốc, nhiệt độ chuẩn theo hệ số pha loãng $DF$; mô phỏng trường dòng chảy không ổn định (transient).
- Should have: Mô hình hóa tương tác dòng rối $k-\varepsilon$; giải phương trình truyền pha và thành phần khí thải (Species Transport); trích xuất mặt cắt kiểm định $xOy$.
- Could have: Xây dựng thuật toán đa thức nội suy Lagrange tự động hóa tính toán thời gian trễ theo lưu lượng.
- Won't have: Mô hình hóa phản ứng cháy thứ cấp trong đường ống và sự ngưng tụ sâu của hơi nước.
Thách thức kỹ thuật lớn nhất là hiện tượng sụt giảm áp suất cục bộ và dòng xoáy phân ly mạnh tại góc ngoặt của ống hòa trộn chữ T và sự giãn nở dòng đột ngột trong bình chứa dung tích lớn ($\varnothing 300\text{ mm} \times 600\text{ mm}$), gây ra hiện tượng lưu giữ khí và kéo dài dải thời gian trễ.
Thiết kế hệ thống
Kiến trúc vật lý của hệ thống thử nghiệm bao gồm: Băng thử động lực học Chassis Dynamometer con lăn $20\text{ inch}$ $\to$ Khớp nối mềm inox $\varnothing 34\text{ mm}$ ($L = 770\text{ mm}$) $\to$ Ống nạp khí sạch pha loãng $\varnothing 150\text{ mm}$ $\to$ Ống hòa trộn pha loãng $\varnothing 100\text{ mm}$ ($L = 1815\text{ mm}$) $\to$ Bình chứa cân bằng $\varnothing 300\text{ mm}$ ($H = 600\text{ mm}$) $\to$ Mặt cắt lấy mẫu $xOy$ $\to$ Ống Venturi kiểm soát lưu lượng $\to$ Quạt hút lưu lượng $10\text{ m}^3/\text{phút}$ $\to$ Tủ phân tích khí xả CEBII.
Hệ thống phân tích CEBII sử dụng các module đo chuyên dụng:
- CO, $\text{CO}_2$: Phương pháp hấp thụ hồng ngoại không tán xạ (NDIR - Non-Dispersive Infrared).
- $\text{HC}$: Phương pháp ion hóa ngọn lửa (FID - Flame Ionization Detector).
- $\text{NO}_x$: Phương pháp phát quang hóa học (CLD - Chemiluminescence Detector).
+-------------------------------------------------------------------------+
| HỆ THỐNG THỬ NGHIỆM KHÍ THẢI CVS |
+-------------------------------------------------------------------------+
[ Xe máy trên Bệ Dyno ]
|
v (Khí xả thô: 300°C)
+----------------------+ +--------------------+
| Ống nối mềm Inox | | Không khí sạch |
| Ø34mm, L=770mm | | Ø150mm, 28°C |
+----------------------+ +--------------------+
\ /
\ /
v v
+--------------------------------+
| Ống hòa trộn chữ T (Mixing-T) |
+--------------------------------+
|
v
+--------------------------------+
| Ống pha loãng (Dilution Tube)|
| Ø100mm, L=1815mm |
+--------------------------------+
|
v
+--------------------------------+
| Bình chứa (Mixing Chamber) |
| Ø300mm, H=600mm |
+--------------------------------+
|
v
+--------------------------------+
| Mặt cắt lấy mẫu xOy | ----> [ Tủ phân tích CEBII ]
+--------------------------------+ (NDIR, FID, CLD)
|
v
+--------------------------------+
| Ống chuẩn Venturi |
| (1.0 hoặc 2.5 m3/phút) |
+--------------------------------+
|
v
+--------------------------------+
| Quạt hút (10 m3/phút) |
+--------------------------------+
Công nghệ và công cụ sử dụng:
- Mô hình CAD 3D: AutoCAD 2022.
- Bộ tiền xử lý & Chia lưới: ANSYS Meshing (Lưới phi cấu trúc Tetrahedral).
- Bộ giải CFD: ANSYS Fluent v2021/v2022 (Pressure-Based Solver, Second-Order Implicit Time Integration, Second-Order Upwind Discretization).
- Hậu xử lý: ANSYS CFD-Post.
Cơ sở toán học và hệ phương trình điều khiển:
-
Phương trình bảo toàn khối lượng (liên tục):
$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \vec{v}) = S_m$$
-
Phương trình bảo toàn động lượng (Navier-Stokes):
$$\frac{\partial}{\partial t}(\rho \vec{v}) + \nabla \cdot (\rho \vec{v} \vec{v}) = -\nabla p + \nabla \cdot (\bar{\bar{\tau}}) + \rho \vec{g} + \vec{F}$$
-
Phương trình năng lượng:
$$\frac{\partial}{\partial t}(\rho E) + \nabla \cdot (\vec{v}(\rho E + p)) = \nabla \cdot \left( k_{\text{eff}} \nabla T - \sum_j h_j \vec{J}j + (\bar{\bar{\tau}}{\text{eff}} \cdot \vec{v}) \right) + S_h$$
-
Phương trình bảo toàn thành phần (Species Transport):
$$\frac{\partial}{\partial t}(\rho Y_i) + \nabla \cdot (\rho \vec{v} Y_i) = -\nabla \cdot \vec{J}_i + R_i + S_i$$
-
Lưu lượng qua họng Venturi:
$$Q_v = K_v \cdot \frac{P_v}{\sqrt{T_v}}$$
Trong đó: $K_v$ là hệ số hiệu chuẩn ống Venturi; $P_v, T_v$ lần lượt là áp suất tuyệt đối và nhiệt độ tuyệt đối phía trước họng đo.
-
Hệ số pha loãng (Dilution Factor - $DF$):
$$DF = \frac{13.4}{C_{\text{CO}2} + C{\text{CO}} + C_{\text{HC}}}$$
Methodology
Phương pháp luận nghiên cứu tuân thủ quy trình 4 giai đoạn chuẩn mực kỹ thuật:
[Khảo sát & Đo thực nghiệm] --> [Mô hình hóa CAD/Mesh 3D] --> [Mô phỏng ANSYS Fluent] --> [Hồi quy toán học Lagrange]
- Giai đoạn 1: Thu thập số liệu thực nghiệm: Đo đạc kích thước thực tế hệ thống CVS phòng C15; ghi nhận chu trình thử nghiệm TCVN 9726:2013 (chu trình 1200 giây gồm 6 pha, lấy mẫu 4 pha sau từ giây thứ 430, vận tốc tối đa $50\text{ km/h}$).
- Giai đoạn 2: Tiền xử lý (Pre-processing): Xây dựng hình học 3D, thiết lập biên lưới Tetrahedral mịn tại các vùng chuyển tiếp dòng chảy.
- Giai đoạn 3: Tính toán số (Numerical Simulation): Áp dụng chiến lược hội tụ 3 bước: chạy chế độ dừng (steady-state) độ chính xác bậc 1 để khởi tạo trường dòng, sau đó chuyển sang dòng chảy không dừng (transient) bậc 2 với tiêu chuẩn hội tụ $10^{-6}$.
- Giai đoạn 4: Hậu xử lý & Hồi quy dữ liệu: Trích xuất trường nồng độ theo thời gian tại mặt cắt $xOy$, xác định độ trễ và thiết lập phương trình đa thức Lagrange.
Quản trị rủi ro: Hiện tượng phân kỳ nghiệm do số Courant lớn tại các hẹp dòng được kiểm soát bằng cách cố định bước thời gian $\Delta t = 0.1\text{ s}$ và giới hạn số vòng lặp tối đa 20 - 600 iterations/time-step.
Implementation và kết quả
Development process
Quá trình triển khai mô phỏng được thực hiện trên ANSYS Fluent với việc thiết lập các đặc tính vật lý và điều kiện biên cụ thể:
- Môi chất không khí (Air): $\rho = 1.225\text{ kg/m}^3$, $C_p = 1006.43\text{ J/kg}\cdot\text{K}$, hệ số dẫn nhiệt $k = 0.0242\text{ W/m}\cdot\text{K}$, độ nhớt động lực học $\mu = 1.7894 \times 10^{-5}\text{ kg/m}\cdot\text{s}$.
- Khí xả động cơ (Emission Gas): $\rho = 1.3679\text{ kg/m}^3$, $C_p = 1006.43\text{ J/kg}\cdot\text{K}$, $k = 0.0242\text{ W/m}\cdot\text{K}$, $\mu = 1.9738 \times 10^{-5}\text{ kg/m}\cdot\text{s}$.
- Điều kiện biên nhiệt độ: Đầu vào khí xả (Inlet 1) $T = 300^\circ\text{C}$ ($573.15\text{ K}$); Đầu vào không khí sạch (Inlet 2) $T = 28^\circ\text{C}$ ($301.15\text{ K}$); Đầu ra (Outlet) $T = 100^\circ\text{C}$ ($373.15\text{ K}$).
- Diện tích cửa vào khí xả: $A_{\text{inlet_emission}} = 0.00089759\text{ m}^2$ ($\varnothing 34\text{ mm}$).
Cấu hình giải thuật ANSYS Fluent được tự động hóa qua script thiết lập:
; === ANSYS FLUENT SETUP SCRIPT FOR CVS SIMULATION ===
/define/models/solver/pressure-based yes
/define/models/unsteady-1st-order? no
/define/models/unsteady-2nd-order? yes
/define/models/energy? yes
/define/models/viscous/ke-standard? yes
/define/models/multiphase/model/species-transport yes
; Define boundary conditions for DF = 4 (Q_outlet = 1 m3/min)
/define/boundary-conditions/velocity-inlet inlet_emission no no yes 4.643 no 0 no 573.15
/define/boundary-conditions/velocity-inlet inlet_air no no yes 0.707 no 0 no 301.15
/define/boundary-conditions/velocity-inlet outlet no no yes 2.301 no 0 no 373.15
; Solution Methods Configuration
/solve/set/discretization-scheme/pressure 12 ; PRESTO!
/solve/set/discretization-scheme/momentum 1 ; Second-Order Upwind
/solve/set/discretization-scheme/temperature 1 ; Second-Order Upwind
/solve/set/discretization-scheme/k 1 ; Second-Order Upwind
/solve/set/discretization-scheme/epsilon 1 ; Second-Order Upwind
; Transient Execution Parameters
/solve/set/time-step 0.1
/solve/dual-time-iterate 300 20
Testing và validation
Điều kiện biên vận tốc đầu vào được tính toán chặt chẽ theo công thức thủy lực $v = \frac{Q}{A}$ dựa trên lưu lượng pha loãng và hệ số $DF$:
Thông số vận tốc biên với lưu lượng Venturi $Q_{\text{outlet}} = 1.0\text{ m}^3/\text{phút}$ ($0.01667\text{ m}^3/\text{s}$):
| Hệ số pha loãng ($DF$) |
Lưu lượng khí xả ($Q_{\text{emission}}$) $[\text{m}^3/\text{phút}]$ |
Vận tốc $v_{\text{inlet_emission}}$ $[\text{m/s}]$ |
Vận tốc $v_{\text{inlet_air}}$ $[\text{m/s}]$ |
Vận tốc $v_{\text{outlet}}$ $[\text{m/s}]$ |
| 4 |
0.250 |
4.643 |
0.707 |
2.301 |
| 8 |
0.125 |
2.321 |
0.825 |
2.301 |
| 12 |
0.083 |
1.548 |
0.864 |
2.301 |
| 16 |
0.063 |
1.161 |
0.884 |
2.301 |
| 20 |
0.050 |
0.929 |
0.895 |
2.301 |
Thông số vận tốc biên với lưu lượng Venturi $Q_{\text{outlet}} = 2.5\text{ m}^3/\text{phút}$ ($0.04167\text{ m}^3/\text{s}$):
| Hệ số pha loãng ($DF$) |
Lưu lượng khí xả ($Q_{\text{emission}}$) $[\text{m}^3/\text{phút}]$ |
Vận tốc $v_{\text{inlet_emission}}$ $[\text{m/s}]$ |
Vận tốc $v_{\text{inlet_air}}$ $[\text{m/s}]$ |
Vận tốc $v_{\text{outlet}}$ $[\text{m/s}]$ |
| 4 |
0.625 |
11.607 |
1.768 |
5.753 |
| 8 |
0.313 |
5.803 |
2.063 |
5.753 |
| 12 |
0.208 |
3.869 |
2.161 |
5.753 |
| 16 |
0.156 |
2.902 |
2.210 |
5.753 |
| 20 |
0.125 |
2.321 |
2.239 |
5.753 |
Kết quả đạt được
- Trường phân bố nhiệt độ: Khí thải ở nhiệt độ cao $300^\circ\text{C}$ khi đi qua ống mềm inox và tiếp xúc với dòng khí sạch ($28^\circ\text{C}$) tại ống chữ T được hòa trộn nhanh chóng. Khi đến mặt cắt lấy mẫu $xOy$, nhiệt độ hỗn hợp khí giảm mạnh xuống còn $28^\circ\text{C} - 30^\circ\text{C}$, tiệm cận nhiệt độ môi trường, đảm bảo điều kiện làm việc an toàn cho cảm biến NDIR/FID và ngăn ngừa hiện tượng ngưng tụ hơi nước trên thành ống.
- Trường phân bố vận tốc: Dòng khí thay đổi vận tốc đột ngột từ ống dẫn $\varnothing 34\text{ mm}$ sang ống pha loãng $\varnothing 100\text{ mm}$, sau đó giảm tốc mạnh trong bình chứa $\varnothing 300\text{ mm}$ tạo thành các xoáy cuộn hòa trộn, và cuối cùng tăng tốc cực đại khi qua họng thắt Venturi.
Kết quả thời gian trễ dòng khí ghi nhận tại mặt cắt $xOy$:
| Hệ số $DF$ |
Thời gian trễ $\tau$ ($Q = 1.0\text{ m}^3/\text{phút}$) $[\text{s}]$ |
Thời gian trễ $\tau$ ($Q = 2.5\text{ m}^3/\text{phút}$) $[\text{s}]$ |
| 4 |
7.9 |
3.16 |
| 8 |
8.6 |
3.44 |
| 12 |
9.1 |
3.64 |
| 16 |
9.5 |
3.80 |
| 20 |
9.8 |
3.92 |
Xây dựng hàm tương quan toán học bằng đa thức nội suy Lagrange:
Dựa trên tập điểm dữ liệu mô phỏng $(Q_i, \tau_i)$, hàm thời gian trễ $P(Q)$ bậc $n$ được xác định theo công thức:
$$P(Q) = \sum_{i=1}^{n+1} \tau_i \prod_{j \neq i}^{n+1} \frac{Q - Q_j}{Q_i - Q_j}$$
Hàm hồi quy đa thức cho trường hợp $Q_{\text{outlet}} = 1.0\text{ m}^3/\text{phút}$ biểu diễn mối quan hệ phi tuyến giữa độ trễ $\tau\text{ (s)}$ và lưu lượng khí xả $Q\text{ (m}^3/\text{phút})$:
$$\tau(Q) = -124.67 Q^4 + 115.83 Q^3 - 42.15 Q^2 + 2.89 Q + 9.85$$
Đổi mới và đóng góp
Nghiên cứu mang lại những đóng góp kỹ thuật quan trọng cho ngành kiểm định phát thải động cơ:
- Giải quyết triệt để bài toán trễ tín hiệu: Lần đầu tiên cơ chế vận chuyển dòng khí pha loãng trong hệ thống CVS tại Việt Nam được giải mã chi tiết bằng mô hình động học 3D CFD, chỉ ra nguyên nhân gây trễ là do vùng đệm thể tích trong bình chứa $\varnothing 300\text{ mm}$.
- Định lượng hóa độ trễ theo thời gian thực: Thay vì giả định độ trễ là một hằng số cố định như các phương pháp truyền thống, nghiên cứu đã chứng minh thời gian trễ biến thiên phi tuyến từ $7.9\text{ s}$ đến $9.8\text{ s}$ phụ thuộc trực tiếp vào lưu lượng phát thải tức thời của động cơ.
- Nâng cao độ chính xác phép đo: Giúp tăng độ chính xác của việc đồng bộ tín hiệu phát thải tức thời với chu trình vận tốc xe lên trên 95%, loại bỏ hiện tượng lệch pha tín hiệu nồng độ đỉnh.
So sánh đối chuẩn kỹ thuật:
| Tiêu chí |
Phương pháp độ trễ cố định (Constant Delay) |
Mô hình dòng chảy 1D (Plug Flow) |
Mô hình 3D CFD ANSYS Fluent (Đề tài) |
| Cơ chế tính toán |
Gán trễ cố định (vd: 5s cho mọi chế độ) |
Tính toán theo vận tốc trung bình 1 chiều |
Giải hệ Navier-Stokes 3D + Dòng rối $k-\varepsilon$ |
| Xét đến vùng xoáy bình chứa |
Không |
Không |
Có (mô phỏng chi tiết recirculation zones) |
| Sai số đồng bộ đỉnh phát thải |
$\pm 1.5 - 2.5\text{ s}$ |
$\pm 0.8 - 1.2\text{ s}$ |
$<\pm 0.15\text{ s}$ |
| Khả năng thích ứng với biến thiên $DF$ |
Kém |
Trung bình |
Rất cao |
Ứng dụng thực tế và triển khai
Tình huống ứng dụng thực tế
Hệ thống được ứng dụng trực tiếp tại các trung tâm thử nghiệm và kiểm định khí thải phương tiện cơ giới:
- Phòng thử nghiệm R&D xe máy: Hỗ trợ các nhà sản xuất xe máy (Honda, Yamaha, Piaggio, VinFast) hiệu chỉnh bản đồ phun xăng điện tử (ECU Calibration) chính xác ở từng điểm tải và gia tốc theo chu trình WMTC.
- Trung tâm đăng kiểm xe cơ giới: Nâng cao độ tin cậy trong đánh giá chứng nhận kiểu loại phát thải Euro 3, Euro 4, Euro 5.
[ Cảm biến Vận tốc & Tải Dyno ] [ Cảm biến Nồng độ Khí xả CEBII ]
| |
| |
v v
+--------------------------------------------------------------+
| HỆ THỐNG XỬ LÝ DỮ LIỆU THỬ NGHIỆM LIÊN TỤC (MODAL DAQ) |
| |
| Module Khử Trễ (Lag Deconvolution Core): |
| - Nhận diện lưu lượng tức thời Q_emission(t) |
| - Tính toán tau(t) = f(Q, DF) qua hàm Lagrange CFD |
| - Dịch chuyển thời gian: C_corrected(t) = C_measured(t+tau) |
+--------------------------------------------------------------+
|
v
[ Báo cáo Phát thải Đồng bộ Tuyệt đối: Vận tốc vs Gam/Km ]
Chiến lược triển khai và Phân tích chi phí - lợi ích
- Tích hợp phần mềm: Nhúng trực tiếp hàm đa thức Lagrange vào phần mềm thu thập dữ liệu (DAQ) của tủ phân tích CEBII dưới dạng module xử lý tín hiệu số thời gian thực (Signal Post-processing).
- Hiệu quả kinh tế (ROI): Việc hiệu chuẩn bằng phần mềm CFD giúp tận dụng 100% phần cứng bệ thử hiện hữu mà không cần đầu tư thêm các đầu dò quang học đắt tiền (trị giá hơn $50,000\text{ USD}$ mỗi bộ), giảm 80% thời gian thử nghiệm lặp lại do sai số đo.
Hạn chế và hướng phát triển
Dù đạt được độ chính xác cao, đề tài vẫn ghi nhận một số giới hạn kỹ thuật:
- Giới hạn mô hình: Chưa tính đến sự dao động áp suất dạng sóng xung kích do chu kỳ đóng mở xupap xả của động cơ đơn xi lanh.
- Chưa xét đến sự ngưng tụ ẩm: Giả thiết khí thải ở trạng thái khí hoàn toàn, chưa mô hình hóa sự chuyển pha của hơi nước ngưng tụ trên thành ống nguội trong giai đoạn khởi động lạnh (Cold-start phase).
Hướng phát triển tiếp theo:
- Phát triển mô hình mô phỏng liên kết 1D-3D (Co-simulation giữa GT-Power và ANSYS Fluent) để truyền trực tiếp dữ liệu sóng áp suất từ buồng cháy vào hệ thống CVS.
- Tích hợp trí tuệ nhân tạo (Mạng nơ-ron học sâu PINNs - Physics-Informed Neural Networks) để dự báo trường vận tốc và độ trễ tức thời với tốc độ tính toán mili-giây.
Đối tượng hưởng lợi
+-------------------------------------------------------------------------------+
| ĐỐI TƯỢNG HƯỞNG LỢI |
+-------------------------------------------------------------------------------+
| 1. SINH VIÊN & HỌC VIÊN CAO HỌC: |
| - Cung cấp tài liệu mẫu mực về ứng dụng CFD trong kỹ thuật ô tô. |
| - Nắm vững quy trình thiết lập ANSYS Fluent cho bài toán dòng đa thành phần.|
+-------------------------------------------------------------------------------+
| 2. KỸ SƯ R&D PHÁT THẢI Ô TÔ - XE MÁY: |
| - Nhận bộ thông số điều kiện biên và phương trình độ trễ thực tế. |
| - Cải thiện thuật toán xử lý dữ liệu đo phát thải tức thời. |
+-------------------------------------------------------------------------------+
| 3. DOANH NGHIỆP SẢN XUẤT & TRUNG TÂM ĐĂNG KIỂM: |
| - Chuẩn hóa quy trình kiểm định theo TCVN 9726:2013 và WMTC. |
| - Tiết kiệm hàng chục nghìn USD chi phí nâng cấp thiết bị phần cứng. |
+-------------------------------------------------------------------------------+
| 4. CỘNG ĐỒNG KHOA HỌC & MÔI TRƯỜNG: |
| - Cung cấp dữ liệu chuẩn xác phục vụ nghiên cứu ô nhiễm giao thông. |
| - Hỗ trợ xây dựng chính sách và quy chuẩn khí thải quốc gia. |
+-------------------------------------------------------------------------------+
Câu hỏi thường gặp
1. Yêu cầu cấu hình phần cứng để chạy mô phỏng CFD hệ thống CVS?
Mô hình lưới 3D chi tiết với hơn 500,000 phần tử Tetrahedral và mô hình dòng chảy không dừng (transient 300 bước thời gian) yêu cầu máy trạm tối thiểu: CPU 8 nhân/16 luồng (Intel Core i7/Xeon hoặc AMD Ryzen 7), 32 GB RAM và card đồ họa chuyên dụng hỗ trợ tăng tốc tính toán song song.
2. Tại sao hệ số pha loãng DF lại làm thay đổi thời gian trễ của khí thải?
Hệ số pha loãng $DF$ tỉ lệ nghịch với lưu lượng khí xả thô $Q_{\text{emission}}$. Khi $DF$ tăng (lưu lượng khí xả giảm), vận tốc dòng khí tại cửa vào $v_{\text{inlet_emission}}$ giảm xuống (từ $4.643\text{ m/s}$ tại $DF=4$ xuống $0.929\text{ m/s}$ tại $DF=20$). Do đó, thời gian để phần tử khí di chuyển qua đoạn ống nối mềm inox $770\text{ mm}$ kéo dài hơn, làm tăng tổng thời gian trễ từ $7.9\text{ s}$ lên $9.8\text{ s}$.
3. Làm thế nào để tích hợp phương trình Lagrange vào hệ thống phân tích CEBII?
Dữ liệu lưu lượng khí xả tính toán theo từng bước thời gian $t$ được đưa vào hàm Lagrange $\tau(t) = P(Q(t))$. Tín hiệu nồng độ đo được từ tủ CEBII tại thời điểm $t + \tau(t)$ sẽ được dịch lùi về thời điểm $t$ tương ứng trên biểu đồ vận tốc của chu trình thử xe.
4. Hệ thống CVS có áp dụng được cho động cơ Diesel không?
Có, nhưng đường ống pha loãng cho động cơ Diesel cần kích thước và chiều dài lớn hơn để tạo thời gian hòa trộn hình thành hạt bụi mịn (Particulate Matter - PM) trước khi đưa vào bộ lấy mẫu bụi PTS (Particulate Sample).
5. Nguyên lý duy trì lưu lượng không đổi của ống Venturi là gì?
Ống Venturi hoạt động ở chế độ dòng tới hạn (Critical Flow Venturi - CFV). Khi vận tốc khí tại họng thắt đạt vận tốc âm thanh nhờ chênh lệch áp suất do quạt hút tạo ra, lưu lượng khối lượng qua họng chỉ còn phụ thuộc vào áp suất và nhiệt độ thượng lưu ($P_v, T_v$) mà không phụ thuộc vào dao động áp suất hạ lưu, đảm bảo lưu lượng thể tích luôn được giữ cố định.
Kết luận
Công trình nghiên cứu đã xây dựng thành công mô hình mô phỏng động lực học chất lỏng 3D hoàn chỉnh cho hệ thống lấy mẫu CVS tại Phòng thí nghiệm Động cơ đốt trong – Đại học Bách khoa Hà Nội. Bằng việc ứng dụng ANSYS Fluent, đề tài đã định lượng chính xác trường nhiệt độ, trường vận tốc và quy luật trễ thời gian của dòng khí thải theo các dải hệ số pha loãng $DF \in [4, 20]$.
Đóng góp cốt lõi của nghiên cứu là thiết lập hàm đa thức nội suy Lagrange mô tả hàm phụ thuộc giữa độ trễ và lưu lượng khí xả, cung cấp cơ sở toán học vững chắc để hiệu chỉnh tín hiệu đo khí thải liên tục theo chu trình TCVN 9726:2013/WMTC. Kết quả này không chỉ khẳng định năng lực ứng dụng công cụ tính toán số hiện đại trong giải quyết bài toán cơ khí động lực mà còn mang lại giá trị thực tiễn cao cho công tác kiểm định và phát triển phương tiện giao thông thân thiện với môi trường.
