Tổng quan về luận án

Trước bối cảnh cạn kiệt nguồn nhiên liệu hóa thạch toàn cầu và sức ép giảm phát thải khí nhà kính (đặc biệt là $\text{CO}_2$, $\text{NO}_x$ và muội than) theo các quy chuẩn quốc tế khắt khe như ISO 8178, ngành công nghiệp động cơ đốt trong phục vụ nông – lâm – ngư nghiệp tại các quốc gia đang phát triển đối mặt với yêu cầu cấp thiết về nâng cao hiệu suất nhiệt và tối ưu hóa quá trình nhiệt - khí động. 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 cung cấp hơn 45.000 động cơ diesel cỡ nhỏ mỗi năm cho thị trường nội địa và xuất khẩu đến hơn 20 quốc gia. Tuy nhiên, sản phẩm chủ lực là động cơ diesel một xilanh VIKYNO RV165-2 (công suất tối đa 16,5 HP tại 2.400 vòng/phút, công suất định mức 14 HP tại 2.200 vòng/phút, dung tích công tác $V_h = 839\text{ cm}^3$, đường kính $\times$ hành trình piston $105 \times 97\text{ mm}$, tỷ số nén $\varepsilon = 20$, phun trực tiếp DI với áp suất mở vòi phun $220\text{ kg/cm}^2$, suất tiêu hao nhiên liệu định mức $206\text{ g/HP.h}$) vốn được phát triển dựa trên nền tảng kinh nghiệm kế thừa từ các mẫu nguyên bản của Kubota và Yanmar.

Điểm nghẽn học thuật và công nghệ (research gap) cốt lõi nằm ở chỗ: các cải tiến đường nạp trước đây chủ yếu dựa trên phương pháp thử nghiệm mò mẫm (trial-and-error), tiếp cận một chiều thông qua mô phỏng 1D đơn giản hóa hoặc chỉ can thiệp vào đường ống nối bên ngoài nắp xilanh, hoàn toàn thiếu vắng một quy trình khép kín tích hợp 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 3D với các giải thuật trí tuệ nhân tạo để tối ưu hóa biên dạng cụm họng nạp xoắn ốc (helical intake port) nằm sâu bên trong nắp xilanh.

Để giải quyết triệt để bài toán này, luận án của NCS. Lê Việt Hùng (chuyên ngành Cơ kỹ thuật, mã số 12252010105, Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh, dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Anh Thi và PGS.TS. Đỗ Văn Dũng) tập trung giải quyết các câu hỏi nghiên cứu và giả thuyết khoa học mang tính nền tảng:

  • Câu hỏi nghiên cứu 1 (RQ1): Biên dạng hình học không gian 3D của cụm họng nạp xoắn ốc bên trong nắp xilanh chi phối như thế nào đến sự hình thành trường vận tốc, trường áp suất và cấu trúc dòng xoáy (swirl/tumble) trong kỳ nạp - nén của động cơ diesel 1 xilanh?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để thiết lập một mô hình toán học xấp xỉ chính xác cao giữa các thông số hình học họng nạp với hệ số nạp ($\eta_v$) và hệ số xoáy ($S_r$, $T_{\text{ratio}}$) nhằm vượt qua sự tốn kém tài nguyên tính toán của mô phỏng CFD truyền thống?
  • Câu hỏi nghiên cứu 3 (RQ3): Giải thuật tối ưu hóa tiến hóa vi phân (DE) kết hợp mạng nơ-ron nhân tạo (ANN) có thể tìm ra điểm thỏa hiệp tối ưu Pareto giữa lưu lượng nạp cực đại và cường độ xoáy lốc lý tưởng hay không?
  • Giả thuyết khoa học (H1): Tối ưu hóa tham số hình học cụm họng nạp kết hợp cổ nối ngoài sẽ đồng thời nâng cao hệ số nạp $\eta_v$ và hệ số xoáy $S_r$, từ đó tăng công suất có ích ($N_e$), gia tăng moment xoắn ($M_e$) và giảm rõ rệt suất tiêu hao nhiên liệu ($g_e$) trên toàn dải đặc tính ngoài của động cơ.

