Giới thiệu dự án

Nghiên cứu cơ học chất lưu và khí động học sinh học (Biomimetic Aerodynamics) đang trở thành xu hướng mũi nhọn trong việc tối ưu hóa hiệu suất năng lượng của các thiết bị cơ nhiệt. Trong tự nhiên, cơ chế tự quay (autorotation) của các loại hạt cây họ Dầu (Dipterocarpaceae) khi rơi tự do giúp chúng giảm đáng kể vận tốc rơi thẳng đứng và phân tán hạt giống đi xa. Theo các thống kê cơ học chất lưu ứng dụng, việc ứng dụng biên dạng khí động học tự nhiên từ quả chò có thể giúp tăng hiệu suất khí động học của cánh quạt tản nhiệt trong kỹ thuật nhiệt lên từ 8% đến 15%, đồng thời giảm độ ồn và tiết kiệm điện năng tiêu thụ.

Vấn đề cốt lõi đặt ra là tại Việt Nam, các nghiên cứu chuyên sâu về mô phỏng số 3D (3D CFD Simulation) quá trình rơi và tương tác dòng khí của quả chò nâu (Dipterocarpus retusus) còn rất hạn chế. Các nghiên cứu trước đây chủ yếu dừng lại ở việc thả rơi thực nghiệm hoặc phân tích video 2D, chưa thể hiện được chi tiết trường vận tốc (velocity streamline), trường áp suất (pressure distribution) và các hệ số khí động học ($C_d, C_l, C_m$) trên bề mặt 3D phức tạp của quả chò.

Đề tài "Nghiên cứu đặc tính khí động học về quá trình rơi của quả chò bằng phần mềm ANSYS FLUENT" do sinh viên Phan Tấn Chí thực hiện dưới sự hướng dẫn của ThS. Nguyễn Lê Hồng Sơn tại Trường Đại học Sư phạm Kỹ thuật TP.HCM (HCMUTE) giải quyết trực tiếp bài toán này với 3 mục tiêu cụ thể:

  1. Xây dựng mô hình 3D chính xác của quả chò nâu (Dipterocarpus retusus) dựa trên các thông số hình học tối ưu từ quy hoạch thực nghiệm (DOE - Design of Experiments).
  2. Mô phỏng động lực học chất lưu số (CFD) dòng không khí chuyển động thẳng đứng tác động vào quả chò trong môi trường phần mềm thương mại chuyên dụng ANSYS FLUENT.
  3. Phân tích định lượng và trực quan hóa phân bố lực nâng ($F_L$), lực cản ($F_D$), mô-men xoay ($M$) và hệ số áp suất ($C_p$), làm cơ sở khoa học để ứng dụng vào thiết kế cánh quạt tản nhiệt trong ngành Công nghệ Kỹ thuật Nhiệt.

Phạm vi nghiên cứu tập trung vào mô hình quả chò mẫu số 13 (mẫu tối ưu từ 81 mẫu thực nghiệm với thời gian rơi $t = 2.64\text{ s}$ và thời gian trung bình kiểm chứng $t_{avg} = 2.79\text{ s}$). Bài toán được thiết lập trong điều kiện dòng khí chảy tầng (laminar flow) với hệ quy chiếu tương đối (vận tốc gió thổi từ dưới lên mô phỏng vận tốc rơi ổn định).

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

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

Trước khi thực hiện mô phỏng số trên ANSYS FLUENT, nhiều phương pháp tiếp cận khác nhau đã được triển khai trên thế giới và tại Việt Nam nhằm khám phá cơ chế bay của các loại quả chò:

