Nghiên Cứu Ảnh Hưởng Của Xoáy Lốc Trên Đường Ống Nạp Đến Hiệu Suất Động Cơ Xe Máy: Tối Ưu Hóa Khí Động Học Bằng Mô Phỏng CFD và MATLAB/Simulink

Tóm tắt nghiên cứu (Research Summary)

  • Câu hỏi nghiên cứu chính (Research Question): Sự biến thiên hình học của đường ống nạp—cụ thể là góc nghiêng cổ nạp—ảnh hưởng như thế nào đến các hệ số xoáy lốc dọc (Swirl ratio), xoáy lốc ngang (Tumble ratio), và qua đó tác động trực tiếp đến hiệu suất nhiệt, công suất đầu ra, mô-men xoắn cùng suất tiêu hao nhiên liệu của động cơ xe máy 125cc?
  • Phương pháp tiếp cận (Methodology Snapshot): Đề tài tích hợp mô hình hóa hình học 3D bằng CATIA V5, mô phỏng động lực học lưu chất tính toán (CFD - ANSYS Fluent ICE) với mô hình rối $k-\varepsilon$ trên động cơ xe máy Honda Future FI 125cc, kết hợp mô phỏng nhiệt động lực học và đặc tính cháy bằng MATLAB/Simulink. Các kịch bản góc nghiêng cổ nạp được khảo sát ở các dải vận tốc góc trục khuỷu khác nhau (1800, 3000, 4200 vòng/phút).
  • Phát hiện chính (Key Findings): Cấu hình góc nghiêng cổ nạp cải tiến $30^\circ$ mang lại tỷ lệ xoáy lốc tổng hợp tối ưu nhất so với mẫu nguyên bản ($25^\circ$) và biến thể ($35^\circ$). Cấu trúc xoáy lốc tối ưu này thúc đẩy mạnh mẽ động năng dòng nạp, cải thiện độ đồng nhất của hòa khí, nâng cao công suất phát và mô-men xoắn, đồng thời hạ thấp rõ rệt suất tiêu hao nhiên liệu riêng.
  • Hàm ý thực tiễn (Implications): Kết quả khẳng định việc tinh chỉnh thông số hình học đường nạp là giải pháp can thiệp thụ động chi phí thấp nhưng mang lại hiệu quả cao, giúp các nhà sản xuất xe máy và kỹ sư thiết kế tối ưu hóa động cơ đốt trong mà không cần thay đổi phức tạp toàn bộ cụm nắp máy.

Bối cảnh và tầm quan trọng (Context & Significance)

Thực trạng kỹ thuật và môi trường

Tại Việt Nam và nhiều quốc gia đang phát triển ở châu Á, xe máy là phương tiện giao thông cá nhân chủ đạo nhờ tính cơ động, kinh tế và linh hoạt. Tuy nhiên, mật độ phương tiện quá lớn đi kèm với sự phát thải khí nhà kính và các chất ô nhiễm độc hại (HC, CO, $\text{NO}_x$, bụi mịn PM) đang gây sức ép nặng nề lên chất lượng không khí đô thị. Phần lớn động cơ xe máy thương mại dung tích nhỏ hiện nay sử dụng hệ thống nạp có kết cấu đơn giản, chưa được tối ưu hóa về mặt khí động học buồng cháy, dẫn đến quá trình hòa trộn nhiên liệu - không khí kém đồng đều, đặc biệt tại các dải tải cục bộ và tốc độ vòng tua thấp.

flowchart LR
    A["Kết cấu ống nạp đơn giản"] --> B["Dòng khí nạp kém xoáy lốc"]
    B --> C["Hòa khí không đồng nhất"]
    C --> D["Cháy không hoàn toàn"]
    D --> E["Tốn nhiên liệu & Tăng phát thải khí xả"]
    
    style A fill:#fee,stroke:#d33
    style E fill:#fee,stroke:#d33
    style C fill:#ffe,stroke:#fa0

Khoảng trống nghiên cứu (Research Gap)

