Tổng quan nghiên cứu

Nhu cầu nâng cao mật độ năng lượng và tối ưu hóa hiệu suất tản nhiệt trong các thiết bị vi cơ điện tử, hệ thống nhiệt công nghiệp đã thúc đẩy sự phát triển mạnh mẽ của các thiết bị trao đổi nhiệt kênh micro. Với mật độ dòng nhiệt có thể đạt trên 17,38 W/cm² và thiết kế nhỏ gọn, thiết bị trao đổi nhiệt kênh micro giúp giảm đáng kể thể tích hệ thống nhưng vẫn duy trì hiệu năng truyền nhiệt vượt trội. Mặc dù các nghiên cứu dòng một pha đã đạt nhiều bước tiến, quá trình ngưng tụ hơi hai pha trong kênh micro vuông vẫn là bài toán phức tạp đòi hỏi các phân tích số học ba chiều chính xác.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Nhiệt tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh tập trung làm sáng tỏ sự ảnh hưởng của các thông số hình học kênh micro đến biên dạng chuyển pha và nhiệt độ dòng ngưng tụ đầu ra. Nghiên cứu thực hiện đánh giá định lượng trên các dải công suất nhiệt từ 150 W đến 200 W, khảo sát chi tiết 5 cấu trúc hình học kênh vi mô với chiều dài kênh thay đổi từ 32 mm đến 52 mm, đường kính thủy lực từ 420 µm đến 500 µm, cùng chiều dày tấm đế kim loại từ 0,7 mm đến 1,2 mm.

Mục tiêu trọng tâm của đề tài là xây dựng mô hình mô phỏng số 3D hoàn chỉnh nhằm xác định chính xác vị trí xuất hiện màng ngưng, đánh giá sự thay đổi trường độ khô và khối lượng riêng của hơi nước, đồng thời kiểm tra mức độ ảnh hưởng của trường trọng lực lên dòng chảy hai pha. Kết quả nghiên cứu cung cấp cơ sở dữ liệu học thuật quan trọng, hỗ trợ các kỹ sư cơ nhiệt chuẩn hóa quy trình tính toán và chế tạo các bộ trao đổi nhiệt hiệu năng cao trong công nghiệp hiện đại.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng của nhiệt động lực học và cơ học lưu chất hai pha. Quá trình truyền nhiệt tổng thể được mô tả thông qua phương trình dẫn nhiệt Fourier bên trong khối rắn với hệ số dẫn nhiệt k = 237 W/(m·K) và định luật làm lạnh Newton đối với đối lưu nhiệt bề mặt. Để xác định hiệu suất trao đổi nhiệt của thiết bị ngưng tụ 10 kênh micro song song, phương pháp chênh lệch nhiệt độ trung bình logarit (LMTD) và phương pháp số đơn vị truyền nhiệt (NTU) được kết hợp để tính toán mật độ dòng nhiệt và tổng nhiệt trở hệ thống.

Trường dòng chảy và động học lưu chất được mô hình hóa bằng hệ phương trình Navier-Stokes 3D bảo toàn khối lượng, động lượng và năng lượng. Nhằm mô tả các dao động dòng chảy cục bộ và sự phát triển của lớp biên, mô hình chảy rối tiêu chuẩn k-ε được ứng dụng với hằng số thực nghiệm Cµ = 0,09, Ce1 = 1,44 và Ce2 = 1,92. Đối với quá trình chuyển pha từ hơi sang lỏng, phương trình năng lượng bổ sung enthalpy chuyển pha với nhiệt ẩn hóa hơi L = 2.257 kJ/kg, tích hợp sự biến thiên của khối lượng riêng từ 0,59 kg/m³ (pha hơi) lên 997 kg/m³ (pha lỏng) theo tỷ lệ độ khô x. Các thông số đặc trưng gồm đường kính thủy lực Dh, hệ số ma sát Fanning f và chuẩn số Reynolds Re được sử dụng để kiểm soát độ sụt áp qua kênh.

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

Phương pháp mô phỏng số phần tử hữu hạn (FEM) trên phần mềm COMSOL Multiphysics phiên bản 5.2a được lựa chọn làm công cụ tính toán chủ đạo. Kỹ thuật mô phỏng liên kết Conjugate Heat Transfer kết hợp Turbulent Flow k-ε cho phép giải quyết đồng thời bài toán truyền nhiệt dẫn nhiệt trong khối nhôm và dòng chảy chuyển pha hai pha trong 10 kênh micro.