Khung lý thuyết của nghiên cứu được xây dựng dựa trên cơ sở nhiệt động lực học động cơ đốt trong (Internal Combustion Engine Theory của Heywood và Taylor), cơ học lưu chất tính toán (Computational Fluid Dynamics - RANS Navier-Stokes), lý thuyết mạng nơ-ron nhân tạo (Artificial Neural Network - MLP) và giải thuật tối ưu hóa tiến hóa (Differential Evolution Algorithm). Nghiên cứu tạo ra bước đột phá khi lần đầu tiên số hóa, tham số hóa hoàn chỉnh cụm họng nạp động cơ diesel RV165-2 bằng 13 biến số hình học độc lập, tự động hóa chu trình mô phỏng ANSYS Fluent - MATLAB - DE, và kiểm chứng thực nghiệm trên băng thử phanh điện động lực học APA 100 cùng thiết bị đo tiêu hao nhiên liệu tự động AVL 733S. Kết quả được kiểm định độc lập bởi Trung tâm Kỹ thuật Tiêu chuẩn Đo lường Chất lượng 3 (QUATEST 3), chứng minh tính tiên phong và độ tin cậy khoa học tuyệt đối.

Literature Review và Positioning

Quá trình nạp của động cơ đốt trong 4 kỳ là một quá trình khí động học không dừng (unsteady), chịu sự chi phối mạnh mẽ bởi quy luật biến thiên tiết diện lưu thông của xúpap và đặc trưng hình học ba chiều phức tạp của họng nạp. Hiệu quả của quá trình này được định lượng thông qua hệ số nạp $\eta_v$, xác định theo biểu thức:

$$\eta_v = \frac{m_a}{\rho_{a,i} V_h} = \frac{2 \dot{m}a}{\rho{a,i} V_h n} = \frac{4 \dot{m}a}{\rho{a,i} A_p S_p}$$

Trong đó, $m_a$ là khối lượng không khí nạp thực tế vào xilanh trong một chu trình; $\dot{m}a$ là lưu lượng khối lượng không khí trung bình; $\rho{a,i}$ là khối lượng riêng không khí trước họng nạp; $V_h$ là thể tích công tác; $n$ là số vòng quay trục khuỷu; $A_p$ là diện tích đỉnh piston; $S_p$ là vận tốc trung bình của piston.

Tổng quan y văn quốc tế ghi nhận nhiều nỗ lực nghiên cứu dòng chảy qua họng nạp động cơ. Puzinauskas và Midkiff khi nghiên cứu tối ưu hóa họng nạp cho động cơ SI bằng phần mềm KIVA và STAR-CD kết hợp tiền xử lý ICEM CFD đã chỉ ra rằng việc triệt tiêu dòng xoáy cục bộ không mong muốn ở đỉnh họng nạp giúp gia tăng độ rối và nâng mức độ bay hơi nhiên liệu từ 68,9% lên 90,1% tại 540° góc quay trục khuỷu (GQTK). Tuy nhiên, công trình này chưa xây dựng được hàm mục tiêu toán học tường minh cho bài toán tối ưu đa biến.

Nhóm nghiên cứu tại D.Y. Patil College of Engineering (Ấn Độ) sử dụng mô phỏng CFD khảo sát họng nạp đôi dạng xoắn ốc và đối chứng thực nghiệm trên băng thổi khí tĩnh, thu được 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 khi độ nâng xúpap tăng từ 1 mm đến 8 mm. Dẫu vậy, hạn chế cốt tử là nghiên cứu chỉ khảo sát ở các độ nâng xúpap tĩnh, rời rạc, không phản ánh tính liên tục theo thời gian thực của chu trình động cơ.

Trong một nghiên cứu khác, S. Sanap tập trung thiết kế họng nạp xoắn ốc cho động cơ diesel nhằm đạt hệ số xoáy mục tiêu 1,8, chỉ ra mối quan hệ nghịch đảo giữa tỷ số $T/R$ (chiều cao họng nạp $T$ trên chiều rộng góc xoắn $R$) với cường độ xoáy: khi giảm tỷ số $T/R$, luồng khí bám sát thành xilanh tạo ra xoáy lốc mạnh hơn, nhưng nghiên cứu này lại bỏ qua sự suy giảm nghiêm trọng của hệ số nạp $\eta_v$. Padmesh Mandloi và Gunjan Verma khảo sát 3 thông số hình học (góc xoắn 50°–70°, bán kính xoắn 30–50 mm, chiều dài họng nạp 50–65 mm) bằng ANSYS Fluent và thiết lập biểu đồ độ nhạy cục bộ (local sensitivity charts), song mô hình hình học còn quá đơn giản hóa. Jun-ichi Kawashima nghiên cứu động cơ diesel 4 xúpap phun trực tiếp tốc độ cao, so sánh hiệu ứng xoáy giữa họng nạp tiếp tuyến và họng nạp xoắn ốc, xác nhận sự cạnh tranh gay gắt giữa hệ số lưu lượng và tỷ số xoáy.

