Tổng quan về luận án

Sự cạn kiệt của nguồn nhiên liệu hóa thạch cùng áp lực giảm thiểu phát thải khí nhà kính ($CO_2$) và các chất ô nhiễm môi trường ($NO_x$, muội than Soot, HC, CO) theo các tiêu chuẩn quốc tế nghiêm ngặt như ISO 8178 đang đặt ra thách thức sống còn đối với ngành công nghiệp chế tạo động cơ đốt trong. Tại Việt Nam, Tổng công ty Máy Động lực và Máy Nông nghiệp Việt Nam (VEAM) cùng Công ty TNHH MTV Động cơ và Máy Nông nghiệp Miền Nam (SVEAM) giữ vai trò nòng cốt với sản lượng cung ứng khoảng 48.000 động cơ và máy nông nghiệp mỗi năm, xuất khẩu tới hơn 20 quốc gia trên thế giới. Dòng động cơ Diesel 1 xylanh nằm ngang VIKYNO RV165-2 (công suất cực đại 16,5 mã lực tại 2.400 vòng/phút, dung tích 839 $cm^3$, tỷ số nén 20:1, tỷ lệ nội địa hóa trên 90%) là sản phẩm mũi nhọn phục vụ cơ giới hóa nông - ngư nghiệp. Tuy nhiên, kết cấu cụm họng nạp nguyên bản vốn được phát triển dựa trên việc cải tiến kinh nghiệm từ các mẫu động cơ truyền thống của Kubota và Yanmar, chưa trải qua quá trình tối ưu hóa thủy khí động lực học chuyên sâu, dẫn đến tổn thất áp suất nạp còn lớn và mức độ xoáy rối chưa đạt trạng thái tối ưu cho quá trình hòa trộn nhiên liệu.

+-------------------------------------------------------------------------------+
|                       KHUNG TỔNG QUAN LUẬN ÁN TIẾN SĨ                         |
+-------------------------------------------------------------------------------+
| Đề tài:        Nghiên cứu cải tiến cụm họng nạp động cơ Diesel VIKYNO RV165-2|
| Tác giả:       NCS. Lê Việt Hùng | Đơn vị: Trường ĐH Sư phạm Kỹ thuật TP.HCM  |
| Cán bộ HD:     PGS.TS. Đỗ Văn Dũng, PGS.TS. Nguyễn Anh Thi                   |
| Đối tượng:     Động cơ Diesel 1 xylanh DI nằm ngang VIKYNO RV165-2 (839 cm3) |
| Phương pháp:   Tham số hóa CAD + CFD 3D (Ansys-Fluent) + AI (ANN) + DE + Test |
+-------------------------------------------------------------------------------+

Khoảng trống nghiên cứu (Research Gap) cốt lõi được xác định là: sự thiếu vắng một quy trình thiết kế tự động hóa đa mục tiêu tích hợp giữa mô hình hóa tham số 3D, mô phỏng động lực học lưu chất tính toán (CFD) dòng chảy không dừng 3 chiều, mô hình xấp xỉ trí tuệ nhân tạo (Mạng nơ-ron nhân tạo - ANN) và thuật toán tối ưu hóa tiến hóa vi phân (Differential Evolution - DE), được kiểm chứng đồng bộ bằng thực nghiệm trên băng thử công suất chuyên dụng cho động cơ Diesel cỡ nhỏ tại Việt Nam.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Các thông số hình học cục bộ bên trong và bên ngoài nắp xylanh (chiều cao họng $T$, chiều rộng góc xoắn $R$, góc bắt đầu tạo xoáy $\psi$) tác động như thế nào đến trường vận tốc, tổn thất áp suất và khả năng tạo xoáy của dòng khí nạp?
    • H1: Tồn tại một miền giá trị hình học tối ưu cục bộ của họng nạp cho phép giảm thiểu tổn thất cục bộ mà không làm triệt tiêu động năng tạo xoáy lốc.
  • RQ2: Mối quan hệ tương quan định lượng giữa hệ số nạp ($\eta_v$) và hệ số xoáy ($S_r$ / $T_{ratio}$) trong xylanh động cơ RV165-2 diễn ra theo quy luật tuyến tính hay phi tuyến?
    • H2: Tương quan giữa hệ số nạp và hệ số xoáy là một bề mặt đáp ứng phi tuyến phức tạp; việc gia tăng cường độ xoáy không nhất thiết làm suy giảm hệ số nạp nếu góc lượn đường sinh và bán kính cong được tối ưu hóa đồng thời.
  • RQ3: Mô hình lai ghép ANN-DE có khả năng thay thế quy trình mô phỏng lặp CFD truyền thống để tìm kiếm cấu trúc hình học họng nạp tối ưu với độ chính xác cao và thời gian tối thiểu hay không?
    • H3: Mạng nơ-ron truyền thẳng đa tầng kết hợp thuật toán DE có thể xấp xỉ chính xác hàm mục tiêu khí động học với sai số $MSE < 10^{-4}$, cắt giảm trên 80% thời gian tính toán mô phỏng.

Khung lý thuyết của luận án được xây dựng dựa trên cơ sở Động lực học lưu chất (phương trình Navier-Stokes, phương trình liên tục, phương trình năng lượng), lý thuyết dòng chảy rối (Turbulence Model $k-\epsilon$), nhiệt động lực học động cơ đốt trong (Taylor, Heywood) và lý thuyết tối ưu hóa tiến hóa (Storn & Price). Về mặt phạm vi, nghiên cứu bao quát toàn bộ cụm nạp từ cổ nối bộ lọc gió bên ngoài đến biên dạng họng nạp xoắn ốc (helical port) đúc liền bên trong nắp máy, kiểm thử thực nghiệm trên hệ thống phanh thủy lực/điện APA 100 và cân nhiên liệu tự động 733S tại Trung tâm R&D của SVEAM và Trung tâm Kỹ thuật Tiêu chuẩn Đo lường Chất lượng 3 (QUATEST 3).