Dù các giải pháp van biến thiên xoáy lốc như T-VIS của Toyota, SCV của Nissan hay STCS của SYM đã chứng minh tính hiệu quả trên ô tô và một số dòng xe tay ga cao cấp, việc áp dụng các cơ cấu chấp hành điện tử - chân không phức tạp này lên các dòng xe máy phổ thông phân khúc 125cc gặp trở ngại lớn về chi phí chế tạo, không gian lắp đặt và bảo trì. Nhiều công trình nghiên cứu trước đây tại Việt Nam chủ yếu tiếp cận theo hướng đo đạc thực nghiệm tốn kém hoặc chỉ phân tích tĩnh mà thiếu vắng chuỗi mô phỏng liên kết đa miền giữa dòng chảy động học 3D (CFD) và chu trình nhiệt động toàn diện của động cơ.

Tính cấp thiết và tác động kỳ vọng

Nghiên cứu này giải quyết bài toán giao thoa giữa:

  1. Kỹ thuật: Tối đa hóa công suất và mô-men xoắn tại mọi dải tốc độ vận hành.
  2. Kinh tế: Giảm thiểu suất tiêu hao nhiên liệu định mức cho người sử dụng.
  3. Môi trường: Giảm phát thải khí xả thông qua việc hoàn thiện quá trình cháy sạch từ gốc.

Phương pháp luận và cách tiếp cận (Methodology & Approach)

Nghiên cứu áp dụng quy trình tính toán số tích hợp liên thông, chuyển tiếp từ đo đạc kích thước thực tế, mô hình hóa CAD 3D, mô phỏng khí động học CFD đến phân tích đặc tính chu trình động cơ bằng Simulink.

flowchart TD
    M["Thu thập thông số động cơ thực tế\n(Honda Future FI 125cc)"] --> CAD["Mô hình hóa hình học 3D\n(CATIA V5 - Part Design, Rib, Pocket)"]
    CAD --> MESH["Rời rạc hóa không gian & Tạo lưới\n(ANSYS Workbench / Mesh Generator)"]
    MESH --> CFD["Mô phỏng động lực học chất lưu CFD\n(ANSYS Fluent ICE - k-epsilon Turbulence)"]
    CFD --> EXTRACT["Trích xuất hệ số Swirl, Tumble\n& Năng lượng động năng dòng chảy"]
    EXTRACT --> SIM["Mô phỏng chu trình cháy & Đặc tính động cơ\n(MATLAB/Simulink)"]
    SIM --> EVAL["Đánh giá Công suất, Mô-men xoắn,\nSuất tiêu hao nhiên liệu"]

1. Thông số đối tượng nghiên cứu

Nghiên cứu thực hiện trên động cơ xe máy Honda Future FI 125cc với các thông số kết cấu định mức:

  • Đường kính xi lanh ($B$): $52.4\text{ mm}$; Bán kính trục khuỷu ($R$): $28.95\text{ mm}$; Chiều dài thanh truyền: $101.5\text{ mm}$.
  • Dung tích buồng cháy: $15\text{ cm}^3$; Tỉ số nén ($\varepsilon$): $9.3:1$; Số lượng van: 2 xú-pap (1 nạp, 1 xả).
  • Đường kính danh định cổ nạp: $23.2\text{ mm}$.
  • Góc pha phân phối khí: Góc mở sớm van nạp (IVO) $82^\circ\text{ CA}$, góc đóng muộn van nạp (IVC) $212^\circ\text{ CA}$; Góc mở sớm van xả (EVO) $521^\circ\text{ CA}$, góc đóng muộn van xả (EVC) $640^\circ\text{ CA}$.
  • Hệ thống nhiên liệu: Phun xăng điện tử đa điểm, áp suất dàn phun $10\text{ MPa}$, góc đánh lửa sớm $345^\circ\text{ CA}$.

2. Kỹ thuật chia lưới và giải thuật CFD

  • Mô hình hình học nắp máy, buồng đốt, van nạp - xả và cổ nạp được thiết kế chính xác trên CATIA V5 bằng các lệnh nâng cao như Profile, Rib, Shaft, Pocket.
  • Định dạng trung gian .isg được chuyển vào mô-đun ANSYS ICE (Internal Combustion Engine).
  • Quá trình chia lưới động (Dynamic Meshing) thích ứng theo hành trình chuyển động của piston và độ nâng hạ van xú-pap.
  • Giải thuật mô phỏng sử dụng mô hình dòng chảy rối $k-\varepsilon$ ($k$-epsilon turbulence model) giải phương trình Navier-Stokes cho dòng nén được, tái hiện chính xác phổ vận tốc, phân bố áp suất và xoáy lốc xuyên suốt toàn chu trình nạp - nén.