Tại Việt Nam, Nguyễn Hữu Hường và Vương Như Long (Đại học Bách Khoa TP.HCM) sử dụng phần mềm AVL BOOST 1D để cải tiến động cơ RV195, đạt mức tăng hệ số nạp 8,92% và tăng công suất 10,2%, nhưng mô hình 1D hoàn toàn bất lực trong việc mô tả trường xoáy không gian 3D. Võ Danh Toàn, Nguyễn Thanh Tuấn, Lê Việt Hùng et al. (2015) khảo sát cải tiến họng nạp RV165-2 bằng AVL BOOST nhưng mới dừng lại ở đường ống bên ngoài nắp xilanh và thiếu kiểm chứng thực nghiệm.

Vị thế học thuật (positioning) của luận án được xác lập rõ ràng: Đây là công trình tiên phong giải quyết cuộc tranh luận kinh điển giữa trường phái "ưu tiên lưu lượng nạp cực đại" và "ưu tiên cường độ xoáy lốc phục vụ hòa trộn" bằng cách tích hợp mô phỏng động lực học 3D Cold Flow trong ANSYS ICE/Fluent với mạng nơ-ron nhân tạo (ANN) và giải thuật tiến hóa vi phân (DE). Nghiên cứu không chỉ tối ưu hóa riêng lẻ từng bộ phận mà thực hiện cải tiến đồng bộ từ cổ nối lọc gió bên ngoài đến họng xoắn ốc 3D bên trong nắp xilanh cho động cơ diesel một xilanh cỡ nhỏ.

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

Đó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 chuyển động xoáy và lý thuyết nạp của Heywood và Taylor trong điều kiện dòng chảy nén được, không dừng (unsteady compressible flow) qua cụm xúpap nạp có biên dạng phức tạp. Nghiên cứu chứng minh rằng mối quan hệ giữa hệ số nạp ($\eta_v$) và hệ số xoáy ngang ($T_{\text{ratio}}$)/xoáy dọc ($S_r$) không phải là mối quan hệ tuyến tính đơn điệu hay đối kháng tuyệt đối, mà là một bề mặt đáp ứng phi tuyến bậc cao có các vùng cực trị cục bộ.

Bằng cách phân tích định lượng trường biến dạng trung bình $E_{ij}$, độ nhớt rối $\mu_t$ và tốc độ tiêu tán động năng rối $\varepsilon$, luận án làm sáng tỏ cơ chế vật lý: sự thay đổi độ dốc của vách xoắn ốc và góc nghiêng của đường sinh họng nạp có thể tái định hướng các vector vận tốc dòng khí, tạo ra chuyển động xoáy lốc cưỡng bức mạnh mà không làm phát sinh hiện tượng nghẽn dòng (flow choking) hay gia tăng tổn thất áp suất cục bộ $\Delta p_a$.

Mô hình lý thuyết mới được đề xuất cho phép dự báo chính xác diễn biến của quá trình nạp - nén thông qua hệ thống phương pháp luận lai (hybrid CFD - AI paradigm), thay thế các công thức kinh nghiệm bán thực nghiệm vốn có sai số lớn đối với động cơ diesel cỡ nhỏ có buồng cháy thống nhất trên đỉnh piston.

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

Khung phân tích của luận án là một bước tiến mang tính tích hợp liên ngành sâu sắc giữa Cơ học lưu chất tính toán (CFD), Trí tuệ nhân tạo (Mạng nơ-ron nhân tạo - ANN) và Toán học tối ưu hóa toàn cục (Giải thuật tiến hóa vi phân - DE):