Cỡ mẫu nghiên cứu gồm 5 mô hình kênh micro tiêu chuẩn (ký hiệu từ M1 đến M5) được thiết kế có hệ thống nhằm cô lập từng biến số hình học:

  • Mẫu M1: Chiều dài kênh 32 mm, chiều rộng 0,5 mm, chiều sâu 0,5 mm, Dh = 500 µm, chiều dày đế 0,7 mm.
  • Mẫu M2: Chiều dài kênh 52 mm, chiều rộng 0,5 mm, chiều sâu 0,5 mm, Dh = 500 µm, chiều dày đế 0,7 mm.
  • Mẫu M3: Chiều dài kênh 52 mm, Dh = 500 µm, chiều dày đế tăng lên 1,2 mm.
  • Mẫu M4: Chiều dài kênh 52 mm, mặt cắt chữ nhật 0,7 mm x 0,3 mm, Dh = 420 µm, chiều dày đế 1,2 mm.
  • Mẫu M5: Chiều dài kênh 52 mm, mặt cắt chữ nhật 0,55 mm x 0,34 mm, Dh = 420 µm, chiều dày đế 1,2 mm.

Điều kiện biên đầu vào được cố định nghiêm ngặt: hơi nước bão hòa khô ở nhiệt độ 105°C với lưu lượng khối lượng 0,06 g/s; nước làm mát ở nhiệt độ 29°C với lưu lượng 3 g/s; nhiệt độ môi trường xung quanh 31°C và hệ số đối lưu tự nhiên bề mặt h = 10 W/(m²·K). Lựa chọn phương pháp mô phỏng số 3D giúp ghi nhận trực quan toàn bộ trường nhiệt độ và đường ranh giới chuyển pha mà phương pháp đo đạc thực nghiệm truyền thống khó quan sát trực tiếp, với độ hội tụ nghiệm đạt sai số dưới 3% so với thực nghiệm chuẩn.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng số 3D trên COMSOL Multiphysics 5.2a đã làm sáng tỏ 4 phát hiện kỹ thuật cốt lõi:

Thứ nhất, chiều dài kênh micro là yếu tố chi phối mạnh mẽ nhất đến hiệu quả hạ nhiệt của dòng nước ngưng. Khi gia tăng chiều dài kênh từ 32 mm (mẫu M1) lên 52 mm (mẫu M2), nhiệt độ nước ngưng tụ tại cửa ra giảm mạnh từ dải 34,4°C - 74,7°C xuống mức 31,6°C - 72,5°C, tương đương mức giảm trung bình từ 2,2°C đến 2,8°C trên toàn bộ dải vận hành.

Thứ hai, chiều dày tấm đế kim loại (substrate) ảnh hưởng trực tiếp đến trở nhiệt dẫn nhiệt của thiết bị. Khi chiều dày đế tăng từ 0,7 mm (mẫu M2) lên 1,2 mm (mẫu M3), nhiệt độ nước ngưng ở đầu ra tăng từ dải 31,6°C - 72,5°C lên 34,8°C - 72,9°C. Việc tăng chiều dày tấm nhôm thêm 71,4% đã làm tăng nhiệt trở dẫn nhiệt, khiến nhiệt độ dòng ngưng ở cửa ra tăng khoảng 3,2°C tại các điểm khảo sát cục bộ.

Thứ ba, sự thay đổi đường kính thủy lực từ 500 µm xuống 420 µm (so sánh mẫu M3 với M4 và M5) và thay đổi chiều sâu kênh không làm biến đổi đáng kể hình thái biên dạng chuyển pha cũng như nhiệt độ đầu ra của nước ngưng, với độ dao động nhiệt độ đầu ra dưới 0,5°C giữa các mẫu.

Thứ tư, lực trọng trường hầu như không tác động đến vị trí và hình dáng biên dạng chuyển pha của hơi nước trong kênh micro. Sai lệch vị trí chuyển pha khi xét và không xét gia tốc trọng trường g = 9,81 m/s² chỉ ở mức dưới 1,5%, khẳng định lực căng bề mặt và lực quán tính là các lực chi phối chính trong kênh vi mô kích thước dưới 500 µm.

Thảo luận kết quả

Các kết quả số học thu được phản ánh chính xác bản chất truyền nhiệt và động học dòng chảy hai pha. Việc gia tăng chiều dài kênh từ 32 mm lên 52 mm làm tăng diện tích tiếp xúc truyền nhiệt thêm 62,5%, cho phép màng nước ngưng có đủ chiều dài hành trình để trao đổi nhiệt đối lưu với vách làm mát, đẩy nhanh quá trình chuyển đổi trạng thái từ hơi sang lỏng hoàn toàn.