Tiêu chí so sánh Phân tích Video Tracker (Dipterocarpus alatus - Thái Lan) Đo sóng siêu âm & Quang học (Triplaris samaras - Ladera et al.) Mô phỏng OpenFOAM (Dryobalanops aromatica - Malaya) Giải pháp ANSYS FLUENT 3D (Đề tài nghiên cứu)
Công cụ chính Phần mềm Tracker 2D Cảm biến siêu âm & Chùm tia quang học SALOME + OpenFOAM + ParaView AutoCAD 2020 + ANSYS Fluent 2020 R2
Bản chất nghiên cứu Thực nghiệm động học 2D Thực nghiệm chuyển động hạt rơi Mô phỏng CFD cánh tuabin thủy triều Mô phỏng CFD 3D trường áp suất & vận tốc
Ưu điểm Dễ triển khai, chi phí thấp Đo chính xác tốc độ góc $\Omega$ và độ cao Mã nguồn mở, ứng dụng vào tuabin biển Phân tích chi tiết lực nâng, cản, mô-men và hội tụ chuẩn xác
Nhược điểm Bỏ qua dòng xoáy 3D quanh cánh Không tính được phân bố áp suất cục bộ Mô hình đơn giản hóa dạng tuabin quay Đòi hỏi tài nguyên tính toán và chất lượng lưới cao 
                       Yêu cầu kỹ thuật hệ thống mô phỏng
   [Must-Have]                   [Should-Have]                  [Could-Have]
- Mô hình CAD 3D chuẩn        - Lưới Tetrahedrons           - Mô phỏng chuyển động
  mặt cắt A-A & B-B             Skewness < 0.85               6 bậc tự do (6-DOF)
- Điều kiện biên              - Báo cáo lực hội tụ          - Khảo sát dòng rối
  Inlet, Outlet, Wall           (Cl, Cd, Cm r-plots)          chuyển tiếp k-omega

Dựa trên phương pháp phân loại yêu cầu MoSCoW:

  • Must-Have (Bắt buộc): Tái tạo chính xác hình học hạt và 2 cánh quả chò trên AutoCAD; chia lưới thể tích điều khiển chất lưu với độ méo (Skewness) nằm trong ngưỡng tin cậy; thiết lập đầy đủ điều kiện biên (Velocity Inlet, Pressure Outlet, No-slip Wall).
  • Should-Have (Nên có): Báo cáo hội tụ đồ thị lực kháng ($C_d$), lực nâng ($C_l$), và mô-men ($C_m$) theo từng bước lặp tính toán (iteration).
  • Could-Have (Có thể mở rộng): Mở rộng sang mô hình chuyển động 6 bậc tự do (6-DOF) và mô hình hóa dòng rối chuyển tiếp.

Thiết kế hệ thống mô phỏng

Quy trình kỹ thuật được thiết kế theo kiến trúc module tích hợp thông qua nền tảng ANSYS Workbench:

  1. CAD Engine (AutoCAD 2020): Thiết lập biên dạng cánh 3D, mặt cắt ngang A-A, mặt cắt dọc B-B, góc mở cánh và phần thân hạt hình cầu elip.
  2. Pre-processing (DesignModeler): Nhập hình học (Import Geometry), tạo khối bao lưu chất không khí (Air Enclosure) và thực hiện phép toán Boolean Subtraction để trừ thể tích vật thể quả chò khỏi miền tính toán.
  3. Meshing Module (ANSYS Meshing): Rời rạc hóa miền tính toán bằng phần tử tứ diện (Tetrahedrons), kiểm soát hệ số Skewness và Aspect Ratio.
  4. Solver Core (ANSYS Fluent 2020 R2): Bộ giải dựa trên phương pháp Thể tích hữu hạn (Finite Volume Method - FVM) giải hệ phương trình bảo toàn Navier-Stokes ở trạng thái dừng (steady-state).

Cơ sở phương pháp luận (Methodology)

Phương pháp số sử dụng hệ phương trình bảo toàn tổng quát (Transport Equation) tích phân trên một thể tích điều khiển $V$:

$$\frac{\partial}{\partial t} \int_V \rho \phi , dV + \oint_A \rho \phi \vec{v} \cdot d\vec{A} = \oint_A \Gamma_\phi \nabla \phi \cdot d\vec{A} + \int_V S_\phi , dV$$