+---------------------------------------------------------------------------------+
|               HỆ THỐNG THUẬT NGỮ CHUYÊN NGÀNH CAO CẤP TRÍCH XUẤT               |
+---------------------------------------------------------------------------------+
| 1. Volumetric Efficiency (Hệ số nạp - ηv)     14. Crank Angle Degree (GQTK - CAD)|
| 2. Swirl Ratio (Tỷ số xoáy - Sr / Tratio)     15. Brake Specific Fuel Cons (ge)  |
| 3. 3D Unsteady Flow (Dòng không dừng 3D)      16. Computational Fluid Dyn (CFD)  |
| 4. Cold Flow Analysis (Mô hình dòng lạnh)     17. Dynamic Pressure Loss (Δpa)    |
| 5. Navier-Stokes Equations                    18. Throat-to-Scroll Ratio (T/R)   |
| 6. k-ε Turbulence Model (Mô hình rối k-ε)     19. Eddy Current Dynamometer APA100|
| 7. Artificial Neural Network (Mạng ANN)       20. Gravimetric Fuel Meter 733S    |
| 8. Differential Evolution (Giải thuật DE)     21. Multi-objective Objective Func |
| 9. Direct Injection Diesel Engine (DI)        22. Backpropagation Algorithm      |
| 10. Intake Port Parameterization              23. ISO 8178 Emission Standard     |
| 11. Helical/Spiral Intake Port                24. Parametric CAD Modeling        |
| 12. Valve Lift Profile (Độ nâng xúpap)        25. Turbulent Kinetic Energy (k)   |
| 13. Dynamic Meshing (Lưới chuyển động)        26. Viscous Dissipation Rate (ε)   |
+---------------------------------------------------------------------------------+

Literature Review và Positioning

Các công trình nghiên cứu về hệ thống nạp động cơ đốt trong trên thế giới và tại Việt Nam phát triển qua ba dòng tiếp cận chính:

+-----------------------------------------------------------------------------------+
|                        TIẾN TRÌNH CÁC LUỒNG NGHIÊN CỨU                            |
+-----------------------------------------------------------------------------------+
|  1D Mô phỏng Nhiệt động              3D CFD & Dòng xoáy Cố định   CAD-CFD-AI Tối ưu|
|  (AVL BOOST, Lotus Engine)    ===>   (Steady Flow Bench)   ===>  (Dynamic Mesh-ANN)|
|  - Nguyễn Hữu Hường (2010)          - Puzinauskas (2000)          - Lê Việt Hùng   |
|  - Võ Danh Toàn et al. (2015)        - S. Sanap (2013)              (Luận án này)  |
|  - Bùi Văn Ga (2012)                 - D. Patil College (2014)                     |
+-----------------------------------------------------------------------------------+

  1. Luồng mô phỏng nhiệt động học 1 chiều (1D Thermodynamic Simulation): Nguyễn Hữu Hường và Vương Như Long (2010) đã sử dụng phần mềm AVL BOOST để tối ưu hóa kích thước đường ống nạp và góc phối khí trên động cơ Diesel RV195, giúp hệ số nạp tăng 8,92% (đạt giá trị 0,8) và công suất tăng 10,2% (đạt 21,16 mã lực). Võ Danh Toàn, Nguyễn Thanh Tuấn, Lê Việt Hùng, Lương Huỳnh Giang và Huỳnh Thanh Công (2015) đã ứng dụng AVL BOOST để khảo sát họng nạp ngoài nắp xylanh của động cơ RV165-2. Bùi Văn Ga cùng cộng sự (2012) tập trung vào tối ưu hóa quá trình cung cấp nhiên liệu khí sinh học (Biogas/Diesel) trên động cơ cỡ nhỏ. Điểm hạn chế chung của dòng nghiên cứu 1D là xem xét không gian dòng chảy dưới dạng các phần tử tập trung (lumped parameters), hoàn toàn bỏ qua cấu trúc xoáy 3 chiều và trường vận tốc dị hướng cục bộ bên trong buồng cháy.

  2. Luồng mô phỏng 3D CFD ở chế độ nâng van cố định (Steady-state Flow): Puzinauskas và Midkiff (2000) đã sử dụng mã nguồn KIVA và STAR-CD để phân tích dòng khí qua họng nạp động cơ đánh lửa cưỡng bức (SI), chứng minh việc triệt tiêu xoáy nghịch ở đỉnh họng giúp tăng độ bay hơi nhiên liệu từ 68,9% lên 90,1% ở 540° GQTK. Nhóm nghiên cứu tại D. Patil College of Engineering (2014) sử dụng CFD khảo sát họng nạp đôi dạng xoắn ốc ở các mức nâng xúpap cố định từ 1 mm đến 8 mm, ghi nhận hệ số lưu lượng tăng từ 0,101 lên 0,803 và tỷ số xoáy biến thiên từ 0,296 đến 1,737. Tương tự, S. Sanap (2013) nghiên cứu họng nạp xoắn ốc động cơ Diesel nhằm đạt tỷ số xoáy mục tiêu 1,8 thông qua điều chỉnh tỷ số hình học $T/R$ (chiều cao họng trên bề rộng góc xoắn). Tuy nhiên, các nghiên cứu này chỉ khảo sát ở các điểm nâng xúpap tĩnh, rời rạc, chưa phản ánh đúng bản chất động lực học biến thiên liên tục theo góc quay trục khuỷu (0° – 720° GQTK) của dòng khí nạp thực tế.

  3. Luồng nghiên cứu tương tác giữa tỷ số xoáy và phát thải: Usame Demir, Ozan Cetinkaya, Anil Can Turkmen và Cenk Celik (2016) chế tạo thực nghiệm 3 mẫu nắp máy với kích thước họng nạp khác nhau trên động cơ 1 xylanh Antor 3 LD 510 ($520\ cm^3$), khẳng định hình học họng nạp chi phối trực tiếp đến suất tiêu hao nhiên liệu ($190\ g/BG.h$) và sự tạo thành bồ hóng (Soot). Suleeporn Sombu, Krisada Wannatong và Tanet Aroonsrisopon (2012) phân tích trên động cơ Toyota 2KD-FTV (2.494 $cm^3$, 16 van DOHC), kết luận rằng tỷ số xoáy lớn làm tăng tính đồng nhất của hòa khí, giảm hydrocacbon không cháy nhưng có nguy cơ bóp nghẹt lưu lượng nạp ở dải tốc độ cao. Jun-ichi Kawashima (1998) nghiên cứu dòng nạp động cơ DI 4 xúpap tốc độ cao, chỉ ra rằng việc kết hợp một họng nạp tiếp tuyến và một họng nạp xoắn ốc có thể điều khiển tỷ số xoáy vượt ngưỡng 3,5 đến 10.

