Bối cảnh và vấn đề nghiên cứu

Sự phát triển nhanh chóng của công nghệ vi điện tử đặt ra yêu cầu nghiêm ngặt về hiệu năng tính toán và mức độ thu nhỏ kích thước của các vi mạch tích hợp mật độ siêu cao (VLSI). Các mạch xử lý tốc độ cao tạo ra mật độ công suất nhiệt cục bộ rất lớn. Các giải pháp tản nhiệt truyền thống bằng không khí dần chạm tới ngưỡng giới hạn vật lý khi mật độ dòng nhiệt vượt ngưỡng $20\text{ W/cm}^2$, trong khi các cấu kiện hiện đại có thể đòi hỏi giải tỏa nguồn nhiệt lên tới hơn $1000\text{ W/cm}^2$. Nhằm giải quyết thách thức này, giải pháp tản nhiệt vi kênh (microchannel heat sink) đã được đề xuất lần đầu bởi Tuckerman và Pease vào năm 1981, ứng dụng chất lỏng làm mát cưỡng bức chảy qua các rãnh dẫn có kích thước hiển vi.

Thực tiễn nghiên cứu cho thấy hiệu quả truyền nhiệt và phân bố dòng chảy trong kênh vi mô chịu sự chi phối phức tạp từ nhiều yếu tố: hình học mặt cắt kênh, cách bố trí kênh, lưu chất giải nhiệt, lưu lượng và điều kiện biên nhiệt. Việc nghiên cứu tối ưu hóa thiết kế trước đây thường dựa trên phương pháp thử sai từng biến riêng lẻ, gây tốn kém thời gian và chi phí chế tạo mẫu thử. Đồ án "Nghiên cứu và mô phỏng số các yếu tố ảnh hưởng tới đặc tính dòng chảy trong kênh micro dạng lai" được nhóm sinh viên ngành Công nghệ Kỹ thuật Nhiệt, Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh thực hiện nhằm giải quyết bài toán tối ưu hóa hình học cho bộ tản nhiệt vi kênh, nâng cao hiệu quả làm mát và giảm tổn thất áp suất trước khi đưa vào sản xuất thực nghiệm.

Mục tiêu và nhiệm vụ nghiên cứu của đồ án bao gồm:

  1. Nghiên cứu tổng quan về cơ chế dòng chảy, sự sụt giảm áp suất và đặc tính truyền nhiệt trong các kênh tản nhiệt vi mô.
  2. Xây dựng mô hình hình học 3D của cụm tản nhiệt vi kênh với các dạng biên dạng mặt cắt khác nhau (chữ nhật, hình thang, hình tam giác).
  3. Thiết lập ma trận quy hoạch thực nghiệm theo phương pháp Taguchi kết hợp với phân tích phương sai (ANOVA) nhằm đánh giá mức độ tác động của các thông số hình học: chiều rộng ($W$), chiều cao ($H$) và hình dạng tiết diện ($S$).
  4. Thực hiện mô phỏng số động lực học chất lưu (CFD) trên phần mềm ANSYS Fluent để thu thập các trường nhiệt độ, trường áp suất và trường vận tốc.
  5. Xử lý số liệu mô phỏng, xác định cấu hình kích thước tối ưu theo các tiêu chí vận hành và đối chiếu kết quả với mẫu gia công cơ khí chính xác.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Cụm tản nhiệt kênh micro (microchannel heat sink) dạng lai với ba kiểu hình dạng tiết diện rãnh/gân: hình chữ nhật (Rectangular), hình tam giác (Triangular) và hình thang (Trapezoidal); cùng các đại lượng vật lý dòng chảy gồm nhiệt độ ($T$), áp suất ($P$) và vận tốc ($V$).
  • Phạm vi nghiên cứu:
    • Vật liệu chế tạo khối vi kênh là nhôm nguyên chất; cấu trúc bao che gồm nhựa acrylic/glass (PVC) và tấm cách nhiệt bakelite.
    • Lưu chất giải nhiệt sử dụng là nước lỏng (water liquid) chảy tầng một pha.
    • Các thông số đầu vào được khảo sát trong phạm vi ba mức kích thước hình học và được mô phỏng tĩnh trong môi trường ANSYS Fluent.

Cơ sở lý thuyết và phương pháp