Trong đó:

  • $\phi$ đại diện cho các biến trạng thái (1 đối với khối lượng, $u, v, w$ đối với động lượng).
  • $\Gamma_\phi$ là hệ số khuếch tán.
  • $S_\phi$ là số hạng nguồn sinh (Generation term).

Trong giai đoạn quả chò tự quay ổn định, lực nâng ($F_L$), lực cản ($F_D$) và mô-men quay ($M_z$) được xác định bằng các biểu thức tích phân áp suất và ứng suất cắt trên bề mặt cánh:

$$F_D = \frac{1}{2} \rho v^2 A C_D(\alpha), \quad F_L = \frac{1}{2} \rho v^2 A C_L(\alpha)$$

$$M_z = \frac{1}{4} l^2 \rho v^2 C_L - \frac{1}{2} l^2 \rho v^2 C_D$$

Trong đó $\rho$ là mật độ không khí ($1.225\text{ kg/m}^3$), $A$ là diện tích bề mặt cánh quả chò, $l$ là chiều dài cánh, $v$ là vận tốc dòng khí tương đối và $\alpha$ là góc tấn (angle of attack).

Implementation và kết quả

Quy trình triển khai thực nghiệm và số trị

Quá trình mô phỏng được thực thi qua 4 giai đoạn logic chặt chẽ:

[Giai đoạn 1: Đo lường & Xử lý số liệu DOE]
[Giai đoạn 2: Tạo hình CAD 3D]
[Giai đoạn 3: Phân chia lưới & Thiết lập điều kiện biên]
[Giai đoạn 4: Tính toán lặp FVM & Hậu xử lý]
      - Solver: Pressure-Based, Laminar Model
      - Hội tụ tại vòng lặp thứ 1686
      - Trích xuất: Streamlines, Contours, Cl, Cd, Cm

Các thông số điều kiện biên và vật liệu được cấu hình trực tiếp trong ANSYS Fluent:

/* Pseudocode cấu hình điều kiện biên và tính chất vật lý trong ANSYS Fluent */
Domain_Setup() {
    Fluid_Material("air") {
        Density = 1.225;             // kg/m3 (Mật độ không khí tiêu chuẩn)
        Viscosity = 1.7894e-05;      // kg/(m.s) (Độ nhớt động học)
    }
    Boundary_Conditions() {
        Zone("Inlet")  = VELOCITY_INLET(Velocity = 2.5 m/s, Direction = Normal_to_Boundary);
        Zone("Outlet") = PRESSURE_OUTLET(Gauge_Pressure = 0 Pa);
        Zone("Samara_Wall") = WALL(Motion = Stationary, Shear = No_Slip);
    }
    Solver_Controls() {
        Scheme = SIMPLE;
        Spatial_Discretization_Momentum = Second_Order_Upwind;
        Convergence_Criteria = 1e-04;
    }
}

Đánh giá chất lượng lưới và kiểm tra hội tụ

Chất lượng hình học của miền lưới là yếu tố quyết định tính hội tụ của bài toán CFD. Lưới được chia theo phương pháp Tetrahedrons (tứ diện) cho phép thích ứng tối đa với các bề mặt cong phức tạp ở mép cánh quả chò.

Thông số kiểm tra lưới Giá trị đạt được Ngưỡng tiêu chuẩn ANSYS Đánh giá chất lượng
Max Skewness (Độ méo tối đa) 0.78 $< 0.85$ (Good) Đạt chuẩn hội tụ
Average Skewness 0.23 $< 0.25$ (Excellent) Xuất sắc
Aspect Ratio (Tỷ số hình học) 3.45 $< 10.0$ Đảm bảo tính đẳng hướng
Tổng số phần tử (Elements) ~450,000 Phù hợp tài nguyên tính toán Đảm bảo tính độc lập lưới 
Tiêu chí Skewness:
  [Excellent]      [Good]          [Fair]        [Poor]     [Bad]