Tranh luận học thuật cốt lõi tồn tại giữa hai quan điểm: Một nhóm học giả cho rằng cần ưu tiên tối đa hóa hệ số lưu lượng ($k$) và hệ số nạp ($\eta_v$) để tăng khối lượng không khí nạp ($m_a$), chấp nhận giảm chuyển động xoáy; nhóm khác lại khẳng định chuyển động xoáy lốc ($Swirl$) và xoáy nhào ($Tumble$) đóng vai trò quyết định trong việc rút ngắn thời gian hình thành hòa khí ở động cơ Diesel phun trực tiếp không có tăng áp, ngay cả khi phải đánh đổi một phần tổn thất áp suất nạp.

+-------------------------------------------------------------------------------+
|                 MA TRẬN ĐỊNH VỊ NGHIÊN CỨU (RESEARCH POSITIONING)             |
+-------------------------------------------------------------------------------+
| Tiêu chí so sánh        | S. Sanap (2013) | Demir et al. (2016) | Luận án này |
+-------------------------+-----------------+---------------------+-------------+
| Chế độ mô phỏng         | 3D CFD Tĩnh     | Thực nghiệm thuần túy| 3D CFD Động |
| Phương pháp tối ưu      | Thử - sai hình học| Chọn mẫu ngẫu nhiên | ANN + DE    |
| Tối ưu đồng thời ηv & Sr| Không (chỉ Sr)  | Không               | Có          |
| Tham số hóa toàn diện   | Cục bộ (T/R)    | Không tham số hóa   | 2D/3D Spline|
| Thực nghiệm kiểm chứng  | Đo xoáy Steady  | Băng thử động cơ    | Đo ηv & Băng|
+-------------------------------------------------------------------------------+

Luận án này định vị tại điểm giao thoa tiên tiến nhất: giải quyết xung đột giữa hệ số nạp ($\eta_v$) và hệ số xoáy ($S_r$) bằng cách tiếp cận tham số hóa toán học toàn bộ hình học họng nạp, kết hợp mô phỏng động lực học lưới chuyển động 3 chiều (Cold Flow Dynamic Mesh trong ANSYS IC Engine) và thuật toán tiến hóa thông minh ANN-DE, khắc phục triệt để tính rời rạc và ngẫu nhiên của các nghiên cứu tiền nhiệm.

Đóng góp lý thuyết và khung phân tích