Nghiên cứu được xây dựng trên nền tảng cơ học chất lưu và truyền nhiệt trong các kênh kích thước nhỏ, kết hợp với các công cụ quy hoạch thực nghiệm thống kê hiện đại:

  • Khung lý thuyết và mô hình tính toán:

    • Động học dòng chảy và sụt áp trong vi kênh: Lý thuyết sụt áp một pha qua vi kênh chỉ ra rằng tổn thất áp suất tỷ lệ nghịch với đường kính thủy lực ($D_h$). Do $D_h$ rất nhỏ, sụt áp trên một đơn vị chiều dài kênh lớn hơn nhiều so với kênh thông thường, dẫn đến công suất bơm tiêu thụ tăng cao.
    • Hiện tượng bất ổn định dòng chảy: Văn bản tham chiếu các nghiên cứu của Ledinegg, Bogojevic, Kaun, Kandlikar và Szczukiewicz về hiện tượng đảo chiều dòng chảy, dao động áp suất và phân bố không đều giữa các kênh song song.
    • Kế thừa các công trình trong nước: Nghiên cứu của Đặng Thành Trung về dòng chảy bay hơi vi kênh; Đặng Hùng Sơn về thiết kế sáu sigma (DFSS) cho tản nhiệt vi kênh hai lớp; Đoàn Minh Hùng về quá trình ngưng tụ; Lê Bá Tân về lưu lượng dòng trong kênh vuông; Nguyễn Trọng Hiếu về vi kênh môi chất $\text{CO}_2$; Nguyễn Huy Bích về phương pháp phần tử hữu hạn phân tích vi giọt trên COMSOL Multiphysics.

  • Phương pháp nghiên cứu:

    • Phương pháp Taguchi (Taguchi Methods): Sử dụng bảng ma trận trực giao (Orthogonal Arrays) $L_9$ để giảm thiểu số lượt thử nghiệm mô phỏng cần thiết. Phương pháp đánh giá dựa trên giá trị trung bình (Mean) và tỷ số tín hiệu trên nhiễu ($S/N$) nhằm xác định cấu hình ít nhạy cảm nhất với các biến động nhiễu môi trường.
    • Phân tích phương sai (Analysis of Variance - ANOVA): Phân tách tổng biến thiên của các biến phụ thuộc ($\Delta T$, $\Delta P$, $\Delta V$) thành các phần tương ứng với các biến giải thích độc lập (chiều rộng $W$, chiều cao $H$, tiết diện $S$) và phần dư không giải thích được, từ đó lượng hóa tỷ lệ đóng góp phần trăm và mức độ tin cậy thống kê của từng biến.
    • Phương pháp mô phỏng số (CFD): Giải hệ phương trình vi phân đạo hàm riêng bảo toàn khối lượng (liên tục), bảo toàn động lượng (Navier-Stokes) và bảo toàn năng lượng trong môi trường ba chiều bằng ANSYS Fluent.

  • Nguồn dữ liệu và công cụ triển khai:

    • Autodesk Inventor: Thiết kế mô hình hình học 3D chi tiết của cụm vi kênh, nắp kênh, điện trở cấp nhiệt, các khối chất lỏng và lớp vỏ cách nhiệt.
    • MATLAB: Lập trình hàm ngẫu nhiên randn để mô phỏng sai số đầu vào thực tế (nhiễu 1%) cho nhiệt độ môi trường xung quanh giá trị $300\text{ K}$, mật độ dòng nhiệt quanh giá trị $2.000.000\text{ W/m}^2$ và vận tốc dòng vào quanh giá trị $1{,}5\text{ m/s}$ (hoặc $1{,}015\text{ m/s}$).
    • ANSYS Fluent: Chia lưới sai số hữu hạn, kiểm tra chỉ số độ méo lưới (Skewness) và tỷ lệ khung hình (Aspect Ratio), thiết lập điều kiện biên và giải số trường dòng chảy.
    • Minitab 19: Tính toán các thông số thống kê ANOVA, đồ thị tác động Mean và $S/N$.

Nội dung chính theo từng chương

CHƯƠNG 1: TỔNG QUAN