+-------------------------------------------------------------------------+
|                  1. THAM SỐ HÓA HÌNH HỌC 3D (CAD/SolidWorks)            |
|       13 thông số độc lập: Chiều cao T, Bề rộng R, Tỷ số T/R, Góc xoắn  |
+-------------------------------------------------------------------------+
+-------------------------------------------------------------------------+
|             2. TỰ ĐỘNG HÓA MÔ PHỎNG CFD 3D (ANSYS ICE / Fluent)         |
|  Phương trình RANS, Mô hình rối k-ε, Dynamic Mesh (0° - 720° GQTK)       |
|  Xuất dữ liệu: Hệ số nạp (η_v), Hệ số xoáy (S_r, T_ratio), Trường P, V  |
+-------------------------------------------------------------------------+
+-------------------------------------------------------------------------+
|                3. XẤP XỈ BỀ MẶT ĐÁP ỨNG BẰNG MẠNG NƠ-RON (ANN)          |
|  Kiến trúc MLP (13 đầu vào - 10 nơ-ron ẩn - 2 đầu ra: η_v, S_r)          |
|  Huấn luyện Backpropagation, 100 điểm dữ liệu, 1.000 epochs, MSE/MAPE   |
+-------------------------------------------------------------------------+
+-------------------------------------------------------------------------+
|             4. TỐI ƯU HÓA TOÀN CỤC BẰNG TIẾN HÓA VI PHÂN (DE)           |
|  Khởi tạo quần thể -> Đột biến (F) -> Lai tạo (CR) -> Chọn lọc Pareto   |
|  Tìm bộ tham số tối ưu thỏa mãn: Maximize(η_v) & Target(S_r)            |
+-------------------------------------------------------------------------+
+-------------------------------------------------------------------------+
|           5. CHẾ TẠO MẪU NẮP XILANH MỚI & THỰC NGHIỆM ĐỐI CHỨNG         |
|  Băng thử APA 100, Cân 733S, Cảm biến áp suất, Kiểm định QUATEST 3      |
+-------------------------------------------------------------------------+

Khung phân tích này thiết lập điều kiện biên xác định (boundary conditions): áp suất cửa vào $p_0 = 101.325\text{ Pa}$, nhiệt độ môi trường $T_0 = 300\text{ K}$, quy luật nâng hạ xúpap nạp/thải theo góc quay trục khuỷu thực tế từ $0^\circ$ đến $720^\circ$ GQTK, bỏ qua hiện tượng truyền nhiệt qua vách họng nạp trong mô hình Cold Flow để tập trung giải quyết tương tác khí động thuần túy.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ thế giới quan thực chứng luận (positivism) kết hợp chặt chẽ với chủ nghĩa hiện thực phê phán (critical realism), đặt trọng tâm vào việc lượng hóa chính xác các đại lượng vật lý thông qua các phương trình bảo toàn và kiểm chứng nghiêm ngặt bằng thực nghiệm vật lý. Thiết kế nghiên cứu là sự kết hợp đa tầng (multi-level design):

  • Tầng 1 (Cấp độ vĩ mô - Ngoại vi): Khảo sát thực nghiệm và tối ưu hóa hình học đường dẫn khí nạp bên ngoài nắp xilanh (cổ nối bộ lọc gió) với 4 phương án thiết kế khác nhau.
  • Tầng 2 (Cấp độ vi mô - Nội tại): Tham số hóa 3D cụm họng nạp xoắn ốc bên trong nắp xilanh với 13 biến hình học chi tiết, giải bài toán động lực học lưu chất không dừng 3D.
  • Tầng 3 (Cấp độ hệ thống): Tổng hợp tối ưu hóa toàn diện họng nạp trong và ngoài, đúc tạo phôi nắp xilanh thực tế, gia công chính xác và đo kiểm toàn bộ đặc tính ngoài của động cơ.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thiết kế với độ chuẩn xác học thuật cao:

Tham số hóa hình học

Số hóa biên dạng nắp xilanh hiện hữu bằng máy quét 3D và phần mềm SolidWorks. Xác định 13 tham số hình học trọng yếu điều khiển đường sinh xoắn ốc: bán kính xoắn $R$, chiều cao họng nạp $T$, góc định vị tạo xoáy $\psi$, tọa độ tâm họng nạp $(X, Y)$, góc nghiêng cửa vào $\theta$, và độ dốc thành xoắn $\beta$.

Mô hình hóa toán học lưu chất

Giải hệ phương trình vi phân đạo hàm riêng phi tuyến bảo toàn khối lượng, động lượng (Navier-Stokes) và năng lượng:

$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \vec{u}) = 0$$

$$\frac{\partial (\rho \vec{u})}{\partial t} + \nabla \cdot (\rho \vec{u} \vec{u}) = -\nabla p + \nabla \cdot (\overline{\bar{\tau}}) + \rho \vec{g}$$

Sử dụng mô hình rối hai phương trình $k-\varepsilon$ tiêu chuẩn (Standard $k-\varepsilon$ Turbulence Model) với các phương trình truyền tải động năng rối $k$ và tỷ lệ tiêu tán rối $\varepsilon$:

$$\frac{\partial (\rho k)}{\partial t} + \frac{\partial (\rho k u_i)}{\partial x_i} = \frac{\partial}{\partial x_j} \left[ \left( \mu + \frac{\mu_t}{\sigma_k} \right) \frac{\partial k}{\partial x_j} \right] + G_k - \rho \varepsilon$$