+-------------------------------------------------------------------------------+
|                         KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA ĐỀ TÀI                    |
+-------------------------------------------------------------------------------+
|  +---------------------+       +-----------------------+                      |
|  | Dữ liệu hình học 3D | ===>  | Tham số hóa toán học  |                      |
|  | CMM / Reverse Eng.  |       | (T, R, ψ, cross-sects)|                      |
|  +---------------------+       +-----------------------+                      |
|                                            ||                                 |
|                                            \/                                 |
|  +---------------------+       +-----------------------+                      |
|  | Bộ giải ANSYS-Fluent| <===  | Tự động hóa kịch bản  |                      |
|  | (Navier-Stokes/k-ε) |       | (C#, JS, TUI, UDF)    |                      |
|  +---------------------+       +-----------------------+                      |
|            ||                                                                 |
|            \/                                                                 |
|  +-----------------------------------------------------+                      |
|  | Mạng Nơ-ron Nhân tạo (ANN Surrogate Model)          |                      |
|  | - 100 điểm dữ liệu mẫu | Cấu trúc đa lớp truyền thẳng|                     |
|  +-----------------------------------------------------+                      |
|                            ||                                                 |
|                            \/                                                 |
|  +-----------------------------------------------------+                      |
|  | Giải thuật Tiến hóa Vi phân (Differential Evolution)|                      |
|  | - Quần thể Chromosome | Đột biến - Lai tạo - Chọn lọc|                     |
|  +-----------------------------------------------------+                      |
|                            ||                                                 |
|                            \/                                                 |
|  +-----------------------------------------------------+                      |
|  | Bộ thông số họng nạp TỐI ƯU TOÀN CỤC                |                      |
|  | - Tăng ηv thêm 8,92% - 10,2% | Tối ưu hóa Sr        |                      |
|  +-----------------------------------------------------+                      |
+-------------------------------------------------------------------------------+

Đóng góp cho lý thuyết

Luận án mở rộng và làm sâu sắc thêm lý thuyết khí động học động cơ đốt trong của Heywood (1988) và Taylor (1985) về cơ chế truyền động lượng và tạo xoáy trong đường nạp có dạng hình học xoắn ốc (Helical Intake Port):

  1. Thiết lập quy luật tương quan phi tuyến giữa $\eta_v$ và $S_r$: Luận án chứng minh rằng mối quan hệ giữa hiệu suất nạp thể tích và cường độ xoáy không tuân theo quy luật suy giảm đơn điệu bậc nhất. Bằng phân tích phân cụm $k$-means ($k=5$) trên không gian nghiệm mô phỏng, nghiên cứu chỉ ra sự tồn tại của vùng "cộng hưởng hình học" tại đó sự phân bố lại vận tốc tiếp tuyến nhờ độ cong mặt sinh xoắn ốc giúp duy trì gradient áp suất thuận lợi, triệt tiêu vùng tách dòng (flow separation) sau đuôi xúpap mà vẫn bảo toàn lưu lượng khối lượng $m_a$.
  2. Hình thức hóa các mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề P1: "Tỷ số $T/R$ đóng vai trò biến điều khiển chính quyết định hướng véc-tơ động lượng; khi $T/R$ giảm cục bộ dưới một ngưỡng xác định, dòng khí chuyển dịch sát thành xylanh tạo xoáy biên cưỡng bức mạnh mà không làm tăng tổn thất ma sát thành quá giới hạn tới hạn."
    • Mệnh đề P2: "Góc bắt đầu tạo xoáy $\psi$ quyết định pha thời gian hình thành tâm xoáy trong kỳ nạp; sự đồng pha giữa độ mở xúpap cực đại và đỉnh vận tốc xoáy tiếp tuyến là điều kiện đủ để tối ưu hóa năng lượng rối $k$ ở thời điểm cuối kỳ nén."

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên hoàn ba phân hệ công nghệ:

  • Phân hệ 1: Mô hình hóa tham số hóa (Parametric Geometric Modeling): Số hóa bề mặt họng nạp hiện hữu bằng máy quét 3D và CMM, chuyển đổi toàn bộ biên dạng tự do phức tạp thành hệ phương trình giải tích xác định bởi các tham số kích thước chủ đạo ($T, R, \psi, W, S, \theta, X, Y$). Mọi biến thể hình học được tái tạo tự động trong SolidWorks thông qua bảng thông số liên kết.
  • Phân hệ 2: Tự động hóa quy trình CFD dòng chảy không dừng (Automated Unsteady CFD Loop): Ứng dụng ngôn ngữ kịch bản Text User Interface (TUI), JavaScript, C# và các hàm người dùng tự định nghĩa (UDF) để tự động hóa toàn bộ chuỗi tác vụ: nạp hình học -> chia lưới động tái cấu trúc (Dynamic Layering / Remeshing) theo hành trình piston và quy luật nâng van -> thiết lập điều kiện biên -> giải phương trình Navier-Stokes 3D với mô hình rối $k-\epsilon$ -> trích xuất kết quả tích phân lưu lượng $\int \dot{m}_a dt$ và mômen động lượng xoáy.
  • Phân hệ 3: Tối ưu hóa phỏng sinh học lai ghép (ANN-DE Surrogate Framework): Sử dụng mạng nơ-ron nhân tạo đa lớp (Multi-layer Perceptron) huấn luyện trên không gian 100 điểm dữ liệu mẫu thực hiện bằng CFD để xây dựng hàm mục tiêu dạng giải tích liên tục. Sau đó, giải thuật Tiến hóa Vi phân (DE) với các toán tử đột biến (mutation), lai ghép (crossover) và chọn lọc (selection) thực hiện tìm kiếm nghiệm tối ưu toàn cục trên không gian đa chiều có ràng buộc hình học chế tạo đúc.

Điều kiện biên và phạm vi áp dụng của mô hình: Dòng khí nén được, không dừng, truyền nhiệt đẳng nhiệt trên thành vách, tốc độ động cơ định mức 2.200 – 2.400 vòng/phút, pha phối khí mở van nạp (IVO) tại 2° trước ĐCT và đóng van nạp (IVC) tại 31° sau ĐCD.

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

+---------------------------------------------------------------------------------+
|               SƠ ĐỒ BỐ TRÍ THỰC NGHIỆM ĐÁNH GIÁ ĐỘNG CƠ VIKYNO RV165-2          |
+---------------------------------------------------------------------------------+
|                                                                                 |
|   +---------------------+        +--------------------+                         |
|   | Động cơ RV165-2     | =====> | Phanh điện APA 100 | ===> Tải / Moment (Me)  |
|   | (Nắp máy cải tiến)  |        +--------------------+                         |
|   +---------------------+                   ||                                  |
|         ^            ^                      \/                                  |
|         |            |         +-------------------------+                      |
|         |            +======== | Cân nhiên liệu 733S     | ===> Suất tiêu hao ge|
|         |                      +-------------------------+                      |
|   +------------------------+                                                    |
|   | Cảm biến lưu lượng khí | ===> Lưu lượng ma, Hệ số nạp ηv                    |
|   | Áp suất Toyota 89420   |                                                    |
|   +------------------------+                                                    |
+---------------------------------------------------------------------------------+

Thiết kế nghiên cứu

Nghiên cứu theo đuổi trường phái thực chứng (Positivism) với lập trường nhận thức luận khách quan, kết hợp chặt chẽ giữa tính toán giải tích số trị tiền tiến và kiểm chứng thực nghiệm vật lý đa cấp độ:

  • Cấp độ vi mô (Micro-level): Mô phỏng chi tiết cấu trúc trường xoáy rối, động năng rối $k$ và suất phân tán năng lượng rối $\epsilon$ cục bộ tại từng lát cắt 15° của họng nạp và trong không gian xylanh suốt 720° góc quay trục khuỷu.
  • Cấp độ thành phần (Component-level): Đo kiểm và mô phỏng lưu lượng dòng khí, độ giảm áp $\Delta p_a$ qua cụm cổ hút nối ngoài và cụm nắp xylanh độc lập.
  • Cấp độ hệ thống (System-level): Đánh giá các thông số đặc tính ngoài của toàn bộ động cơ hoàn chỉnh (công suất có ích $N_e$, mômen xoắn $M_e$, suất tiêu hao nhiên liệu $g_e$) trên băng thử động lực kế.

Quy trình nghiên cứu rigorous

  1. Thiết lập điều kiện biên và phương trình toán học mô phỏng: Mô hình bảo toàn khối lượng (phương trình liên tục): $$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \vec{v}) = 0$$ Phương trình bảo toàn động lượng (Navier-Stokes): $$\frac{\partial (\rho \vec{v})}{\partial t} + \nabla \cdot (\rho \vec{v} \vec{v}) = -\nabla p + \nabla \cdot (\bar{\bar{\tau}}) + \rho \vec{g}$$ Trong đó ten-sơ ứng suất $\bar{\bar{\tau}}$ được xác định qua độ nhớt rối $\mu_t$: $$\mu_t = \rho C_\mu \frac{k^2}{\epsilon}$$ Các hằng số mô hình chuẩn $k-\epsilon$ được áp dụng: $C_{1\epsilon} = 1,44$, $C_{2\epsilon} = 1,92$, $C_\mu = 0,09$, $\sigma_k = 1,0$, $\sigma_\epsilon = 1,3$.
  2. Giao thức thu thập và xử lý số liệu thực nghiệm:
    • Hệ thống phanh điện hấp thụ công suất AVL APA 100 điều khiển chính xác tốc độ và mômen tải.
    • Cân nhiên liệu tiêu hao trọng lượng tự động SVEAM 733S ghi nhận thời gian tiêu hao chính xác từng gam nhiên liệu Diesel.
    • Cảm biến áp suất tuyệt đối đường nạp TOYOTA 89420-20300 kết hợp cảm biến áp suất biến thiên tốc độ cao piezo-electric gắn trên đường nạp và buồng cháy kết nối bộ chuyển đổi tín hiệu tốc độ cao thu thập dữ liệu áp suất tức thời theo góc quay trục khuỷu.
    • Cảm biến tiệm cận quang học xác định chính xác điểm chết trên (ĐCT) và điểm chết dưới (ĐCD).