Chương này trình bày cơ sở lịch sử và sự cần thiết của tản nhiệt vi kênh trong lĩnh vực kỹ thuật hiện đại. Tác giả nêu bật việc vi kênh được phát triển từ năm 1981 bởi Tuckerman và Pease để giải quyết bài toán tản nhiệt cho mạch tích hợp mật độ rất lớn (VLSI). Ứng dụng thực tiễn của công nghệ được minh họa qua các hệ thống làm mát card đồ họa máy tính (VGA Block EK-FC1080, Asus GTX-650), bộ tản nhiệt CPU và tản nhiệt ổ cứng SSD M.2. Tác giả phân tích tính cấp thiết của việc tìm hiểu các quy luật dòng chảy và truyền nhiệt khi kích thước kênh thu nhỏ xuống cấp micro, đồng thời xác lập đối tượng nghiên cứu (vi kênh nhôm tiết diện chữ nhật, tam giác, hình thang) và phạm vi nghiên cứu (lưu chất nước, mô phỏng số trên các biến hình học).

CHƯƠNG 2: CƠ SỞ LÝ THUYẾT

Chương 2 tổng hợp bức tranh học thuật chuyên sâu về vi kênh từ các công trình nghiên cứu trong nước (Đặng Thành Trung, Đặng Hùng Sơn, Đoàn Minh Hùng, Lê Bá Tân, Nguyễn Trọng Hiếu, Nguyễn Huy Bích) và các nghiên cứu quốc tế (Sobierska, Megahed, Zhang & Fu, Cornwell & Kew, Qu & Mudawar, Ergu, Peiyi & Little, Peng & Wang, v.v.). Nội dung tập trung vào bốn chủ đề: cấu trúc các dạng dòng chảy trong vi kênh, đặc tính sụt áp, hệ số truyền nhiệt đối lưu và các dạng bất ổn định nhiệt - thủy lực (như mất ổn định Ledinegg, dao động áp suất, đảo chiều dòng chảy).

Phần tiếp theo của chương đi sâu vào cơ sở toán học của phương pháp ANOVA và kỹ thuật Taguchi. Tác giả làm rõ triết lý hàm tổn thất chất lượng, kiểm soát chất lượng ngoại tuyến, nguyên lý trực giao của các bảng quy hoạch thực nghiệm và quy trình triển khai phân tích tỷ số $S/N$. Chương cũng giới thiệu chức năng và vai trò của hệ sinh thái phần mềm được ứng dụng: Autodesk Inventor, ANSYS Fluent, MATLAB và Minitab.

CHƯƠNG 3: XÂY DỰNG MODEL VÀ THIẾT LẬP THÔNG SỐ MÔ PHỎNG

Chương 3 chi tiết hóa quy trình thiết lập bài toán mô phỏng. Bắt đầu từ việc xác định 3 yếu tố khảo sát gồm Chiều rộng kênh $P_1$ ($W$), Chiều cao kênh $P_2$ ($H$) và Dạng mặt cắt $P_3$ ($S$ gồm Rec, Trape, Tri). Bảng quy hoạch thực nghiệm trực giao Taguchi $L_9$ được thiết lập với 9 mô hình thử nghiệm:

Số thứ tự (TN) Chiều rộng kênh $P_1$ ($W$) Chiều cao kênh $P_2$ ($H$) Dạng tiết diện $P_3$ ($S$)
1 Mức 1 ($0{,}7\text{ mm}$) Mức 1 Mức 1 (Rec)
2 Mức 1 ($0{,}7\text{ mm}$) Mức 2 Mức 2 (Trape)
3 Mức 1 ($0{,}7\text{ mm}$) Mức 3 Mức 3 (Tri)
4 Mức 2 Mức 1 Mức 2 (Trape)
5 Mức 2 Mức 2 Mức 3 (Tri)
6 Mức 2 Mức 3 Mức 1 (Rec)
7 Mức 3 Mức 1 Mức 3 (Tri)
8 Mức 3 Mức 2 Mức 1 (Rec)
9 Mức 3 Mức 3 Mức 2 (Trape)

Mô hình 3D được dựng trên Autodesk Inventor bao gồm các chi tiết: khối vi kênh bằng nhôm, nắp kênh nhôm, khối điện trở gia nhiệt, khối chất lỏng nước, tấm nhựa trong suốt acrylic/glass (PVC) và tấm đế cách nhiệt bằng nhựa Phenolic Bakelite. Nhóm tác giả sử dụng phần mềm MATLAB với hàm randn(1,3) để tạo sai số ngẫu nhiên 1% cho các điều kiện vận hành nhằm đưa vào mô hình tính toán:

  • Nhiệt độ môi trường: $x = 300 + 3 \times \text{randn}(1,3)\text{ [K]}$
  • Mật độ dòng nhiệt: $x = 2.000.000 + 20.000 \times \text{randn}(1,3)\text{ [W/m}^2\text{]}$
  • Vận tốc dòng vào: $x = 1{,}015 \times \text{randn}$ hoặc quanh giá trị danh định $1{,}5\text{ m/s}$.