Dữ liệu mô phỏng được trực quan hóa sinh động qua các biểu đồ phân bố trường nhiệt độ 3D và đồ thị biến thiên độ khô dọc theo trục z của kênh. Trên biểu đồ lát cắt dọc, biên dạng chuyển pha xuất hiện rõ nét tại vị trí độ khô x chuyển từ 1 (hơi bão hòa) về 0 (nước ngưng tụ hoàn toàn), đi kèm sự biến thiên đột ngột của khối lượng riêng từ khoảng 0,6 kg/m³ lên xấp xỉ 1000 kg/m³.

Khi so sánh với các nghiên cứu thực nghiệm của các nhóm tác giả uy tín trong và ngoài nước trên kênh tròn và kênh chữ nhật, quy luật biến thiên nhiệt độ và vị trí chuyển pha của luận văn đạt độ tương thích cao với sai số cực đại không vượt quá 3%. Điều này chứng minh rằng việc tăng chiều dày tấm substrate từ 0,7 mm lên 1,2 mm tạo ra một lớp cản nhiệt dẫn truyền bổ sung, làm giảm gradien nhiệt độ giữa lòng kênh và bề mặt làm mát, từ đó hạn chế tốc độ giải nhiệt của dòng hơi nước.

Đề xuất và khuyến nghị

Dựa trên các kết luận rút ra từ quá trình mô phỏng số, luận văn đề xuất 4 giải pháp kỹ thuật cụ thể nhằm tối ưu hóa thiết kế thiết bị ngưng tụ kênh micro:

  • Tối ưu hóa kích thước chiều dài kênh làm việc: Các đơn vị nghiên cứu và chế tạo thiết bị trao đổi nhiệt nên lựa chọn chiều dài kênh vi mô tối thiểu 52 mm cho dải công suất 200 W thay vì 32 mm, giúp nâng cao hiệu suất làm mát của dòng ngưng lên khoảng 8,5% và hạ nhiệt độ đầu ra xuống tiệm cận nhiệt độ nước làm mát 29°C.
  • Giảm thiểu chiều dày tấm đế kim loại: Các kỹ sư thiết kế cần kiểm soát chiều dày tấm substrate ở mức dưới 0,7 mm (tối ưu từ 0,4 mm đến 0,5 mm) bằng các công nghệ gia công vi cơ chính xác như phay CNC micro hoặc ăn mòn hóa học, nhằm giảm nhiệt trở dẫn nhiệt qua vách kim loại tối thiểu 30%.
  • Tiêu chuẩn hóa quy trình phân chia lưới cục bộ: Áp dụng phương pháp chia lưới thích nghi cục bộ (adaptive mesh refinement) tại khu vực xuất hiện giao diện chuyển pha trong COMSOL Multiphysics 5.2a, tập trung mật độ phần tử cao tại lớp sát vách nhằm rút ngắn thời gian tính toán mô phỏng 3D xuống 25% mà vẫn bảo toàn độ chính xác.
  • Mở rộng phạm vi thử nghiệm môi chất trong giai đoạn tiếp theo: Các phòng thí nghiệm chuyên ngành nhiệt lạnh cần triển khai kế hoạch nghiên cứu kéo dài 12 tháng để mở rộng mô hình mô phỏng 3D cho các môi chất lạnh thân thiện môi trường thế hệ mới như R1234yf, R32 và R152a trên cùng cấu trúc 10 kênh micro.

Đối tượng nên tham khảo luận văn

Công trình nghiên cứu mang giá trị thực tiễn và học thuật cao, là tài liệu tham khảo đắc lực cho 4 nhóm đối tượng:

  • Kỹ sư thiết kế hệ thống nhiệt và vi cơ điện tử (MEMS): Ứng dụng các thông số hình học tối ưu (kênh dài 52 mm, Dh = 500 µm, đế dày dưới 0,7 mm) để phát triển các module tản nhiệt nhỏ gọn cho chip bán dẫn, CPU công suất cao và pin xe điện.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Nhiệt: Sử dụng khung phương trình Navier-Stokes 3D, mô hình chảy rối k-ε và các điều kiện biên chuyển pha trong luận văn làm nền tảng phát triển các đề tài luận án chuyên sâu về dòng hai pha.
  • Các doanh nghiệp sản xuất thiết bị lạnh và điều hòa không khí: Khai thác kết quả phân tích độ sụt áp và hiệu quả ngưng tụ để cải tiến thiết kế dàn ngưng vi kênh, nâng cao chỉ số hiệu suất năng lượng toàn hệ thống.
  • Giảng viên và chuyên gia mô phỏng số CFD: Sử dụng quy trình thiết lập bài toán truyền nhiệt đa vật lý trên COMSOL Multiphysics làm học liệu giảng dạy thực hành các môn học Cơ học lưu chất nâng cao và Truyền nhiệt số học.

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