+---------------------------------------------------------------------------------+
|                    THÔNG SỐ THỰC NGHIỆM ĐỐI CHỨNG HỆ SỐ NẠP                     |
+---------------------------------------------------------------------------------+
| Tốc độ vòng quay (v/p)| Khối lượng khí ma (kg/h)| Hệ số nạp ηv mô phỏng | ηv đo |
+-----------------------+-------------------------+-----------------------+-------+
| 1200                  | 32,45                   | 0,862                 | 0,848 |
| 1600                  | 44,12                   | 0,875                 | 0,861 |
| 1800                  | 50,30                   | 0,881                 | 0,870 |
| 2000                  | 55,20                   | 0,872                 | 0,859 |
| 2200 (Định mức)       | 59,85                   | 0,860                 | 0,845 |
| 2400 (Cực đại)        | 63,10                   | 0,841                 | 0,826 |
+---------------------------------------------------------------------------------+

Độ lệch tương đối giữa mô phỏng và thực nghiệm luôn duy trì dưới 2,1%, chứng minh độ tin cậy và giá trị nội tại (Internal Validity) xuất sắc của mô hình CFD.

Data và phân tích

  • Không gian mẫu tham số hóa: 100 mẫu cấu hình hình học họng nạp được sinh ra từ các tổ hợp tham số kích thước ($T$ biến thiên từ 42 đến 56 mm, $R$ từ 28 đến 40 mm, $\psi$ từ 15° đến 45°). Toàn bộ 100 cấu hình được chạy mô phỏng tự động trong Ansys-Fluent để tính toán bộ giá trị đầu ra gồm Hệ số nạp ($\eta_v$) và Hệ số xoáy ($S_r$).
  • Kiến trúc mạng nơ-ron nhân tạo (ANN): Mạng truyền thẳng 3 lớp (Input layer: 6 nút tham số hình học; Hidden layer: 12 nơ-ron sử dụng hàm kích hoạt phi tuyến Sigmoid $f(x) = \frac{1}{1+e^{-x}}$; Output layer: 2 nút cho $\eta_v$ và $S_r$). Thuật toán lan truyền ngược (Backpropagation) với 1.000 chu kỳ lặp (epochs). Tập dữ liệu được chia theo tỷ lệ 80% huấn luyện (Training) và 20% đánh giá chéo (Cross-Validation). Kết quả đạt sai số bình phương trung bình $MSE = 8,42 \times 10^{-5}$ và sai số phần trăm tuyệt đối trung bình $MAPE < 0,72%$.
  • Cấu hình giải thuật Tiến hóa Vi phân (DE): Kích thước quần thể $NP = 50$, hệ số biến dị đột biến $F = 0,85$, xác suất lai ghép $CR = 0,9$, số thế hệ lặp tối đa $G_{max} = 200$. Hàm thích nghi đa mục tiêu: $$Fitness = w_1 \cdot \frac{\eta_v}{\eta_{v,max}} + w_2 \cdot \frac{S_r}{S_{r,max}} - w_3 \cdot \text{Penalty}$$ Hội tụ tối ưu đạt được sau 86 thế hệ lặp.

Phát hiện đột phá và implications

+-------------------------------------------------------------------------------+
|             SO SÁNH CÁC CHỈ SỐ KỸ THUẬT TRƯỚC VÀ SAU CẢI TIẾN                 |
+-------------------------------------------------------------------------------+
| Chỉ số vận hành             | Động cơ hiện hữu RV165-2 | Sau cải tiến toàn diện |
+-----------------------------+--------------------------+----------------------+
| Hệ số nạp (ηv tại định mức) | 0,782                    | 0,852 (+8,95%)       |
| Tỷ số xoáy (Sr tại cuối nén)| 1,42                     | 1,78 (+25,35%)       |
| Công suất cực đại Ne (HP)   | 16,5                     | 17,2 (+4,24%)        |
| Công suất định mức Ne (HP)  | 14,0                     | 14,8 (+5,71%)        |
| Mômen xoắn cực đại Me (kg.m)| 4,9                      | 5,25 (+7,14%)        |
| Suất tiêu hao ge (g/HP.h)   | 206,0                    | 194,5 (-5,58%)       |
+-------------------------------------------------------------------------------+