Mô hình sau đó được đưa vào module Fluent Meshing để khởi tạo lưới phần tử thể tích, kiểm soát nghiêm ngặt các chỉ số chất lượng lưới gồm độ lệch méo (Skewness) và tỷ lệ cạnh (Aspect Ratio). Các điều kiện biên được gán cụ thể: bề mặt nguồn nhiệt (Heatflux), ngõ vào chất lỏng (Inlet), ngõ ra tự do (Outlet), các bề mặt vách bên và vách đáy cách nhiệt (Bakewall, Bakebottom, Bakehole, Glasswall, Glasstop, Portwall).

CHƯƠNG 4: KẾT QUẢ VÀ THẢO LUẬN

Chương 4 trình bày các dữ liệu thu nhận được sau quá trình giải lặp bằng phương pháp khởi tạo lai (Hybrid Initialization). Đối với từng trường hợp, tác giả trích xuất các phân bố đường đẳng nhiệt, trường áp suất và trường vector vận tốc (tại mặt phẳng đầu vào, bên trong kênh và mặt phẳng đầu ra Outlet).

Số liệu được phân tích theo 5 chỉ số đáp ứng: độ chênh lệch nhiệt độ bề mặt $\Delta T$, độ chênh nhiệt độ tổng thể $\Delta T_{\text{overall}}$, độ giảm áp suất $\Delta P$, độ biến thiên vận tốc dòng $\Delta V_{\text{overall}}$ và độ biến thiên vận tốc tại ngõ ra $\Delta V_{\text{outlet}}$. Thông qua Minitab 19, các biểu đồ phân tích phương sai ANOVA và đồ thị tỷ số $S/N$ được thiết lập nhằm lượng hóa thứ tự tác động của từng thông số hình học. Kết quả kiểm định cho thấy mô hình tối ưu được lựa chọn là cấu hình hình học thỏa mãn hài hòa giữa khả năng tản nhiệt cao và tổn thất áp suất thấp. Chương 4 cũng ghi nhận hình ảnh kênh micro thực tế sau khi gia công cơ khí chính xác tại Công ty TNHH TTL Precision.

CHƯƠNG 5: KẾT LUẬN VÀ KIẾN NGHỊ

Chương cuối tóm tắt các kết luận cốt lõi rút ra từ quá trình tính toán mô phỏng và phân tích thống kê; khẳng định tính khả thi của quy trình tích hợp giữa phương pháp Taguchi, ANOVA và ANSYS Fluent trong việc dự đoán đặc tính dòng nhiệt trong vi kênh. Tác giả đồng thời nêu các đề xuất, kiến nghị kỹ thuật để hoàn thiện mô hình trong các giai đoạn thử nghiệm thực tế tiếp theo.

Kết quả và đóng góp

Nghiên cứu đã lượng hóa sự ảnh hưởng của ba yếu tố hình học ($P_1$: Chiều rộng $W$, $P_2$: Chiều cao $H$, $P_3$: Dạng mặt cắt $S$) lên các đặc tính thủy lực và nhiệt của kênh vi mô. Kết quả phân tích thống kê xác định mức độ tác động và cấu hình tốt nhất theo từng tiêu chí đơn lẻ được tổng hợp như sau:

Tiêu chí khảo sát Yếu tố ảnh hưởng mạnh nhất Cấu hình mô hình được chọn
Chênh lệch nhiệt độ vi kênh ($\Delta T$) $P_2$ (Chiều cao kênh $H$) $W_2H_1S_2$
Biến thiên nhiệt độ tổng thể ($\Delta T_{\text{overall}}$) Không có yếu tố đơn lẻ vượt trội $W_2H_3S_1$
Độ giảm áp suất qua kênh ($\Delta P$) $P_1, P_3$ (Chiều rộng $W$ và Tiết diện $S$) $W_1H_1S_3$
Biến thiên vận tốc tổng thể ($\Delta V_{\text{overall}}$) $P_1, P_2, P_3$ (Cả ba yếu tố) $W_3H_2S_1$
Biến thiên vận tốc tại cửa ra ($\Delta V_{\text{outlet}}$) $P_1, P_2$ (Chiều rộng $W$ và Chiều cao $H$) $W_1H_1S_2$

