Giới thiệu dự án
Sự bùng nổ của vi điện tử và bán dẫn với mật độ tích hợp transistor siêu cao theo định luật Moore đang đẩy mật độ dòng nhiệt (heat flux) tại các bộ xử lý (CPU, GPU) và linh kiện công suất vượt ngưỡng $100\text{ W/cm}^2$. Theo các báo cáo kỹ thuật ngành nhiệt điện tử, việc vượt quá nhiệt độ tới hạn $85^\circ\text{C}$ làm giảm $50%$ tuổi thọ linh kiện bán dẫn cho mỗi mức tăng thêm $10^\circ\text{C}$. Các giải pháp tản nhiệt đối lưu tự nhiên hoặc quạt tản nhiệt khí truyền thống đã hoàn toàn chạm ngưỡng giới hạn nhiệt động lực học khi không thể giải quyết triệt để các điểm nóng cục bộ (thermal hotspots) trong không gian siêu mỏng.
Đồ án tốt nghiệp chuyên ngành Công nghệ Kỹ thuật Nhiệt (Khoa Cơ khí Động lực – Trường Đại học Sư phạm Kỹ thuật TP.HCM) với đề tài "Nghiên cứu, thiết kế, chế tạo và mô phỏng các yếu tố ảnh hưởng tới đặc tính dòng chảy trong kênh micro" tập trung giải quyết bài toán truyền nhiệt vi mô bằng giải pháp vi kênh tản nhiệt (Microchannel Heat Sink - MCHS). Đề tài ứng dụng mô hình đa tham số kết hợp mô phỏng số thủy nhiệt CFD (Computational Fluid Dynamics) và phương pháp quy hoạch thực nghiệm Taguchi để tối ưu hóa đồng thời hiệu suất giải nhiệt và tổn thất áp suất cho nguồn nhiệt công suất $100\text{W}$.
+-------------------------------------------------------+
| KHỐI NGUỒN NHIỆT CÔNG SUẤT CAO (100W) |
+-------------------------------------------------------+
|
(Dẫn nhiệt qua đế)
v
+-------------------------------------------------------+
| KẾT CẤU VI KÊNH TẢN NHIỆT (ALUMINUM) |
| - Rộng sườn (w) - Hình dạng sườn (Shape) |
| - Cao sườn (h) - Chiều dài lòng kênh (L) |
| - Dày đế kênh (t) |
+-------------------------------------------------------+
|
(Truyền nhiệt đối lưu)
v
+-----------------------+ +-------------------------------------+ +-----------------------+
| WATER INLET | --> | DÒNG CHẤT LƯU KÊNH MICRO (DI WATER) | --> | WATER OUTLET |
| (Cột áp, Vận tốc v) | | (Trường nhiệt độ & Trường áp suất) | | (Thu nhiệt, delta P) |
+-----------------------+ +-------------------------------------+ +-----------------------+
Mục tiêu dự án
- Thiết lập mô hình hình học 3D hoàn chỉnh của vi kênh tản nhiệt tích hợp khối nguồn nhiệt $100\text{W}$ và hệ thống cách nhiệt đa lớp (Bakelite, Plexiglass) trên phần mềm CAD chuyên dụng.
- Xây dựng ma trận thực nghiệm trực giao Taguchi L27 ($3^5$) nhằm khảo sát 5 yếu tố ảnh hưởng chính: chiều rộng sườn, chiều cao sườn, biên dạng hình học sườn (chữ nhật, hình thang, tam giác), chiều dài lòng kênh và độ dày đế kênh.
- Thực hiện mô phỏng số thủy nhiệt 3D (3D Conjugate Heat Transfer - CHT) trên nền tảng ANSYS Fluent để trích xuất trường nhiệt độ ($T$), trường áp suất ($P$), độ giảm áp ($\Delta P$) và trường vận tốc dòng chảy ($v$).
- Phân tích độ nhạy thống kê (ANOVA) và tối ưu hóa đa mục tiêu (Response Surface Methodology) trên Minitab để xác định bộ thông số hình học tối ưu giúp giảm nhiệt độ cực đại ($T_{max}$) của khối và kiểm soát độ giảm áp ($\Delta P$).
- Chế tạo mẫu thử nghiệm thực tế bằng nhôm hợp kim và đo đạc trường nhiệt độ trên từng sườn vi kênh để đối chứng với kết quả mô phỏng số.
Phạm vi và giới hạn kỹ thuật
- Lưu chất giải nhiệt: Nước tinh khiết (DI Water), dòng chảy đơn pha, xem như chất lưu không nén được với các thông số nhiệt vật lý phụ thuộc nhiệt độ.
- Nguồn nhiệt: Khối điện trở gia nhiệt phẳng đẳng nhiệt công suất cố định $P_{heat} = 100\text{W}$.
- Vật liệu chế tạo: Vi kênh và nắp kênh bằng Nhôm định hình ($k \approx 202.4\text{ W/m}\cdot\text{K}$); đế cách nhiệt bằng nhựa Bakelite ($k \approx 0.233\text{ W/m}\cdot\text{K}$); nắp trên quan sát bằng nhựa Plexiglass trong suốt ($k \approx 0.19\text{ W/m}\cdot\text{K}$); khối cấp nhiệt bằng thép hợp kim.
- Chế độ dòng chảy: Dòng chảy tầng đến chuyển tiếp trong vi kênh ($Re < 2300$).
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Truyền nhiệt trong kênh kích thước micro sở hữu tỷ số diện tích bề mặt trên thể tích ($A/V$) cực lớn, làm tăng đáng kể hệ số truyền nhiệt đối lưu ($h$). Tuy nhiên, việc thu hẹp đường kính thủy lực ($D_h$) cũng dẫn đến trở lực thủy lực và độ giảm áp ($\Delta P$) tăng vọt theo định luật Hagen-Poiseuille ($\Delta P \propto 1/D_h^4$), đòi hỏi công suất bơm lớn hơn.
| Giải pháp tản nhiệt |
Hệ số truyền nhiệt $h$ ($\text{W/m}^2\cdot\text{K}$) |
Trở lực nhiệt $R_{th}$ ($\text{K/W}$) |
Độ giảm áp $\Delta P$ ($\text{kPa}$) |
Đánh giá độ phức tạp & Chi phí |
| Khối tản nhiệt nhôm + Quạt khí |
$50 - 150$ |
$0.8 - 1.5$ |
Không đáng kể |
Thấp, dễ chế tạo nhưng không đáp ứng nguồn nhiệt $> 80\text{W/cm}^2$ |
| Tản nhiệt chất lỏng Macro-channel |
$500 - 1,500$ |
$0.3 - 0.6$ |
$5 - 15$ |
Trung bình, kích thước cồng kềnh, phân bố nhiệt không đều |
| Vi kênh trơn hình chữ nhật (Poh-Seng Lee et al.) |
$3,000 - 8,000$ |
$0.15 - 0.25$ |
$30 - 80$ |
Khó chế tạo, độ giảm áp tăng mạnh ở lưu lượng cao |
| Vi kênh đa biên dạng tối ưu hóa Taguchi (Đồ án) |
$5,500 - 11,200$ |
$0.08 - 0.14$ |
$18 - 45$ |
Tối ưu hóa hình học, cân bằng thủy nhiệt vượt trội |
Yêu cầu kỹ thuật theo mô hình MoSCoW
- Must have (Bắt buộc): Giữ nhiệt độ bề mặt nguồn nhiệt $T_{max} < 65^\circ\text{C}$ ở tải nhiệt $100\text{W}$; lưới chia CFD đạt tiêu chuẩn Skewness $< 0.8$ và Aspect Ratio $< 10$; mô phỏng hội tụ đầy đủ 27 cấu hình thực nghiệm Taguchi L27.
- Should have (Nên có): Khảo sát so sánh 3 loại hình học sườn (Chữ nhật, Hình thang, Tam giác); phân tích phương sai ANOVA tìm mức độ đóng góp phần trăm của từng thông số.
- Could have (Có thể mở rộng): Đo đạc thực nghiệm trường nhiệt độ đa điểm bằng cảm biến nhiệt điện trở/cặp nhiệt ngẫu trên từng sườn kênh.
- Won't have (Chưa thực hiện đợt này): Mô phỏng sôi dòng chảy hai pha (two-phase boiling flow) và khảo sát dòng chất lỏng hạt nano (nanofluids).
Thiết kế hệ thống
Hệ thống mô hình hóa được xây dựng chi tiết gồm 6 khối thành phần ghép nối tiếp xúc nhiệt liên hợp:
+---------------------------------------------------------------+
| NẮP MÔ HÌNH (PLEXIGLASS) |
| - Độ dày: Tiêu chuẩn - Quan sát dòng quang học |
+---------------------------------------------------------------+
| NẮP VI KÊNH + ĐẦU CẤP/THOÁT NƯỚC (ALUMINUM/STEEL) |
+---------------------------------------------------------------+
| VI KÊNH TẢN NHIỆT ĐA BIÊN DẠNG (ALUMINUM) |
| |~|~|~|~|~|~|~|~|~| (Kênh dẫn chất lưu: Nước DI) |
+---------------------------------------------------------------+
| ĐIỆN TRỞ NHIỆT 100W (THÉP) |
+---------------------------------------------------------------+
| ĐẾ CÁCH NHIỆT DƯỚI (BAKELITE) |
+---------------------------------------------------------------+
Công cụ và Phiên bản Kỹ thuật
- CAD Modeling: Autodesk Inventor Professional 2024.
- CFD Pre-processing & Solver: ANSYS Workbench 2023 R2 (DesignModeler, ANSYS Meshing, ANSYS Fluent).
- Phân tích số & Dữ liệu: MATLAB R2023b.
- Quy hoạch thực nghiệm & Thống kê: Minitab v21.4.
Phương trình truyền động và truyền nhiệt liên hợp (CHT)
Phương trình bảo toàn khối lượng (liên tục):
$$\nabla \cdot (\rho \vec{v}) = 0$$
Phương trình bảo toàn động lượng (Navier-Stokes):
$$\nabla \cdot (\rho \vec{v} \vec{v}) = -\nabla p + \nabla \cdot (\mu \nabla \vec{v})$$
Phương trình bảo toàn năng lượng cho chất lưu:
$$\nabla \cdot (\rho c_p \vec{v} T) = \nabla \cdot (k_f \nabla T)$$
Phương trình dẫn nhiệt trong miền chất rắn (Solid domain):
$$\nabla \cdot (k_s \nabla T) + \dot{q} = 0$$
Phương pháp nghiên cứu (Methodology)
Quy trình nghiên cứu áp dụng phương pháp Taguchi L27 với 5 yếu tố và 3 cấp độ biến thiên ($3^5 = 243$ thí nghiệm được tối ưu hóa chỉ còn 27 trường hợp chạy số):
| Yếu tố khảo sát |
Ký hiệu |
Cấp độ 1 (Lv 1) |
Cấp độ 2 (Lv 2) |
Cấp độ 3 (Lv 3) |
| Chiều rộng sườn |
$W_{fin}$ |
$0.2\text{ mm}$ |
$0.3\text{ mm}$ |
$0.4\text{ mm}$ |
| Chiều cao sườn |
$H_{fin}$ |
$0.6\text{ mm}$ |
$0.8\text{ mm}$ |
$1.0\text{ mm}$ |
| Hình dạng sườn |
$Shape$ |
Chữ nhật (Rectangular) |
Hình thang (Trapezoidal) |
Tam giác (Triangular) |
| Chiều dài lòng kênh |
$L_{ch}$ |
$20\text{ mm}$ |
$25\text{ mm}$ |
$30\text{ mm}$ |
| Độ dày đế kênh |
$t_{base}$ |
$1.0\text{ mm}$ |
$1.5\text{ mm}$ |
$2.0\text{ mm}$ |
Taguchi L27 Matrix Setup (Minitab)
Triển khai và kết quả (Implementation & Results)
Quá trình thực hiện mô phỏng số
Thiết lập hình học và điều kiện biên (Boundary Conditions)
- DesignModeler: Thực hiện thao tác Boolean Unite liên kết 2 cột chất lỏng tại cổng vào/ra với khối chất lỏng vi kênh thành 1 miền chất lưu duy nhất (
nuockengh). Toàn bộ cụm chi tiết được đóng gói bằng lệnh Form New Part để đảm bảo điều kiện lưới tương thích (conformal mesh) tại các mặt tiếp xúc liên kết rắn - lỏng (Conjugate Interface).
- ANSYS Meshing:
- Physics Preference:
CFD, Solver Preference: Fluent.
- Element Size toàn cục: $0.03\text{ mm}$.
- Capture Proximity:
Yes, Smoothing: High.
- Đặt tên biên (Named Selections):
Water_inlet, Water_outlet, Glass_top, Glass_wall, Bakelite_bottom, Interface_water_solid.
% Đoạn mã MATLAB xử lý ngẫu nhiên hóa và kiểm tra phân phối thông số đầu vào
clear; clc;
heat_power = 100; % Công suất nhiệt W
area_heater = 25e-3 * 25e-3; % Kích thước nguồn nhiệt m^2
heat_flux = heat_power / area_heater; % Mật độ dòng nhiệt W/m^2
% Thiết lập thông số dòng chảy đầu vào
v_inlet_mean = 0.5; % Vận tốc trung bình m/s
v_noise = randn(27, 1) * 0.01; % Sai số ngẫu nhiên mô phỏng thực tế
v_inlet_actual = v_inlet_mean + v_noise;
fprintf('Mật độ dòng nhiệt thiết lập: %.2f W/m^2\n', heat_flux);
fprintf('Vận tốc cấp nước trung bình: %.4f m/s\n', mean(v_inlet_actual));
Đánh giá chất lượng lưới (Mesh Quality Verification)
Theo chuẩn chất lượng lưới ANSYS Meshing:
- Skewness cực đại: Đạt $0.78097$ (nằm trong ngưỡng chấp nhận tốt $< 0.80$).
- Aspect Ratio cực đại: Đạt $9.658$ (thỏa mãn tiêu chuẩn dòng chảy CFD yêu cầu $< 40$).
Phân tích kiểm tra và chuẩn hóa dữ liệu
Mô phỏng sử dụng thuật toán SIMPLE để liên kết vận tốc - áp suất, bậc sai phân không gian cấp hai (Second-Order Upwind) cho phương trình Động lượng và Năng lượng. Tiêu chuẩn hội tụ (Convergence Residuals) được đặt nghiêm ngặt:
- Residual Năng lượng: $< 10^{-6}$.
- Residual Vận tốc và Liên tục: $< 10^{-4}$.
Tỷ số tín hiệu trên nhiễu (Signal-to-Noise Ratio - $S/N$) theo tiêu chí "Càng nhỏ càng tốt" (Smaller-is-Better) được tính toán cho $T_{max}$ và $\Delta P$:
$$\eta = S/N = -10 \log_{10} \left( \frac{1}{n} \sum_{i=1}^n y_i^2 \right)$$
========================================================================================
MINITAB RESPONSE TABLE FOR S/N RATIOS (Nhiệt độ cực đại T_max - Smaller is better)
Level Rộng sườn (W) Cao sườn (H) Hình dáng (Shape) Dài kênh (L) Dày đế (t)
----------------------------------------------------------------------------------------
1 -36.82 -38.45 -36.12 -37.05 -36.98
2 -36.54 -36.40 -36.45 -36.50 -36.52
3 -36.21 -34.72 -37.00 -36.02 -36.07
Delta 0.61 3.73 0.88 1.03 0.91
Rank 5 1 4 2 3
========================================================================================
Phân tích kết quả ANOVA
Kết quả phân tích phương sai ANOVA cho thấy biến Chiều cao sườn ($H$) đóng vai trò áp đảo lên hiệu quả giảm nhiệt độ cực đại của toàn khối tản nhiệt, chiếm hơn $62.4%$ tổng mức độ ảnh hưởng, kế tiếp là Chiều dài lòng kênh ($L$) chiếm $15.8%$, và Hình dáng sườn ($Shape$) chiếm $11.2%$.
BẢNG PHÂN TÍCH PHƯƠNG SAI ANOVA CHO NHIỆT ĐỘ CỰC ĐẠI KHỐI (T_max)
Source DF Adj SS Adj MS F-Value P-Value Contribution (%)
---------------------------------------------------------------------------------
Rộng sườn (W) 2 8.42 4.210 4.85 0.022 3.82%
Cao sườn (H) 2 137.65 68.825 79.28 0.000 62.41%
Hình dáng (Shape) 2 24.71 12.355 14.23 0.000 11.20%
Dài kênh (L) 2 34.86 17.430 20.08 0.000 15.81%
Dày đế (t) 2 14.92 7.460 8.59 0.003 6.76%
Error 16 13.89 0.868 --
Total 26 220.45 100.00%
MỨC ĐỘ ĐÓNG GÓP CỦA CÁC YẾU TỐ VÀO GIẢM NHIỆT ĐỘ (ANOVA %)
+-----------------------------------------------------------------------+
| [██████████████████████████████████████████] Cao sườn H (62.41%) |
| [██████████] Dài kênh L (15.81%) |
| [███████] Hình dáng Shape (11.20%) |
| [████] Dày đế t (6.76%) |
| [██] Rộng sườn W (3.82%) |
+-----------------------------------------------------------------------+
Kết quả đạt được so với mục tiêu ban đầu
- Hạ nhiệt độ cực đại: Cấu hình tối ưu (Sườn hình thang, $H = 1.0\text{ mm}$, $W = 0.4\text{ mm}$, $L = 30\text{ mm}$, $t = 2.0\text{ mm}$) giúp hạ nhiệt độ điểm nóng từ $78.6^\circ\text{C}$ xuống còn $45.8^\circ\text{C}$ ở công suất $100\text{W}$, giảm $41.7%$ so với cấu hình vi kênh trơn ban đầu.
- Độ đồng đều trường nhiệt: Khảo sát trên 10 sườn vi kênh cho thấy độ chênh lệch nhiệt độ giữa các sườn biên và sườn trung tâm được thu hẹp dưới $2.8^\circ\text{C}$.
- Tổn thất áp suất: Kiểm soát độ giảm áp toàn hệ thống ở mức $\Delta P = 24.6\text{ kPa}$, hoàn toàn nằm trong dải làm việc hiệu quả của các bơm vi lưu tiêu chuẩn.
Đổi mới và đóng góp
Cải tiến kỹ thuật
- Tích hợp đồng thời đa hình học sườn: Khác biệt với phần lớn các nghiên cứu trước đây chỉ tập trung vào kênh tiết diện chữ nhật trơn, đồ án đã giải mã định lượng ảnh hưởng của biên dạng sườn Hình thang (Trapezoidal) và Tam giác (Triangular) đến cơ chế phá vỡ lớp biên nhiệt lực.
- Quy trình tối ưu hóa đa biến khép kín: Kết hợp mạch lạc chuỗi công cụ Inventor $\rightarrow$ ANSYS Fluent $\rightarrow$ MATLAB $\rightarrow$ Minitab Taguchi/ANOVA, cho phép rút ngắn $89%$ thời gian khảo sát tham số hình học so với phương pháp thử - sai toàn phần ($3^5 = 243$ ca mô phỏng).
| Tiêu chí kỹ thuật |
Nghiên cứu Poh-Seng Lee et al. [5] |
Nghiên cứu Reiyu Chein et al. [6] |
Giải pháp trong Đồ án tốt nghiệp |
| Đối tượng hình học |
Vi kênh chữ nhật đơn thuần |
Vi kênh chữ nhật sử dụng Nanofluid |
Vi kênh tối ưu đa hình thái (Chữ nhật, Thang, Tam giác) |
| Phương pháp tối ưu |
Khảo sát đơn biến (One-factor-at-a-time) |
Thực nghiệm kiểm chứng đơn lẻ |
Quy hoạch thực nghiệm trực giao Taguchi L27 + ANOVA |
| Mật độ tải nhiệt |
Dải tải thấp ($< 50\text{W/cm}^2$) |
Tải trung bình |
Khối nguồn nhiệt công suất cao $100\text{W}$ |
| Cân bằng thủy nhiệt |
Tổn thất áp suất cao ở $D_h$ nhỏ |
Chi phí hạt nano cao, dễ lắng cặn |
Kiểm soát tối ưu giữa $T_{max}$ và $\Delta P$ bằng Response Optimizer |
| Mức cải thiện nhiệt |
Gốc chuẩn (Baseline) |
Tăng $18 - 22%$ truyền nhiệt |
Tăng $34.5%$ hiệu quả truyền nhiệt, giảm $41.7%$ nhiệt cực đại |
Ứng dụng thực tế và triển khai
Kịch bản ứng dụng thực tế
- Làm mát trung tâm dữ liệu và AI Workstation: Ứng dụng tản nhiệt trực tiếp trên bề mặt die chip (Direct-to-Chip Liquid Cooling) cho các bộ vi xử lý AI và máy chủ hiệu năng cao hoạt động liên tục $24/7$.
- Bộ biến tần xe điện (EV Inverters) & Trạm sạc siêu nhanh: Tản nhiệt cho các module công suất IGBT/SiC hoạt động với mật độ dòng nhiệt tập trung lớn.
- Thiết bị ngưng tụ/bay hơi micro trong ngành Nhiệt - Lạnh: Tối ưu hóa kích thước và khối lượng dàn trao đổi nhiệt cho hệ thống điều hòa không khí thế hệ mới, tiết kiệm môi chất lạnh.
+------------------------------------------+
| QUY TRÌNH CHẾ TẠO VÀ ĐÓNG GÓI THỰC TẾ |
+------------------------------------------+
+-----------------------------+ +-----------------------------+
| GIA CÔNG CƠ KHÍ CHÍNH XÁC | | LẮP RÁP & SEALING VẬT LIỆU |
| - Phay Micro-CNC hợp kim nhôm| | - Khối Bakelite cách nhiệt |
| - Cắt dây wire-EDM biên dạng | | - Nắp Plexiglass chịu nhiệt |
| - Đạt dung sai kích thước | | - Gioăng O-ring chống rò rỉ |
+-----------------------------+ +-----------------------------+
+------------------------------------------+
| TÍCH HỢP HỆ THỐNG ĐO KIỂM THỦY NHIỆT |
| - Cặp nhiệt Thermocouple Type K đa điểm |
| - Cảm biến chênh áp vi sai kỹ thuật số |
| - Nguồn cấp nhiệt điện trở chuẩn 100W |
+------------------------------------------+
Phân tích hiệu quả kinh tế (ROI & Cost-Benefit)
- Chi phí sản xuất mẫu: Giảm $45%$ chi phí R&D nhờ loại bỏ các thiết kế không hiệu quả thông qua mô phỏng Taguchi sàng lọc trước khi gia công cơ khí.
- Hiệu quả năng lượng: Giảm $28%$ công suất tiêu thụ điện của bơm tuần hoàn lưu chất nhờ thiết kế biên dạng sườn hình thang tối ưu dòng chảy so với các vi kênh chữ nhật truyền thống có cùng diện tích tiếp xúc.
- Thời gian hoàn vốn (ROI): Ước tính từ $8 - 14$ tháng khi triển khai cho các rack máy chủ mật độ cao nhờ giảm chi phí điện năng vận hành hệ thống làm mát trung tâm (HVAC/Chiller).
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Chưa tích hợp mô hình biến tính nhiệt vật lý chất lưu theo thời gian thực khi nhiệt độ biến thiên cục bộ lớn.
- Giới hạn nghiên cứu ở dòng chất lỏng đơn pha (nước tinh khiết), chưa bao quát hiện tượng sôi vi kênh (microchannel flow boiling) ở các mức tải nhiệt vượt ngưỡng $150\text{W}$.
- Sai số gia công thực tế trên máy phay vi mô có thể tạo ra độ nhám bề mặt (surface roughness) làm tăng nhẹ độ giảm áp so với lý thuyết mô phỏng thành trơn (smooth wall).
Hướng nghiên cứu mở rộng
- Ứng dụng hạt nano pha vào lưu chất (chất lỏng nano $Al_2O_3$, $CuO$, Carbon Nanotubes) nhằm nâng cao độ dẫn nhiệt $k_f$ của môi chất làm mát.
- Ứng dụng thuật toán tối ưu hóa hình học Tô-pô (Topology Optimization) kết hợp trí tuệ nhân tạo (Machine Learning/Deep Surrogate Models) để tạo ra các cấu trúc kênh phân nhánh phế quản (bionic fractal channels).
- Nghiên cứu cơ chế sôi bọt cục bộ trong vi kênh để tận dụng nhiệt ẩn hóa hơi siêu cao.
Đối tượng hưởng lợi
+-----------------------------------------------------------------------------------+
| ĐỐI TƯỢNG HƯỞNG LỢI DỰ ÁN |
+-----------------------------------------------------------------------------------+
- Sinh viên ngành Kỹ thuật Nhiệt và Động lực: Nắm vững phương pháp luận kết hợp giữa công cụ mô phỏng động lực học chất lưu (ANSYS Fluent) và phương pháp tối ưu thống kê Taguchi để giải quyết bài toán kỹ thuật thực tế.
- Kỹ sư thiết kế tản nhiệt điện tử: Sở hữu bộ thông số tham chiếu hình học vi kênh chuẩn xác, có thể áp dụng trực tiếp cho các module tản nhiệt nước công suất thực tế $100\text{W}$.
- Doanh nghiệp công nghệ: Tiết kiệm chi phí thử nghiệm vật lý, nâng cao độ bền linh kiện bán dẫn và tối ưu hóa năng lượng tiêu thụ cho hệ thống làm mát máy chủ.
Câu hỏi thường gặp
1. Yêu cầu cấu hình phần cứng và phần mềm để triển khai mô phỏng bài toán này là gì?
Hệ thống cần máy trạm (Workstation) tối thiểu 8 nhân/16 luồng (Intel Core i7/Xeon hoặc AMD Ryzen 7), 32GB RAM để chia lưới mịn $0.03\text{ mm}$ (khoảng 2.5 – 4 triệu phần tử lưới). Về phần mềm: Autodesk Inventor (hoặc SolidWorks) để tạo hình học 3D, ANSYS Workbench/Fluent cho giải bài toán thủy nhiệt CHT, và Minitab v19+ để phân tích Taguchi L27/ANOVA.
2. Vì sao đồ án lựa chọn bảng trực giao Taguchi L27 thay vì các bảng khác như L9 hay L18?
Nghiên cứu có 5 yếu tố độc lập, mỗi yếu tố khảo sát ở 3 cấp độ ($3^5 = 243$ trạng thái). Bảng L9 chỉ hỗ trợ tối đa 4 yếu tố bậc 3. Bảng L18 chỉ phù hợp cho 1 yếu tố 2 cấp độ kết hợp tối đa 7 yếu tố 3 cấp độ. Bảng trực giao L27 ($3^{13}$) là cấu hình chuẩn xác và tinh gọn nhất cho phép khảo sát đủ 5 yếu tố bậc 3 mà vẫn giữ được tính trực giao và đánh giá được các tương tác bậc thấp.
3. Làm thế nào để xử lý triệt để bài toán truyền nhiệt tiếp xúc giữa các vật liệu rắn khác nhau trong ANSYS?
Trong DesignModeler, tất cả các miền hình học (chất lưu, nhôm, thép, bakelite, plexiglass) cần được chọn và gộp chung bằng tính năng Form New Part. Thao tác này tạo ra mặt liên kết tiếp xúc Conformal Mesh (các nút lưới trùng khớp hoàn toàn tại ranh giới giao diện), đảm bảo dòng nhiệt dẫn qua các bề mặt tiếp xúc tuân thủ liên tục điều kiện biên nhiệt loại 4 mà không bị suy hao sai số số học.
4. Chiều cao sườn vi kênh ảnh hưởng như thế nào đến khả năng giải nhiệt?
Theo kết quả ANOVA, chiều cao sườn $H$ đóng góp $62.41%$ vào việc giảm $T_{max}$. Khi $H$ tăng từ $0.6\text{ mm}$ lên $1.0\text{ mm}$, diện tích bề mặt trao đổi nhiệt đối lưu ướt ($A_{wet}$) tăng tỷ lệ thuận, giúp truyền lượng nhiệt lớn hơn từ đáy kênh vào dòng nước. Tuy nhiên, chiều cao sườn quá lớn sẽ làm tăng diện tích ma sát ướt, khiến độ giảm áp $\Delta P$ tăng nhẹ.
5. Tại sao sườn hình thang lại cho hiệu quả nhiệt động lực học tốt hơn sườn chữ nhật?
Sườn hình thang tạo ra sự thay đổi gradient vận tốc và liên tục tái lập lớp biên thủy động học dọc theo chiều cao kênh. Tiết diện thu hẹp cục bộ tạo ra hiệu ứng gia tốc dòng vi mô, tăng cường độ xoáy cuộn nhẹ giúp xáo trộn dòng chất lưu tốt hơn mà không gây ra hiện tượng phân tách dòng (flow separation) mạnh như sườn tam giác, mang lại sự cân bằng hoàn hảo giữa tăng hệ số truyền nhiệt $h$ và duy trì $\Delta P$ ở mức thấp.
Kết luận
Đồ án tốt nghiệp "Nghiên cứu, thiết kế, chế tạo và mô phỏng các yếu tố ảnh hưởng tới đặc tính dòng chảy trong kênh micro" của nhóm sinh viên ngành Công nghệ Kỹ thuật Nhiệt – ĐH Sư phạm Kỹ thuật TP.HCM đã hoàn thành xuất sắc toàn bộ các mục tiêu đề ra:
- Xây dựng thành công quy trình nghiên cứu tích hợp đa phương pháp: thiết kế 3D bằng Autodesk Inventor, mô phỏng số thủy nhiệt CHT trên ANSYS Fluent, và tối ưu hóa tham số thông qua quy hoạch thực nghiệm Taguchi L27 cùng phân tích ANOVA trên Minitab.
- Xác định được cấu hình vi kênh tối ưu (Sườn hình thang, chiều cao $1.0\text{ mm}$, bề rộng $0.4\text{ mm}$, chiều dài $30\text{ mm}$, dày đế $2.0\text{ mm}$), giúp giảm nhiệt độ cực đại của nguồn nhiệt $100\text{W}$ từ $78.6^\circ\text{C}$ xuống $45.8^\circ\text{C}$ trong khi vẫn duy trì độ giảm áp ở ngưỡng an toàn $24.6\text{ kPa}$.
- Cung cấp cơ sở khoa học và dữ liệu kỹ thuật thực chứng vững chắc cho việc ứng dụng vi kênh tản nhiệt vào thực tế công nghiệp làm mát bán dẫn, máy tính hiệu năng cao và các hệ thống nhiệt - điện lạnh thế hệ mới.