3. Mô hình toán nhiệt động học trên MATLAB/Simulink

Dữ liệu động năng dòng chảy và các hệ số xoáy từ ANSYS Fluent được trích xuất để làm biến đầu vào cho mô hình toán học trên MATLAB/Simulink, tính toán các chỉ tiêu cơ nhiệt:

  • Vận tốc trung bình và mô-men động lượng dòng xoáy: $$\vec{\Omega}_A = \vec{\Omega}_S + \vec{\Omega}_T$$
  • Hệ số xoáy lốc dọc ($C_S$) và xoáy lốc ngang ($C_T$): $$C_S = \frac{8 \cdot \Omega_S}{m \cdot B \cdot v_0}, \quad C_T = \frac{8 \cdot \Omega_T}{m \cdot B \cdot v_0}$$
  • Tích hợp hàm toả nhiệt và phương trình cân bằng năng lượng để xuất ra đường đặc tính ngoài của động cơ.

Phát hiện chính của nghiên cứu (Key Findings)

Qua việc khảo sát 3 biến thể góc nghiêng đường ống nạp ($25^\circ$ - mẫu nguyên bản, $30^\circ$ - cải tiến 1, $35^\circ$ - cải tiến 2) tại các dải tốc độ 1800, 3000 và 4200 vòng/phút, nghiên cứu rút ra các kết quả định lượng nổi bật:

Thông số đánh giá Cổ nạp nguyên bản ($25^\circ$) Cổ nạp cải tiến ($30^\circ$) Cổ nạp thực nghiệm ($35^\circ$) Nhận xét xu hướng
Cường độ xoáy Swirl ($C_S$) Trung bình thấp Tối ưu cao nhất Giảm do va đập thành ống $30^\circ$ tạo dòng cuộn dọc đều quanh thân nấm xú-pap
Cường độ xoáy Tumble ($C_T$) Trung bình Tăng trưởng mạnh Phân tán rối cục bộ Dòng khí hướng vào lòng xilanh tạo xoáy lốc đứng
Năng lượng động năng dòng Bị tổn thất do cản áp Duy trì ổn định Suy giảm ở cuối kỳ nạp Cải thiện hiệu suất nạp $\eta_v$
Mô-men xoắn cực đại Giá trị cơ sở Tăng từ 4.2% - 6.8% Thấp hơn mức $30^\circ$ Nâng cao lực kéo ở dải tua thấp và trung bình
Công suất động cơ ($N_e$) Giá trị cơ sở Cải thiện đáng kể Tương đương nguyên bản Tốc độ cháy tăng, rút ngắn góc cháy trễ
Suất tiêu hao nhiên liệu ($g_e$) Cao nhất Thấp nhất toàn dải Tăng nhẹ Hòa khí đồng nhất giúp cháy kiệt nhiên liệu 
flowchart TD
    subgraph GOC25 ["Góc nghiêng 25° (Nguyên bản)"]
        A1["Xoáy Swirl/Tumble thấp"] --> B1["Hòa trộn hòa khí chậm"] --> C1["Hiệu suất nhiệt trung bình"]
    end

    subgraph GOC30 ["Góc nghiêng 30° (Tối ưu)"]
        A2["Tỷ lệ Swirl & Tumble cực đại hài hòa"] --> B2["Tốc độ lan tràn màng lửa cực đại"] --> C2["Công suất & Moment cao nhất, BSFC thấp nhất"]
    end

    subgraph GOC35 ["Góc nghiêng 35° (Quá dốc)"]
        A3["Dòng nạp cản trở, ma sát tăng"] --> B3["Năng lượng động năng dòng phân tán"] --> C3["Hiệu suất giảm nhẹ"]
    end

    style GOC30 fill:#dfd,stroke:#2b2,stroke-width:2px
    style GOC25 fill:#f9f9f9,stroke:#999
    style GOC35 fill:#fee,stroke:#f88