Sau khi tổng hợp, cân bằng đa mục tiêu giữa tối thiểu hóa tổn thất áp suất và tối ưu hóa khả năng giải nhiệt, mô hình tối ưu tổng thể được nhóm nghiên cứu xác định bao gồm các thông số:

  • Chiều rộng kênh ($P_1$): $W = 0{,}8\text{ mm}$
  • Chiều cao kênh ($P_2$): $H = 0{,}32\text{ mm}$
  • Dạng mặt cắt tiết diện kênh ($P_3$): Dạng tam giác (Triangular)

Về mặt đóng góp và giải pháp kỹ thuật:

  • Đồ án đã xây dựng quy trình thiết kế kỹ thuật dựa trên chất lượng ngay từ giai đoạn mô hình hóa ảo: thay vì quét thông số rời rạc tốn kém tài nguyên máy tính, việc kết hợp ma trận trực giao Taguchi với ANOVA giúp rút gọn còn 9 ca mô phỏng nhưng vẫn nhận diện chính xác vai trò chi phối của chiều cao kênh đối với tản nhiệt và tiết diện đối với sụt áp.
  • Việc tích hợp hàm phân phối chuẩn từ MATLAB để đưa nhiễu ngẫu nhiên 1% vào điều kiện biên nhiệt và vận tốc ban đầu giúp phản ánh sát hơn các dao động vật lý trong môi trường làm việc thực tế của chip điện tử.
  • Đồ án đã liên hệ thực tế và bước đầu chuyển giao mô hình tối ưu sang công đoạn chế tạo thử nghiệm thông qua việc gia công cơ khí chính xác mẫu vi kênh bằng nhôm tại Công ty TNHH TTL Precision.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu:

  • Phạm vi mô phỏng chỉ mới dừng lại ở chế độ dòng chảy một pha đối với lưu chất nước lỏng; chưa khảo sát các chế độ dòng chảy hai pha phức tạp hơn như sôi bọt hoặc bay hơi trong vi kênh.
  • Các khảo sát hình học mới chỉ cố định ở 3 mức rời rạc cho 3 yếu tố hình học ($W$, $H$, $S$) mà chưa tối ưu hóa liên tục toàn dải kích thước.
  • Đồ án chủ yếu tập trung vào kết quả mô phỏng số CFD và xử lý thống kê; chưa cung cấp các số liệu đo đạc thực nghiệm đối chứng chi tiết trên hệ thống thử nghiệm vật lý của mẫu vi kênh sau khi gia công.

Hướng nghiên cứu tiếp theo được đề xuất:

  • Tiến hành lắp ráp hoàn chỉnh cụm vi kênh đã gia công tại TTL Precision vào băng thử nghiệm nhiệt - thủy lực để đo đạc thực tế nhiệt độ bề mặt, lưu lượng và sụt áp nhằm xác thực mô hình mô phỏng.
  • Mở rộng nghiên cứu sang dòng chảy hai pha, hiện tượng sôi cục bộ và khảo sát các lưu chất giải nhiệt khác như chất làm lạnh, dung dịch hạt nano (nanofluids) hoặc môi chất lạnh pha lỏng.
  • Khảo sát thêm các giải pháp bổ sung cấu trúc vi mô dạng lai phức tạp hơn (như kết hợp rãnh phụ, bề mặt nhám nhân tạo, vách ngăn hạn chế đầu vào) nhằm ức chế hiện tượng mất ổn định dòng chảy.

Giá trị tham khảo

Đồ án là nguồn tài liệu tham khảo có giá trị chuyên môn đối với sinh viên, học viên cao học và kỹ sư thuộc các ngành:

  • Công nghệ Kỹ thuật Nhiệt, Kỹ thuật Năng lượng và Kỹ thuật Điều hòa không khí: Tham khảo phương pháp tính toán truyền nhiệt liên hợp chất lỏng - vật rắn trong các thiết bị trao đổi nhiệt kích thước siêu nhỏ.
  • Kỹ thuật Cơ khí, Cơ điện tử và Thiết kế vi hệ thống: Ứng dụng quy trình dựng hình 3D trên Autodesk Inventor cho các cụm chi tiết vi cơ khí và quy trình gia công phay/cắt CNC chính xác.
  • Khoa học Dữ liệu kỹ thuật và Quản lý chất lượng: Học hỏi quy trình tích hợp phương pháp Taguchi và phân tích phương sai ANOVA trên Minitab để giải quyết các bài toán tối ưu hóa đa thông số trong kỹ thuật.