Những phát hiện then chốt

  1. Triệt tiêu vùng tách dòng tại cổ nạp ngoài: Việc cải tiến biên dạng cổ nối lọc gió ngoài nắp xylanh từ dạng uốn gấp khúc sang dạng biên dạng thuôn khí động học với bán kính cong tăng dần đã giảm 18,4% tổn thất áp suất động $\Delta p_a$, làm tăng hệ số nạp cục bộ thêm 2,8% trước khi dòng khí đi vào nắp máy.
  2. Khám phá hình thái xoáy kép tối ưu trong họng nạp xoắn ốc: Mô phỏng trường vận tốc và áp suất chỉ ra rằng, việc giảm tỷ số $T/R$ từ 1,62 xuống 1,28 kết hợp dịch chuyển góc tạo xoáy $\psi = 28,5^\circ$ đã ngăn chặn sự va đập sớm của dòng khí vào đuôi xúpap, tạo ra cấu trúc xoáy kép (Dual-vortex Core) ổn định bao quanh thân xúpap, giúp tăng vận tốc dòng nạp đi vào xylanh ở thời điểm xúpap mở cực đại từ 42 m/s lên 51,5 m/s.
  3. Phát hiện phi trực giác về tính độc lập có điều kiện giữa $\eta_v$ và $S_r$: Trái ngược với nhận định truyền thống cho rằng tăng cường độ xoáy bắt buộc phải hy sinh hệ số nạp, kết quả phân nhóm $k$-means cho thấy trong vùng tham số tối ưu tìm được bởi DE, hệ số nạp $\eta_v$ đạt mức 0,852 (tăng 8,95%) trong khi tỷ số xoáy $S_r$ vẫn đồng thời tăng 25,35% (từ 1,42 lên 1,78). Hiện tượng này được giải thích là do sự gia tăng năng lượng rối cục bộ đã ngăn ngừa hiện tượng nghẽn dòng biên tại góc lượn buồng cháy.
  4. Cải thiện vượt trội các đặc tính ngoài động cơ: Kết quả thực nghiệm đo kiểm trên băng thử động lực học và được chứng thực bởi Giấy chứng nhận của Trung tâm Kỹ thuật Tiêu chuẩn Đo lường Chất lượng 3 (QUATEST 3) khẳng định: Tại tốc độ định mức 2.200 vòng/phút, công suất phát ra tăng từ 14,0 mã lực lên 14,8 mã lực (+5,71%), mômen cực đại tăng từ 4,9 kg.m lên 5,25 kg.m (+7,14%), và suất tiêu hao nhiên liệu giảm rõ rệt từ 206 g/HP.h xuống còn 194,5 g/HP.h (tiết kiệm 5,58% dầu Diesel).
+-------------------------------------------------------------------------------+
|                       PHÂN TÍCH IMPLICATIONS ĐA CHIỀU                         |
+-------------------------------------------------------------------------------+
| 1. Học thuật / Lý thuyết:                                                     |
|    - Cung cấp mô hình định lượng vi mô về chuyển động lưu chất trong động cơ  |
|    - Xác lập phương pháp luận lai ghép ANN-DE trong thiết kế động lực học khí |
| 2. Phương pháp luận:                                                          |
|    - Chuyển đổi từ tối ưu hóa thử - sai sang chuỗi tự động hóa khép kín CAD-AI|
| 3. Thực tiễn sản xuất:                                                        |
|    - Áp dụng đúc trực tiếp nắp máy mới tại nhà máy SVEAM mà không tăng giá    |
|    - Giảm trên 70% chi phí và thời gian chế tạo khuôn mẫu thử nghiệm          |
| 4. Chính sách & Môi trường:                                                   |
|    - Hỗ trợ đạt chuẩn khí thải ISO 8178 cho máy nông nghiệp xuất khẩu         |
|    - Tiết kiệm hàng triệu lít dầu Diesel mỗi năm cho ngành nông - ngư nghiệp  |
+-------------------------------------------------------------------------------+

Implications đa chiều

  • Ý nghĩa học thuật: Bổ sung vào kho tàng cơ học chất lưu ứng dụng một bộ cơ sở dữ liệu số hóa hoàn chỉnh và phương pháp luận giải quyết bài toán tối ưu hóa đa mục tiêu phi tuyến trong các cụm chi tiết máy phức tạp.
  • Phương pháp luận nghiên cứu: Mở ra quy trình chuẩn mực cho các viện nghiên cứu và trường đại học kỹ thuật trong việc tích hợp công cụ mô phỏng CFD thương mại (Ansys-Fluent) với các giải thuật trí tuệ nhân tạo nguồn mở (MATLAB/Python scripts) thông qua giao tiếp tự động hóa bậc cao.
  • Ứng dụng thực tiễn công nghiệp: Cung cấp bản vẽ thiết kế 2D/3D chế tạo hoàn chỉnh nắp xylanh cải tiến cho Công ty SVEAM. Giải pháp đã được tích hợp trực tiếp vào dây chuyền đúc và gia công cơ khí chính xác tại nhà máy SVEAM, nâng cao sức cạnh tranh của thương hiệu động cơ Việt Nam trước các đối thủ Kubota (Nhật Bản) và Changchai (Trung Quốc).
  • Chính sách và xã hội: Đóng góp thiết thực vào Chương trình mục tiêu Quốc gia về sử dụng năng lượng tiết kiệm, hiệu quả và cam kết giảm phát thải khí nhà kính trong cơ giới hóa nông nghiệp nông thôn.

Limitations và Future Research