Quá trình lặp tính toán giải phương trình Navier-Stokes đạt tiêu chuẩn hội tụ tại vòng lặp thứ 1686 (Iteration 1686), khi các giá trị thặng dư (Residuals) giảm xuống dưới $10^{-4}$ và các đường giám sát báo cáo lực (drag-rplot, lift-rplot, moment-rplot) đạt trạng thái tiệm cận nằm ngang ổn định.

Residual Log Scale
      0       400       800      1200      1686                 Iterations

Kết quả phân tích khí động học chi tiết

Kết quả trích xuất từ mô-đun Post-Processing thể hiện rõ bản chất khí động học của quả chò:

                            Trường vận tốc & Đường dòng
                           (Velocity Streamlines trên XY)
                                  ~~~~|~~~~
                                ~~~ /   \ ~~~
                      Vùng xoáy -> ( ( O ) ) <- Phân tách dòng
                                 ~~~ \ / ~~~
                                  ====|==== (Cánh quả chò)
                                (Inlet: 2.5 m/s)
  1. Trường vận tốc (Velocity Streamline & Vectors): Dòng khí tiếp cận mặt dưới của cánh bị cản lại và chuyển hướng dọc theo độ nghiêng hình học. Vận tốc dòng tăng đột ngột tại mép trước (leading edge) của cánh, đạt giá trị cực đại khi vượt qua đỉnh cong và tạo ra một vùng xoáy cuộn phân tách dòng ở mặt trên của cánh.
  2. Trường áp suất (Pressure Distribution): Phân bố áp suất thể hiện sự chênh lệch rõ rệt giữa mặt đón gió (vùng áp suất dương cao, $P > 0$) và mặt lưng cánh (vùng áp suất âm thấp, hình thành do hiện tượng hút chân không cục bộ). Sự chênh lệch áp suất này theo định luật Bernoulli sinh ra lực nâng khí động học ($F_L$).
  3. Mô-men xoay ($M_z$): Góc nghiêng của hai cánh so với trục thẳng đứng tạo ra lực khí động tiếp tuyến không đối xứng, tạo nên mô-men xoay tự do quanh trục thẳng đứng, duy trì chuyển động quay tít và hãm vận tốc rơi tịnh tiến.

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

Nghiên cứu mang lại những đóng góp kỹ thuật và cải tiến mang tính thực tiễn cao:

                 So sánh hiệu quả phân tích giữa các phương pháp
  • Số hóa biên dạng sinh học phức tạp: Chuyển đổi thành công cấu trúc cánh xoắn tự nhiên của quả chò Dipterocarpus retusus sang mô hình CAD 3D chính xác với các mặt cắt ngang A-A và dọc B-B, phục vụ trực tiếp cho bài toán mô phỏng số.
  • Tối ưu hóa hình học dựa trên dữ liệu thực nghiệm: Kết hợp chặt chẽ giữa ma trận quy hoạch thực nghiệm DOE (Minitab) từ 81 mẫu thực tế với mô phỏng số 3D, chứng minh rằng Góc cánh (Blade Angle) là tham số có trọng số ảnh hưởng lớn nhất đến thời gian rơi và độ ổn định quay.
  • Tiết kiệm chi phí thử nghiệm khí động: Giảm hơn 70% chi phí chế tạo mô hình thử nghiệm ống khí động (Wind Tunnel Testing) nhờ ứng dụng công nghệ mô phỏng số CFD trên ANSYS Fluent.

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

Kết quả nghiên cứu mở ra nhiều định hướng ứng dụng kỹ thuật giá trị:

                       CÁC HƯỚNG ỨNG DỤNG THỰC TẾ
  [Kỹ thuật Nhiệt]        [Thiết bị Bay Mini]      [Năng lượng Tái tạo]
- Cánh quạt tản nhiệt     - Drone đơn cánh tự quay - Tuabin gió trục đứng
  cho hệ thống HVAC & CPU   (Monocopter MAV)         hiệu suất cao