Thiết lập mô phỏng tự động hóa

Xây dựng chuỗi kịch bản điều khiển tự động kết hợp mã nguồn JavaScript, C# UDFs (User-Defined Functions), Text User Interface (TUI) và MATLAB để tự động hóa hoàn toàn các bước: xuất file CAD $\rightarrow$ chia lưới động (dynamic layering và remeshing trong module ANSYS ICE) $\rightarrow$ gán điều kiện biên chuyển động piston/xúpap $\rightarrow$ giải hội tụ $\rightarrow$ trích xuất dữ liệu $\eta_v, S_r, T_{\text{ratio}}$.

Data và phân tích

Bộ dữ liệu nền tảng gồm 100 trường hợp hình học biến thiên được mô phỏng CFD chi tiết. Dữ liệu này được đưa vào huấn luyện mạng nơ-ron nhân tạo đa tầng (MLP) với thuật toán lan truyền ngược (Backpropagation):

  • Cấu trúc mạng: 13 nơ-ron lớp đầu vào, 10 nơ-ron lớp ẩn (hàm kích hoạt Sigmoid) và 2 nơ-ron lớp đầu ra ($\eta_v$ và $S_r$).
  • Quá trình huấn luyện qua 1.000 vòng lặp (epochs) đạt độ hội tụ xuất sắc với sai số toàn phương trung bình $\text{MSE} = 1,24 \times 10^{-4}$ và sai số phần trăm tuyệt đối trung bình $\text{MAPE} < 1,5%$ trên tập dữ liệu kiểm định chéo (cross-validation).

Sau khi mạng ANN xấp xỉ thành công bề mặt đáp ứng, thuật toán Tiến hóa vi phân (DE) được kích hoạt với quy mô quần thể $NP = 50$, hệ số đột biến $F = 0,8$, xác suất lai tạo $CR = 0,7$. Sau 150 thế hệ tiến hóa, giải thuật đã hội tụ về nghiệm tối ưu toàn cục, xác định chính xác bộ kích thước hình học lý tưởng cho cụm họng nạp mới.

Hệ thống thực nghiệm đối chứng được lắp đặt tại Trung tâm R&D của SVEAM:

  • Băng thử công suất phanh điện APA 100 (công suất hấp thụ tối đa 100 kW, sai số đo moment $< \pm 0,2%$).
  • Thiết bị cân đo tiêu hao nhiên liệu tự động AVL 733S (đo lưu lượng theo phương pháp trọng lượng vi sai với độ chính xác $\pm 0,1%$).
  • Cảm biến áp suất khí nạp đáp ứng tần số cao TOYOTA 89420-20300 kết hợp bộ thu thập dữ liệu đa kênh và cảm biến tiệm cận quang học đo chính xác góc quay trục khuỷu với độ phân giải $0,1^\circ$ GQTK.
  • Băng đo dòng thổi khí chuyên dụng đo kiểm chứng hệ số nạp tĩnh và động.

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

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

Kết quả nghiên cứu mang lại những phát hiện đột phá có giá trị khoa học và thực tiễn sâu sắc:

+-------------------------------------------------------------------------------------------------------+
|             BẢNG SO SÁNH CÁC CHỈ TIÊU KỸ THUẬT then chốt TRƯỚC VÀ SAU CẢI TIẾN                        |
+------------------------------------+--------------------+--------------------+------------------------+
| Thông số đánh giá                  | Động cơ hiện hữu   | Động cơ cải tiến   | Mức độ cải thiện (%)   |
+------------------------------------+--------------------+--------------------+------------------------+
| Hệ số nạp thực nghiệm (η_v)        | 0,812              | 0,874              | + 7,64%                |
| Công suất tối đa (N_e max)         | 16,5 HP / 2400 rpm | 17,62 HP / 2400 rpm| + 6,79%                |
| Công suất định mức (N_e đm)        | 14,0 HP / 2200 rpm | 15,08 HP / 2200 rpm| + 7,71%                |
| Moment xoắn cực đại (M_e max)      | 4,90 kg.m          | 5,28 kg.m          | + 7,76%                |
| Suất tiêu hao nhiên liệu (g_e đm)  | 206 g/HP.h         | 194,5 g/HP.h       | - 5,58% (Tiết kiệm)   |
| Hệ số xoáy ngang (T_ratio)         | 1,42               | 1,78               | + 25,35%               |
+------------------------------------+--------------------+--------------------+------------------------+
  • Đột phá về cấu trúc dòng chảy và trường vận tốc: Biên dạng họng nạp xoắn ốc cải tiến đã triệt tiêu hoàn toàn vùng xoáy quẩn tách dòng (flow separation) xuất hiện ở khúc uốn gắt của nắp xilanh nguyên bản. Trường vận tốc tại tiết diện xúpap phân bố đồng đều hơn, vận tốc dòng cực đại đạt $92\text{ m/s}$ tại thời điểm xúpap mở lớn nhất mà không gây nghẽn dòng.
  • Tương quan phi tuyến giữa hệ số nạp và hệ số xoáy: Phân tích cụm dữ liệu (K-means clustering với $k=5$) từ 100 bộ tham số chứng minh rằng việc tăng hệ số xoáy không nhất thiết làm giảm hệ số nạp nếu góc nghiêng vách dẫn hướng $\beta$ và bán kính cung lượn họng xúpap được tối ưu hóa đồng thời.
  • Tính ưu việt của giải pháp tích hợp: Cải tiến độc lập cổ nối ngoài nắp xilanh đóng góp tăng 2,1% hệ số nạp; cải tiến họng nạp trong nắp xilanh đóng góp tăng 5,54%; khi kết hợp đồng bộ cả hai, hệ số nạp tổng thể tăng vọt 7,64%, tạo hiệu ứng cộng hưởng khí động học rõ rệt.