Phân tích chi tiết động học dòng chảy

  1. Góc nghiêng $30^\circ$ là điểm tới hạn khí động học: Khi tăng góc nghiêng từ $25^\circ$ lên $30^\circ$, độ uốn lượn của dòng khí khi đi qua xú-pap nạp giảm bớt góc chết, hướng dòng tiếp tuyến áp sát thành xilanh hiệu quả hơn, gia tăng mô-men động lượng ban đầu $J_i$.
  2. Khắc phục hiện tượng suy thoái xoáy ở $35^\circ$: Khi tăng độ nghiêng lên $35^\circ$, tiết diện lưu thông hiệu dụng tại vùng chuyển tiếp giữa cổ nạp và đế van bị co thắt, gây tổn thất áp suất cục bộ và sinh ra các vùng xoáy quẩn tự do làm triệt tiêu lẫn nhau, làm giảm vận tốc tiếp tuyến vào buồng đốt.
  3. Ảnh hưởng trực tiếp đến chu trình cháy: Việc tăng cường tỷ lệ $C_S$ và $C_T$ ở cấu hình $30^\circ$ giúp các hạt xăng phun sương (đường kính ti phun $0.6\text{ mm}$) bay hơi nhanh hơn, giảm màng bám ướt trên thành xilanh, tạo điều kiện thuận lợi cho bugi kích hoạt tâm cháy lan tỏa nhanh khắp thể tích công tác.

Đóng góp khoa học và giá trị học thuật (Scientific Contributions)

Về mặt lý thuyết (Theoretical Contributions)

  • Hoàn thiện mô hình toán học giải thích sự tương tác phức hợp giữa cấu trúc xoáy dọc (Swirl) và xoáy ngang (Tumble) đối với quá trình hình thành hòa khí trong động cơ xe máy cỡ nhỏ 4 kỳ 2 van.
  • Cung cấp dữ liệu định lượng chuẩn xác về hệ số tổn thất mô-men động lượng do ma sát màng khí trong quá trình hút và nén dưới dạng phương trình giải tích tường minh.

Về mặt phương pháp luận (Methodological Innovations)

  • Xây dựng quy trình liên kết chuẩn giữa CATIA V5 $\rightarrow$ ANSYS Fluent ICE $\rightarrow$ MATLAB/Simulink, rút ngắn thời gian thử nghiệm ảo và cho phép tinh chỉnh nhanh các thông số kết cấu buồng cháy.
  • Giải quyết bài toán chia lưới động thích ứng cho động cơ đốt trong với chi phí tài nguyên tính toán hợp lý nhưng vẫn bảo đảm độ tin cậy của trường vận tốc dòng chảy.

Về mặt thực tiễn (Practical Applications)

  • Cung cấp giải pháp chế tạo/đúc lại cổ nạp cải tiến nghiêng $30^\circ$ cho dòng xe Honda Future FI 125cc và các dòng xe tương đương trên thị trường, nâng cao tính kinh tế nhiên liệu trực tiếp cho phương tiện giao thông phổ thông.
  • Làm tài liệu tham khảo kỹ thuật nền tảng cho công tác giảng dạy, nghiên cứu chuyên ngành Kỹ thuật Cơ khí Động lực tại các trường đại học kỹ thuật.

Đối tượng hưởng lợi và phạm vi ứng dụng (Target Audience & Benefits)

  • Kỹ sư R&D và Nhà sản xuất ô tô - xe máy: Nắm giữ quy trình mô phỏng CFD tối ưu để ứng dụng trực tiếp vào việc thiết kế nắp máy, ống góp nạp cho các thế hệ động cơ mới nhằm thỏa mãn các tiêu chuẩn khí thải khắt khe (Euro 3, Euro 4, Euro 5).
  • Các nhà nghiên cứu và Học viên sau đại học: Tiếp cận bộ cơ sở dữ liệu mô hình toán động cơ và phương pháp thiết lập bài toán trong ANSYS Fluent ICE để mở rộng cho các nghiên cứu về nhiên liệu sinh học (E5, E10, E85) hoặc động cơ lai (Hybrid).
  • Các xưởng kỹ thuật và Đơn vị nâng cấp hiệu năng động cơ: Ứng dụng nguyên lý thay đổi góc nghiêng cổ nạp để nâng cao hiệu suất nạp tự nhiên mà không cần can thiệp sâu vào phần mềm ECU hay độ chế tỷ số nén nguy hiểm.

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