- Tăng lưu lượng khí 12%  - Thiết bị cứu hộ tự rơi - Khởi động ở vận tốc gió
- Giảm tiếng ồn cánh        không cần động cơ        thấp (< 2 m/s)
  1. Ngành Công nghệ Kỹ thuật Nhiệt (HVAC & Cooling Systems):
    • Ứng dụng góc xoắn và độ cong cánh quả chò vào thiết kế cánh quạt dàn ngưng, quạt tản nhiệt máy tính và tháp giải nhiệt.
    • Giúp phân bổ đều luồng khí qua két tản nhiệt, giảm hiện tượng tách dòng ở đầu cánh, nâng cao hiệu suất truyền nhiệt đối lưu cưỡng bức thêm 8-12%.
  2. Hàng không và Thiết bị bay không người lái (UAV/MAV):
    • Thiết kế các cảm biến khí tượng hoặc thiết bị quan trắc môi trường dạng hạt tự quay (biomimetic samara drones) có khả năng tự thả rơi từ độ cao lớn để thu thập dữ liệu trong thời gian dài mà không tiêu hao năng lượng.
  3. Tuabin gió mô phỏng sinh học (Biomimetic Wind Turbines):
    • Thiết kế cánh tuabin gió trục đứng tự khởi động ở dải vận tốc gió thấp ($v < 2\text{ m/s}$), ứng dụng cho các trạm phát điện nhỏ ở đô thị hoặc hải đảo.

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

Hạn chế kỹ thuật

  • Giả định dòng chảy dừng (Steady-State): Bài toán hiện tại xem xét dòng khí chuyển động tương đối qua mô hình quả chò cố định, chưa mô phỏng trực tiếp tương tác đa vật thể động (Dynamic Mesh / 6-DOF) theo thời gian thực (Transient).
  • Mô hình dòng chảy tầng: Giả định dòng laminar ở vận tốc thấp ($v = 2.5\text{ m/s}$) chưa bao quát hết hiện tượng chuyển tiếp dòng rối (transitional turbulence) khi quả chò quay ở tốc độ góc cao trong điều kiện gió giật.

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

  • Tích hợp mô hình chuyển động 6 bậc tự do (6-DOF) & Lưới trượt (Overset Mesh): Mô phỏng trực tiếp quá trình lật ngược ban đầu ($t_0 \to t^*$) và chuyển sang chế độ tự quay ổn định trong không gian 3 chiều.
  • Tương tác kết cấu - lưu chất (FSI - Fluid-Structure Interaction): Nghiên cứu độ biến dạng đàn hồi của cánh quả chò dưới tác dụng của áp lực khí động học để tối ưu hóa vật liệu composite cho cánh quạt công nghiệp.

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

                                  ĐỐI TƯỢNG HƯỞNG LỢI
    [Sinh viên]         [Kỹ sư R&D]       [Doanh nghiệp]     [Nhà nghiên cứu]
- Nắm vững quy trình   - Áp dụng biên     - Rút ngắn chu kỳ   - Bổ sung dữ liệu
  CFD từ CAD đến       dạng tự nhiên      R&D cánh quạt       khí động học sinh
  hội tụ Solver        vào sản phẩm       công nghiệp 40%     học thực nghiệm
  • Sinh viên & Học viên ngành Cơ khí/Nhiệt: Cung cấp tài liệu mẫu hoàn chỉnh về phương pháp luận kết hợp giữa nghiên cứu thực nghiệm (DOE Minitab) và mô phỏng số (ANSYS Fluent CFD).
  • Kỹ sư thiết kế thiết bị nhiệt & cơ khí: Tham khảo quy trình phân tích lực khí động ($C_d, C_l, C_m$) và trường áp suất để tối ưu hóa cánh quạt tản nhiệt, van điều tiết khí.
  • Doanh nghiệp sản xuất quạt công nghiệp: Giảm 40% chi phí thử nghiệm thực tế (prototyping cost) nhờ khả năng dự đoán chính xác qua mô phỏng số.
  • Cộng đồng nghiên cứu Khí động học sinh học: Đóng góp bộ cơ sở dữ liệu số trị về quả chò Dipterocarpus retusus tại khu vực Đông Nam Á vào kho tàng nghiên cứu khí động học thế giới.

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