Implications đa chiều

Ý nghĩa lý thuyết và học thuật

Công trình xác lập một phương pháp luận chuẩn mực trong việc kết hợp CFD 3D với Mạng nơ-ron nhân tạo và Giải thuật tiến hóa vi phân (CFD-ANN-DE Framework). Khung lý thuyết này có khả năng tổng quát hóa cao, mở đường cho việc ứng dụng trên các loại động cơ đốt trong đa xilanh, động cơ nhiên liệu kép (dual-fuel) hoặc động cơ sử dụng nhiên liệu sinh học (biogas, biodiesel).

Ứng dụng thực tiễn trong công nghiệp chế tạo

Cung cấp bộ bản vẽ thiết kế chế tạo nắp xilanh hoàn chỉnh và khuôn đúc kim loại đã được tối ưu hóa cho SVEAM. Doanh nghiệp có thể đưa trực tiếp vào dây chuyền sản xuất hàng loạt mà không cần đầu tư thay đổi toàn bộ kết cấu thân máy hiện hữu.

Tác động môi trường và chính sách

Việc giảm suất tiêu hao nhiên liệu $11,5\text{ g/HP.h}$ (tương đương giảm 5,58%) trực tiếp cắt giảm hàng ngàn tấn phát thải $\text{CO}_2$, $\text{CO}$ và hydrocacbon chưa cháy (HC) ra môi trường nông nghiệp mỗi năm, hỗ trợ đắc lực cho các mục tiêu chuyển đổi xanh và đạt chuẩn khí thải theo lộ trình của Chính phủ.

Limitations và Future Research

Dù đạt được những bước tiến đột phá, luận án thẳng thắn ghi nhận 4 giới hạn học thuật cần tiếp tục mở rộng:

  1. Giới hạn mô hình dòng chảy lạnh (Cold Flow Assumption): Mô phỏng CFD trong luận án sử dụng mô hình dòng chảy lạnh, chưa kết hợp mô hình phun nhiên liệu đa pha (discrete phase model - DPM) và mô hình phản ứng cháy hóa học chi tiết trong xilanh để đánh giá trực tiếp quá trình hình thành muội than (soot) và oxit nitơ ($\text{NO}_x$).
  2. Bỏ qua hiệu ứng biến dạng nhiệt (Thermal-Structural Coupling): Mô hình tính toán chưa xét đến ảnh hưởng của trường nhiệt độ cao gây biến dạng cơ học của họng nạp và xúpap trong điều kiện vận hành thực tế ở chế độ quá tải kéo dài.
  3. Phạm vi đối tượng nghiên cứu: Các thông số tối ưu mới chỉ thực hiện trên dòng động cơ diesel 1 xilanh nằm ngang làm mát bằng nước dung tích $839\text{ cm}^3$ (RV165-2), chưa tự động ngoại suy cho các dòng động cơ công suất nhỏ hơn (RV125, RV145) hoặc lớn hơn (RV195, RV220).
  4. Không gian tham số tối ưu hóa: Giải thuật DE mới chỉ khảo sát 13 thông số hình học tĩnh của họng nạp mà chưa kết hợp tối ưu hóa biến thiên thời gian đóng mở xúpap (Variable Valve Timing - VVT).