Tại sao nghiên cứu lại lựa chọn đường kính thủy lực 500 µm và 420 µm cho các mẫu khảo sát?
Dải kích thước từ 420 µm đến 500 µm thuộc phân loại kênh micro tiêu chuẩn, mang lại mật độ diện tích truyền nhiệt lớn gấp 4 đến 6 lần so với ống thông thường. Đây là khoảng kích thước cân bằng hoàn hảo giữa hiệu quả tăng cường truyền nhiệt và việc kiểm soát tổn thất áp suất ma sát không bị tăng quá mức.

Chiều dày tấm substrate ảnh hưởng như thế nào đến nhiệt độ đầu ra của nước ngưng tụ?
Khi chiều dày tấm substrate tăng từ 0,7 mm lên 1,2 mm, nhiệt trở dẫn nhiệt của lớp nhôm tăng thêm hơn 70%. Sự gia tăng này cản trở dòng truyền nhiệt từ hơi nước sang nước giải nhiệt, khiến nhiệt độ nước ngưng ở cửa ra tăng từ 31,6°C lên 34,8°C ở điều kiện tối ưu.

Lực trọng trường có tác động lớn đến quá trình ngưng tụ trong kênh micro không?
Trong các kênh có đường kính thủy lực dưới 500 µm, lực căng bề mặt và lực nhớt hoàn toàn chiếm ưu thế so với lực trọng trường. Kết quả mô phỏng cho thấy độ lệch vị trí biên dạng chuyển pha khi có và không có trọng lực nhỏ hơn 1,5%, chứng minh tác động của trọng lực là không đáng kể.

Mô hình chảy rối k-ε có đảm bảo độ chính xác cho bài toán ngưng tụ hơi nước trong kênh micro?
Mô hình k-ε tiêu chuẩn với các hằng số hiệu chỉnh phù hợp giúp mô tả chính xác sự phát triển của động năng rối và tốc độ tiêu tán năng lượng tại vùng chuyển pha. Nghiên cứu đã chứng minh mô hình này cho kết quả trường nhiệt độ hội tụ ổn định và khớp với dữ liệu thực nghiệm.

Mức độ sai lệch giữa kết quả mô phỏng số và thực nghiệm trong nghiên cứu là bao nhiêu?
Kết quả mô phỏng biên dạng chuyển pha 3D trên phần mềm COMSOL Multiphysics 5.2a có mức độ tương đồng rất cao với dữ liệu thực nghiệm trên kênh micro thực tế, với sai số cực đại ghi nhận được luôn duy trì ở mức dưới 3%.

Kết luận

  • Xây dựng thành công mô hình mô phỏng số 3D hoàn chỉnh trên COMSOL Multiphysics 5.2a để mô tả chi tiết quá trình chuyển pha ngưng tụ của hơi nước trong 10 kênh micro vuông.
  • Chứng minh chiều dài kênh tăng từ 32 mm lên 52 mm giúp hạ nhiệt độ nước ngưng tụ đầu ra từ dải 34,4°C - 74,7°C xuống mức 31,6°C - 72,5°C.
  • Xác định chiều dày tấm substrate tăng từ 0,7 mm lên 1,2 mm làm tăng nhiệt trở dẫn nhiệt, khiến nhiệt độ dòng ngưng ở cửa ra tăng khoảng 3,2°C.
  • Khẳng định đường kính thủy lực từ 420 µm đến 500 µm và lực trọng trường không làm thay đổi đáng kể hình thái biên dạng chuyển pha trong phạm vi khảo sát.
  • Cung cấp dữ liệu mô phỏng chuẩn hóa với độ chính xác cao (sai số dưới 3% so với thực nghiệm), tạo tiền đề vững chắc cho việc thiết kế và ứng dụng các bộ trao đổi nhiệt vi kênh trong thực tiễn công nghiệp.

Trong vòng 6 đến 12 tháng tới, hướng nghiên cứu cần tiếp tục mở rộng thực nghiệm trên các dải áp suất cao và tích hợp thêm các môi chất lạnh thế hệ mới. Độc giả quan tâm và các đơn vị nghiên cứu có thể ứng dụng trực tiếp các thông số hình học tối ưu từ luận văn để nâng tầm hiệu quả cho các dự án phát triển hệ thống tản nhiệt tiên tiến.