1. Phát hiện quan trọng nhất của đề tài nghiên cứu này là gì?

Kết quả nghiên cứu xác định góc nghiêng cổ nạp $30^\circ$ là cấu hình tối ưu nhất cho động cơ xe máy 125cc. Tại góc nghiêng này, hệ số xoáy lốc dọc (Swirl) và xoáy lốc ngang (Tumble) đạt giá trị cực đại hài hòa, giúp tăng công suất và mô-men xoắn động cơ trong khi giảm thiểu suất tiêu hao nhiên liệu so với thiết kế nguyên bản $25^\circ$.

2. Sự kết hợp giữa CATIA, ANSYS Fluent và MATLAB/Simulink có ưu điểm gì nổi bật?

Phương pháp tiếp cận chuỗi khép kín này giúp liên kết trực tiếp giữa mô hình hóa hình học 3D chi tiết (CATIA), tính toán trường vận tốc động lực học chất lưu đa chiều (ANSYS CFD) và mô phỏng phản ứng nhiệt động học - đặc tính công suất toàn chu trình (Simulink). Quy trình này giúp tiết kiệm tối đa chi phí chế tạo mẫu thử nghiệm thực tế.

3. Kết quả tối ưu góc nghiêng $30^\circ$ có thể áp dụng cho các dòng xe máy khác không?

Nguyên lý tăng cường xoáy lốc thông qua hình học cổ nạp hoàn toàn có thể khái quát hóa. Tuy nhiên, giá trị góc nghiêng tối ưu chính xác ($30^\circ$) phụ thuộc vào tỉ số $S/D$ (hành trình/đường kính xilanh), bố trí góc xú-pap và thể tích buồng cháy của từng dòng động cơ cụ thể, do đó cần chạy lại mô phỏng tinh chỉnh theo từng kích thước buồng đốt.

4. Hướng phát triển tiếp theo của nghiên cứu này là gì?

Nghiên cứu có thể mở rộng khảo sát trên các dòng động cơ 4 van (DOHC/SOHC 4V), tích hợp mô hình phun nhiên liệu trực tiếp (GDI), đánh giá phát thải chi tiết từng thành phần ($\text{NO}_x$, CO, unburned HC), và tiến hành đo đạc trên băng thử công suất động cơ (Dynamometer) thực tế.

5. Việc thay đổi góc nghiêng cổ nạp có gây khó khăn cho việc lắp ráp thực tế không?

Góc nghiêng từ $25^\circ$ lên $30^\circ$ là mức biến thiên nhỏ về mặt không gian hình học. Thiết kế này hoàn toàn có thể gia công đúc hoặc phay CNC trực tiếp trên thân cổ hút mà không làm thay đổi vị trí gá đặt bình xăng con/họng ga hay chạm vào khung sườn xe máy hiện hữu.

Kết luận (Conclusion)

Công trình nghiên cứu khoa học "Nghiên cứu ảnh hưởng xoáy lốc trên đường ống nạp đến hiệu suất động cơ xe máy" (Mã số: T2018-21TĐ) do PGS. Lý Vĩnh Đạt chủ nhiệm đã chứng minh tính ưu việt của việc tối ưu hóa khí động học đường nạp thông qua công cụ mô phỏng số chuyên sâu. Bằng việc xác lập góc nghiêng cổ nạp tối ưu $30^\circ$, động cơ đạt được sự cải thiện vượt bậc về khả năng hòa trộn khí nạp, gia tăng công suất phát và tiết kiệm nhiên liệu rõ rệt. Đây là tiền đề khoa học vững chắc mở ra hướng thiết kế nắp máy hiệu năng cao, chi phí thấp cho ngành công nghiệp sản xuất phương tiện hai bánh tại Việt Nam.