Chương trình nghiên cứu tương lai (Future Agenda):

  • Phát triển mô hình mô phỏng cháy toàn diện 3D (Combustion CFD) tích hợp động học phản ứng rút gọn để khảo sát phát thải chi tiết theo chu trình thử ISO 8178.
  • Ứng dụng giải thuật tối ưu hóa đa mục tiêu nâng cao (NSGA-III, MOEA/D) để đồng thời tối ưu hóa họng nạp và biên dạng buồng cháy trên đỉnh piston (piston bowl geometry).
  • Nghiên cứu ứng dụng họng nạp cải tiến khi động cơ chuyển đổi sang sử dụng nhiên liệu sinh học B10, B20 và hỗn hợp khí sinh học Biogas/Diesel.

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

Luận án tạo ra tác động lan tỏa mạnh mẽ trên nhiều phương diện:

  • Ảnh hưởng học thuật: Thiết lập chuẩn mực nghiên cứu mới trong lĩnh vực Cơ kỹ thuật và Động cơ nhiệt tại Việt Nam. Các công bố quốc tế thuộc hệ thống ISI/Scopus (như Journal of Key Engineering Materials, GTSD, ICFMAS) cùng nhiều bài báo trên các tạp chí chuyên ngành uy tín trong nước khẳng định giá trị trích dẫn học thuật lâu dài.
  • Chuyển dịch công nghệ trong ngành cơ khí chế tạo: Giúp SVEAM và VEAM làm chủ công nghệ thiết kế lõi, chuyển từ vị thế gia công lắp ráp thụ động sang tự chủ thiết kế, gia tăng tỷ lệ nội địa hóa thực chất và nâng cao năng lực cạnh tranh trực diện với các thương hiệu Nhật Bản (Kubota, Yanmar) và Trung Quốc.
  • Lợi ích kinh tế - xã hội: Với sản lượng tiêu thụ hàng chục ngàn máy mỗi năm tại thị trường ĐBSCL, đồng bằng sông Hồng và xuất khẩu quốc tế, việc tiết kiệm 5,58% nhiên liệu diesel mang lại giá trị kinh tế trực tiếp hàng trăm tỷ đồng mỗi năm cho bà con nông dân và ngư dân, đồng thời đóng góp thiết thực vào an ninh năng lượng quốc gia.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một quy trình phương pháp luận chuẩn hóa về liên kết CAD-CFD-ANN-DE, có thể tái lập (replicate) và mở rộng cho các bài toán tối ưu hóa khí - thủy động phức tạp khác.
  • Kỹ sư R&D ngành Động cơ và Ô tô: Nhận chuyển giao trọn vẹn bộ công cụ tự động hóa mô phỏng bằng mã lệnh JavaScript, C# UDF và MATLAB, rút ngắn chu kỳ thiết kế sản phẩm mới từ hàng năm xuống còn vài tuần.
  • Doanh nghiệp sản xuất máy nông nghiệp (SVEAM, VEAM): Sở hữu giải pháp công nghệ có khả năng thương mại hóa ngay lập tức, nâng cao tính năng động cơ RV165-2 lên mức công suất 17,62 HP mà không phát sinh thêm chi phí vật tư sản xuất.
  • Nông dân và Ngư dân: Được sử dụng thế hệ máy móc động lực mạnh mẽ hơn, bền bỉ hơn, tiết kiệm chi phí nhiên liệu vận hành máy cày, máy gặt đập liên hợp, ghe xuồng và trạm bơm tưới tiêu.

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à chứng minh thực nghiệm mô hình tương tác phi tuyến đa chiều giữa cấu trúc dòng xoáy 3D không dừng với hiệu suất nạp thể tích trong động cơ diesel cỡ nhỏ. Công trình đã mở rộng lý thuyết quá trình nạp truyền thống của Heywood và Taylor, chỉ ra rằng việc điều khiển chính xác tỷ số $T/R$ và độ dốc thành xoắn $\beta$ trong họng nạp xoắn ốc có thể đồng thời tối ưu hóa hệ số nạp $\eta_v$ và hệ số xoáy $S_r$, phá vỡ quan niệm kinh điển cho rằng hai đại lượng này luôn triệt tiêu lẫn nhau.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây là gì?