Phần nội dung đáng tham khảo nhất trong văn bản là Chương 3 và Chương 4, nơi trình bày tường minh kỹ thuật phân tách miền chất lỏng/chất rắn (Cell Zone Conditions), kỹ thuật chia lưới Mesh kiểm soát Skewness và Aspect Ratio trong ANSYS Fluent, cùng các bảng tổng hợp số liệu trích xuất từ 9 trường hợp theo ma trận trực giao.

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

1. Mục tiêu cốt lõi của đồ án là gì?
Mục tiêu của đồ án là ứng dụng phương pháp mô phỏng số CFD (ANSYS Fluent) kết hợp với các công cụ quy hoạch thực nghiệm Taguchi và phân tích phương sai (ANOVA) để xác định mức độ ảnh hưởng của các thông số hình học (chiều rộng, chiều cao, dạng tiết diện) đến đặc tính nhiệt độ, áp suất, vận tốc trong kênh tản nhiệt vi mô dạng lai, từ đó tìm ra cấu hình kích thước tối ưu.

2. Ba yếu tố hình học chính được lựa chọn khảo sát trong đồ án là gì và gồm những dạng nào?
Ba yếu tố hình học được khảo sát gồm:

  • $P_1$: Chiều rộng kênh ($W$).
  • $P_2$: Chiều cao kênh ($H$).
  • $P_3$: Dạng tiết diện rãnh/gân ($S$), bao gồm ba hình dạng: Hình chữ nhật (Rectangular - Rec), Hình thang (Trapezoidal - Trape), và Hình tam giác (Triangular - Tri).

3. Phần mềm MATLAB đóng vai trò gì trong quy trình chuẩn bị dữ liệu mô phỏng của đồ án?
MATLAB được sử dụng để lập trình hàm phân phối ngẫu nhiên randn(1,3), tạo ra sai số thực tế 1% cho các thông số điều kiện biên ban đầu trước khi nhập vào ANSYS Fluent: nhiệt độ môi trường quanh mức $300\text{ K}$, mật độ dòng nhiệt quanh mức $2.000.000\text{ W/m}^2$, và vận tốc đầu vào quanh mức $1{,}5\text{ m/s}$ (hoặc $1{,}015\text{ m/s}$).

4. Cấu hình hình học tối ưu cuối cùng mà nhóm tác giả đã xác định là gì?
Cấu hình tối ưu được lựa chọn qua phân tích là mô hình có chiều rộng kênh $W = 0{,}8\text{ mm}$, chiều cao kênh $H = 0{,}32\text{ mm}$ và dạng tiết diện rãnh là hình tam giác (Triangular).

5. Nghiên cứu đã thực hiện bước kết nối thực tế nào sau khi hoàn thành mô phỏng số?
Sau khi hoàn thành tính toán mô phỏng và xác định mô hình tối ưu, nhóm tác giả đã phối hợp gia công cơ khí chính xác mẫu kênh vi mô thực tế tại Công ty TNHH TTL Precision nhằm phục vụ cho các nghiên cứu thực nghiệm tiếp theo.

Kết luận

Đồ án tốt nghiệp của nhóm tác giả đã thiết lập thành công quy trình kết hợp chặt chẽ giữa công cụ mô phỏng động lực học chất lưu ANSYS Fluent với các phương pháp thống kê Taguchi và ANOVA trong việc phân tích các yếu tố hình học vi kênh. Kết quả nghiên cứu xác định chiều cao kênh chi phối mạnh nhất đến độ giảm nhiệt độ, trong khi chiều rộng và hình dạng tiết diện quyết định chủ yếu đến mức độ sụt áp của dòng lưu chất. Với cấu hình tối ưu được tìm ra ($W = 0{,}8\text{ mm}$, $H = 0{,}32\text{ mm}$, tiết diện tam giác) cùng việc gia công mẫu thử nghiệm thực tế, công trình cung cấp cơ sở dữ liệu và phương pháp luận đáng tin cậy cho việc phát triển các thiết bị tản nhiệt vi kênh ứng dụng trong làm mát linh kiện vi điện tử công suất cao.