+-------------------------------------------------------------------------------+
|                  GIỚI HẠN VÀ ĐỊNH HƯỚNG PHÁT TRIỂN 10 NĂM                     |
+-------------------------------------------------------------------------------+
| Giới hạn hiện tại:                                                            |
| [!] Mô hình CFD dừng ở phân tích Cold Flow (chưa ghép nối mô hình cháy phun)  |
| [!] Góc phối khí và độ nâng xúpap cố định theo biên dạng cam nguyên bản       |
| [!] Nhiên liệu thử nghiệm giới hạn ở Diesel khoáng tiêu chuẩn                 |
|                                                                               |
| Định hướng tương lai (Agenda 2026 - 2035):                                   |
| [*] Giai đoạn 1: Ghép nối mô hình phản ứng cháy đa pha (Spray & Combustion)   |
| [*] Giai đoạn 2: Tối ưu hóa đồng thời hệ thống cam nạp biến thiên (VVT)       |
| [*] Giai đoạn 3: Mở rộng thử nghiệm trên nhiên liệu sinh học Biodiesel B20-B100|
| [*] Giai đoạn 4: Ứng dụng Deep Reinforcement Learning trong tối ưu hóa khí động|
+-------------------------------------------------------------------------------+

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận các giới hạn học thuật:

  1. Phạm vi mô hình dòng chảy: Quá trình mô phỏng CFD mới tập trung vào giai đoạn nạp - nén dòng lạnh (Cold Flow Analysis) với môi chất không khí thuần túy, chưa tích hợp mô hình phun tơi nhiên liệu đa pha (Lagrangian Discrete Phase Model), bay hơi và phản ứng động học hóa học của quá trình cháy thực tế.
  2. Ràng buộc cơ cấu phối khí: Nghiên cứu giữ nguyên quy luật nâng xúpap và góc pha phối khí của trục cam nguyên bản nhằm đảm bảo tính lắp lẫn và tiết kiệm chi phí hoán cải dây chuyền sản xuất của doanh nghiệp.
  3. Giới hạn nhiên liệu: Toàn bộ dữ liệu thực nghiệm mới thực hiện trên nhiên liệu dầu Diesel khoáng truyền thống, chưa mở rộng trên các dòng nhiên liệu thay thế như Biodiesel (B5, B10, B20) hoặc cồn sinh học.

Chương trình nghiên cứu tiếp nối giai đoạn 10 năm tới bao gồm:

  • Mở rộng mô hình CFD ghép nối đầy đủ giữa dòng nạp và quá trình cháy - phát thải động học ($NO_x$ / Soot kinetic mechanisms) trên phần mềm ANSYS Forte hoặc CONVERGE CFD.
  • Tối ưu hóa đồng thời biên dạng cam nạp và họng nạp sử dụng mạng nơ-ron sâu kết hợp học tăng cường (Deep Reinforcement Learning).
  • Nghiên cứu thực nghiệm tính năng phát thải độc hại theo chu trình thử 8 chế độ ISO 8178 E2/E3 trên các loại nhiên liệu tái tạo và hỗn hợp nhiên liệu sinh học.

Tác động và ảnh hưởng

+-------------------------------------------------------------------------------+
|                         MA TRẬN LAN TỎA ẢNH HƯỞNG                             |
+-------------------------------------------------------------------------------+
| Lĩnh vực tác động   | Kết quả lượng hóa / Chỉ số đo lường                     |
+---------------------+---------------------------------------------------------+
| Học thuật quốc tế   | 04 bài báo ISI/Scopus & KEM, 04 bài báo tạp chí chuyên  |
| & trong nước        | ngành Giao thông Vận tải, KH&CN ĐHQG, Thủy lợi-Môi trường|
| Chuyển giao công nghệ| Triển khai tại SVEAM, sản lượng tác động >45.000 máy/năm|
| Thị trường xuất khẩu| Nâng cao khả năng tiếp cận 20+ thị trường quốc tế       |
| Môi trường xã hội   | Giảm 5,58% tiêu thụ nhiên liệu => Giảm hàng chục ngàn   |
|                     | tấn CO2 phát thải hàng năm trong nông nghiệp            |
+-------------------------------------------------------------------------------+

Công trình của NCS. Lê Việt Hùng tạo ra chuỗi giá trị tác động sâu rộng:

  • Ảnh hưởng học thuật: Các công bố quốc tế trên Journal of Key Engineering Materials (KEM), Hội nghị Quốc tế ICFMAS và GTSD, cùng các tạp chí đầu ngành trong nước đã khẳng định năng lực tự chủ nghiên cứu khoa học công nghệ đỉnh cao của các nhà khoa học Việt Nam trong lĩnh vực kỹ thuật cơ khí động lực.
  • Tái cấu trúc hoạt động R&D doanh nghiệp: Chuyển đổi căn bản phương thức nghiên cứu phát triển sản phẩm tại VEAM/SVEAM từ phương pháp "kinh nghiệm thử - sai" lạc hậu sang quy trình "số hóa - mô phỏng tối ưu - kiểm nghiệm thực chứng", giảm chu kỳ phát triển sản phẩm mới từ 18 tháng xuống dưới 6 tháng.
  • Lợi ích kinh tế - môi trường: Với hàng trăm ngàn động cơ RV165-2 đang lưu hành trên các cánh đồng và tàu thuyền nhỏ, việc giảm $11,5\ g/HP.h$ suất tiêu hao nhiên liệu mang lại giá trị tiết kiệm kinh tế trực tiếp hàng trăm tỷ đồng cho nông - ngư dân, đồng thời cắt giảm phát thải lượng lớn $CO_2$ và khói đen vào khí quyển.

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

  • Nghiên cứu sinh và học viên cao học ngành Kỹ thuật Động lực & Cơ khí: Tiếp cận một tài liệu tham khảo mẫu mực về cấu trúc luận án, quy trình tham số hóa hình học và phương pháp tích hợp công cụ AI (ANN-DE) vào bài toán kỹ thuật thực tế.
  • Giảng viên và các nhà khoa học chuyên ngành Động cơ đốt trong: Khai thác hệ thống dữ liệu số hóa, mã nguồn tự động hóa kịch bản TUI/C#/MATLAB và các phát hiện lý thuyết để phục vụ giảng dạy các môn học Động cơ đốt trong nâng cao, Động lực học lưu chất tính toán.
  • Kỹ sư R&D tại các doanh nghiệp chế tạo máy nông nghiệp (VEAM, SVEAM, Vikyno): Ứng dụng ngay các thông số hình học tối ưu và quy trình công nghệ đúc để nâng cấp dòng sản phẩm hiện hữu và phát triển các thế hệ động cơ mới.
  • Các nhà hoạch định chính sách cơ giới hóa nông nghiệp và quản lý môi trường: Có thêm căn cứ khoa học xác đáng để xây dựng lộ trình áp dụng tiêu chuẩn phát thải và định mức tiết kiệm năng lượng cho phương tiện cơ giới nông thôn.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc thiết lập và giải thích bản chất vật lý của mối quan hệ phi tuyến giữa Hệ số nạp thể tích ($\eta_v$) và Tỷ số xoáy ($S_r$) trong không gian tham số hóa 3D của họng nạp dạng xoắn ốc. Nghiên cứu đã mở rộng Lý thuyết Dòng xoáy rối trong động cơ đốt trong của Heywood (1988) bằng cách chứng minh rằng: sự suy giảm hệ số nạp không phải là hệ quả tất yếu của việc tăng cường độ xoáy, mà là do sự phát triển mất kiểm soát của các vùng xoáy cục bộ gây nghẽn dòng; việc làm mượt đường sinh và tối ưu hóa tỷ số $T/R$ sẽ duy trì được cả hai chỉ tiêu khí động ở mức cực đại.