So với nghiên cứu của Puzinauskas & Midkiff (chưa có hàm tối ưu toán học) hay D.Y. Patil College & S. Sanap (chỉ mô phỏng tĩnh ở độ nâng xúpap rời rạc), luận án đã:

  1. Xây dựng mô hình động 3D liên tục xuyên suốt kỳ nạp - nén ($0^\circ - 720^\circ$ GQTK) bằng lưới chuyển động trong ANSYS ICE.
  2. Tự động hóa hoàn toàn chu trình tính toán bằng kịch bản lập trình C#/JavaScript/MATLAB.
  3. Sử dụng mô hình xấp xỉ mạng nơ-ron nhân tạo kết hợp giải thuật tiến hóa vi phân DE để tìm điểm cực trị toàn cục, giảm 90% thời gian tính toán so với việc quét lưới tham số CFD truyền thống.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là sự đóng góp phi tuyến tính của góc uốn cổ nối lọc gió bên ngoài nắp xilanh. Việc chỉ thay đổi độ cong và tiết diện của chi tiết ngoài này đã giúp tăng ngay 2,1% hệ số nạp và khi kết hợp với nắp xilanh mới đã tạo ra bước nhảy vọt công suất định mức từ $14,0\text{ HP}$ lên $15,08\text{ HP}$ ($+7,71%$), vượt xa dự đoán ban đầu của các mô hình lý thuyết 1D thuần túy.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không?

Luận án cung cấp chi tiết toàn bộ quy trình tái lập tại hệ thống 13 Phụ lục, bao gồm:

  • Bản vẽ chế tạo cơ khí 2D/3D chi tiết nắp xilanh (Phụ lục 1, 3).
  • Toàn bộ mã nguồn kịch bản JavaScript, C# UDF, TUI điều khiển ANSYS Fluent (Phụ lục 4, 6, 7, 8, 9).
  • Toàn bộ mã nguồn MATLAB huấn luyện mạng ANN và giải thuật tối ưu hóa tiến hóa DE (Phụ lục 10, 13).
  • Hồ sơ kiểm chuẩn thiết bị đo và Giấy chứng nhận kết quả thử nghiệm độc lập của QUATEST 3 (Phụ lục 2, 12).

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm tới định hướng:

  • Giai đoạn 1 (1–3 năm): Mở rộng nền tảng CFD-ANN-DE cho toàn bộ dải động cơ RV từ 5 HP đến 30 HP của SVEAM.
  • Giai đoạn 2 (3–6 năm): Tích hợp tối ưu hóa quá trình cháy đa pha và phát thải độc hại ($\text{NO}_x$, PM) đáp ứng chuẩn Tier 4/Euro 5 cho máy nông nghiệp.
  • Giai đoạn 3 (6–10 năm): Phát triển hệ thống động cơ lai (Hybrid Diesel-Electric) và động cơ sử dụng nhiên liệu không carbon (Hydrogen/Ammonia) phục vụ cơ giới hóa nông nghiệp thông minh.

Kết luận

Luận án tiến sĩ của NCS. Lê Việt Hùng đã giải quyết xuất sắc và trọn vẹn một bài toán kỹ thuật công nghệ trọng điểm có tính thời sự và giá trị ứng dụng thực tiễn to lớn:

  1. Hoàn thiện phương pháp luận tham số hóa 3D cụm họng nạp xoắn ốc động cơ diesel một xilanh với 13 biến số hình học độc lập.
  2. Xây dựng thành công chuỗi công cụ tự động hóa mô phỏng khí động học 3D không dừng trong ANSYS Fluent bằng kịch bản lập trình chuyên sâu.
  3. Ứng dụng đột phá mô hình lai Trí tuệ nhân tạo (ANN) và Giải thuật tiến hóa vi phân (DE), giải quyết triệt để bài toán tối ưu hóa đa mục tiêu giữa hệ số nạp và hệ số xoáy.
  4. Chứng minh thực nghiệm thuyết phục: nâng cao hệ số nạp thêm 7,64%, tăng công suất định mức thêm 7,71% (từ 14 HP lên 15,08 HP), tăng moment xoắn cực đại thêm 7,76% (từ 4,9 kg.m lên 5,28 kg.m) và giảm suất tiêu hao nhiên liệu 5,58% (từ 206 xuống 194,5 g/HP.h).
  5. Mở ra 3 hướng nghiên cứu mới: tối ưu hóa buồng cháy động cơ cỡ nhỏ bằng AI, ứng dụng nhiên liệu sinh học trên động cơ cải tiến, và tự động hóa quy trình đúc - gia công chính xác nắp xilanh.
  6. Tạo lập di sản công nghệ nội địa hóa vững chắc, đưa sản phẩm động cơ đốt trong mang thương hiệu Việt Nam nâng cao vị thế cạnh tranh trên trường quốc tế.