1. Cấu hình máy tính tối thiểu để chạy mô phỏng ANSYS Fluent cho bài toán này là gì?

Hệ thống cần tối thiểu CPU 4 nhân (khuyến nghị 8 nhân Intel Core i7 hoặc AMD Ryzen 7), 16GB RAM (khuyến nghị 32GB RAM để xử lý lưới trên 1 triệu phần tử) và card đồ họa chuyên dụng hỗ trợ OpenGL để hiển thị mượt mà trường vector và streamline trong phần hậu xử lý.

2. Tại sao góc cánh lại là yếu tố quan trọng nhất ảnh hưởng đến thời gian rơi của quả chò?

Theo biểu đồ Pareto từ quy hoạch thực nghiệm DOE, góc cánh trực tiếp quyết định góc tấn ($\alpha$) của biên dạng cánh đối với luồng khí. Góc tấn này chi phối tỷ số lực nâng trên lực cản ($C_L/C_D$) và mô-men quay tự do ($M_z$), từ đó kiểm soát trực tiếp vận tốc rơi ổn định cuối cùng (terminal velocity).

3. Làm thế nào để kiểm tra tính độc lập của lưới (Grid Independence) trong bài toán?

Cần tiến hành chia lưới với ít nhất 3 mật độ phần tử khác nhau (ví dụ: 250.000 phần tử, 450.000 phần tử và 800.000 phần tử). Nếu sai số của lực nâng $F_L$ và lực cản $F_D$ giữa các cấp lưới nhỏ hơn 2%, bài toán đạt tính độc lập của lưới và có thể chọn cấp lưới 450.000 phần tử để tối ưu thời gian tính toán.

4. Kết quả mô phỏng CFD có thể xuất sang các định dạng nào để gia công hoặc tối ưu tiếp?

Biên dạng cánh sau khi tối ưu có thể xuất trực tiếp từ DesignModeler/SpaceClaim sang các định dạng chuẩn CAD/CAM như STEP, IGES, Parasolid (.x_t) để nạp vào máy in 3D tạo mẫu nhanh hoặc máy phay CNC 5 trục gia công khuôn đúc cánh quạt.

5. Tiêu chuẩn hội tụ nào khẳng định kết quả mô phỏng là đáng tin cậy?

Bài toán được xác nhận hội tụ khi đồng thời thỏa mãn 3 điều kiện:

  • Sai số thặng dư (Residuals) của phương trình liên tục và động lượng giảm dưới $10^{-4}$.
  • Đồ thị theo dõi lực cản (drag-rplot) và lực nâng (lift-rplot) ổn định hoàn toàn không đổi qua ít nhất 100 vòng lặp liên tiếp (đạt trạng thái tiệm cận tại vòng lặp 1686).
  • Cân bằng lưu lượng tổng khối lượng giữa ngõ vào (Inlet) và ngõ ra (Outlet) có sai số $< 0.1%$.

Kết luận

Đồ án tốt nghiệp "Nghiên cứu đặc tính khí động học về quá trình rơi của quả chò bằng phần mềm ANSYS FLUENT" của sinh viên Phan Tấn Chí đã hoàn thành xuất sắc các nhiệm vụ khoa học đặt ra. Bằng việc kết hợp hài hòa giữa cơ sở dữ liệu thực nghiệm DOE và công cụ mô phỏng số tiên tiến ANSYS Fluent, nghiên cứu đã làm sáng tỏ cơ chế khí động học tự quay của quả chò nâu (Dipterocarpus retusus).

Kết quả mô phỏng không chỉ giải thích chi tiết trường phân bố áp suất và vận tốc giúp quả chò duy trì thời gian bay tối ưu trong tự nhiên, mà còn mở ra tiềm năng ứng dụng thực tiễn to lớn trong việc cải tiến thiết kế cánh quạt tản nhiệt, nâng cao hiệu quả năng lượng cho ngành Kỹ thuật Nhiệt và hệ thống thiết bị cơ điện hiện đại.