2. Điểm đổi mới phương pháp luận khi so sánh với ít nhất hai nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của S. Sanap (2013) (vốn chỉ khảo sát CFD ở trạng thái tĩnh với mục tiêu duy nhất là tỷ số xoáy $S_r = 1,8$) và nghiên cứu của Padmesh Mandloi & Gunjan Verma (2015) (chỉ khảo sát 3 thông số hình học đơn giản trên mô hình tĩnh), luận án này đã đổi mới vượt bậc khi:

  1. Xây dựng mô hình động 3D lưới chuyển động (Dynamic Mesh Cold Flow) mô phỏng chính xác hành trình 720° trục khuỷu.
  2. Tích hợp chuỗi khép kín tự động từ trích xuất tham số CAD -> giải thuật CFD -> hàm xấp xỉ ANN -> tối ưu hóa toàn cục DE, loại bỏ hoàn toàn yếu tố phỏng đoán chủ quan.

3. Phát hiện bất ngờ nhất trong quá trình phân tích dữ liệu là gì?

Phát hiện bất ngờ nhất là tại dải tốc độ cao (2.200 - 2.400 vòng/phút), việc giảm tỷ số chiều cao họng trên bề rộng góc xoắn ($T/R$) về vùng giá trị $1,28 - 1,35$ không những không làm tăng lực cản ma sát thành mà ngược lại còn làm tăng lưu lượng không khí thực tế qua xúpap thêm 5,3% nhờ việc hình thành dòng dẫn hướng bám sát vách buồng cháy, làm giảm áp suất ngược tại cửa nạp.

4. Luận án có cung cấp đầy đủ quy trình và công cụ để tái lập nghiên cứu (Replication Protocol) không?

Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình tái lập thông qua 13 phần phụ lục chuyên sâu bao gồm:

  • Toàn bộ bản vẽ kỹ thuật 2D/3D kích thước tham số hóa nắp máy.
  • Hệ thống mã nguồn kịch bản tự động hóa bằng C#, JavaScript và Text User Interface (TUI) cho Ansys-Fluent.
  • Mã nguồn lập trình mạng nơ-ron nhân tạo ANN và thuật toán tiến hóa vi phân DE trên môi trường MATLAB.
  • Giấy chứng nhận hiệu chuẩn thiết bị đo và kết quả đo kiểm độc lập của Trung tâm QUATEST 3.

5. Lộ trình nghiên cứu phát triển 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình 10 năm (2026 - 2035) được định hình qua 3 giai đoạn:

  • Giai đoạn ngắn hạn (1-3 năm): Mở rộng mô hình CFD sang giai đoạn cháy cưỡng bức có phun tơi nhiên liệu đa pha và đánh giá phát thải $NO_x$, bồ hóng theo chu trình ISO 8178.
  • Giai đoạn trung hạn (3-5 năm): Thiết kế tối ưu hóa hệ thống phối khí biến thiên thông minh (VVT) tích hợp họng nạp cải tiến trên các dòng động cơ nhiều xylanh.
  • Giai đoạn dài hạn (5-10 năm): Chuyển đổi và tương thích hóa toàn bộ buồng cháy và hệ thống nạp cho các dòng nhiên liệu sạch thế hệ mới (E-fuels, Hydrogen và Biodiesel nguyên chất).

Kết luận

  1. Luận án đã giải quyết trọn vẹn và xuất sắc bài toán khoa học kỹ thuật cấp thiết về nâng cao tính năng làm việc của động cơ Diesel VIKYNO RV165-2 thông qua cải tiến toàn diện cụm họng nạp.
  2. Xây dựng thành công phương pháp tham số hóa toán học chính xác biên dạng hình học phức tạp của họng nạp xoắn ốc bên trong và cổ hút bên ngoài nắp xylanh.
  3. Thiết lập một quy trình tự động hóa mô phỏng CFD 3 chiều dòng không dừng tích hợp mô hình xấp xỉ mạng nơ-ron nhân tạo (ANN) và giải thuật tối ưu hóa tiến hóa vi phân (DE) có tính chính xác cao, hội tụ nhanh.
  4. Xác lập cấu trúc hình học tối ưu giúp tăng hệ số nạp $\eta_v$ thêm 8,95%, tăng tỷ số xoáy $S_r$ thêm 25,35%, từ đó tăng công suất động cơ thêm 5,71% và giảm suất tiêu hao nhiên liệu 5,58% tại chế độ công suất định mức.
  5. Toàn bộ kết quả mô phỏng số đã được kiểm chứng và khẳng định tính chính xác tuyệt đối thông qua thực nghiệm đo kiểm trên băng thử công nghiệp đạt chuẩn đo lường quốc gia.
  6. Mở ra một hướng tiếp cận công nghệ mới, nâng tầm năng lực tự chủ thiết kế và sản xuất động cơ đốt trong của nền cơ khí chế tạo máy Việt Nam hội nhập vững chắc vào chuỗi giá trị toàn cầu.