Tổng quan nghiên cứu

Nhu cầu tản nhiệt mật độ cao trong các thiết bị vi điện tử và hệ thống công nghiệp hiện đại đang đối mặt với thách thức giải phóng mật độ dòng nhiệt vượt ngưỡng 100 W/cm² đến 120 W/cm² trong một không gian giới hạn. Các thiết bị trao đổi nhiệt truyền thống thường gặp nhược điểm về kích thước cồng kềnh, khối lượng lưu chất nạp lớn và hệ số truyền nhiệt đối lưu chưa tối ưu. Để giải quyết triệt để vấn đề này, công nghệ kênh vi mô (microchannel) đã mở ra bước đột phá mới nhờ tỷ lệ diện tích bề mặt trên thể tích vượt trội, giúp tăng cường cường độ trao đổi nhiệt và giảm thiểu không gian lắp đặt.

Luận văn tập trung nghiên cứu định lượng ảnh hưởng của hình dáng hình học đến quá trình bay hơi của dòng chảy hai pha lỏng - hơi trong kênh vi mô bằng việc kết hợp đồng thời phương pháp mô phỏng số và kiểm chứng thực nghiệm. Mục tiêu cụ thể là phân tích, so sánh cơ chế chuyển pha, sự hình thành mầm bọt hơi, nhiệt độ đầu ra và độ khô giữa hai dạng hình học: kênh hình chữ nhật tiêu chuẩn và kênh hình chữ nhật bo tròn đáy.

Công trình được thực hiện tại Bộ môn 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 Thành phố Hồ Chí Minh, hoàn thành vào tháng 10 năm 2017. Các thí nghiệm được thiết lập trong điều kiện kiểm soát nghiêm ngặt với công suất nguồn nhiệt cấp 132 W, lưu lượng khối lượng nước cấp 0,3 g/s, nhiệt độ nước đầu vào 60°C và nhiệt độ môi trường 30°C. Kết quả của nghiên cứu mang ý nghĩa thực tiễn quan trọng, đóng góp nền tảng dữ liệu giúp cắt giảm từ 15% đến 25% kích thước bộ trao đổi nhiệt và nâng cao hiệu suất truyền nhiệt tổng thể cho các hệ thống nhiệt lạnh thế hệ mớ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 vững chắc của ba trụ cột lý thuyết nhiệt động lực học và thủy lực học vi mô:

Thứ nhất, định luật thứ nhất nhiệt động lực học và phương trình dẫn nhiệt Fourier tổng quát. Truyền nhiệt dẫn qua vật liệu đế nhôm có hệ số dẫn nhiệt k = 160 W/mK, khối lượng riêng 2700 kg/m³ và nhiệt dung riêng 900 J/kgK được mô hình hóa qua phương trình vi phân truyền nhiệt kết hợp giữa dẫn nhiệt và đối lưu cưỡng bức.

Thứ hai, lý thuyết thủy động lực học dòng chảy hai pha và mô hình dòng chảy rối k - ε mức Reynolds thấp (Low Re k-ε). Mô hình giải quyết hệ phương trình Navier-Stokes, phương trình liên tục và phương trình năng lượng để mô tả chính xác sự phân bố vận tốc, tổn thất áp suất và sự tiêu tán năng lượng động lực học trong không gian vi kênh.

Thứ ba, lý thuyết mao dẫn và cơ chế truyền nhiệt khi sôi bọt trong kênh vi mô. Nghiên cứu khai thác sâu các khái niệm cốt lõi:

  • Lực căng bề mặt và hiệu ứng mao dẫn: Sức căng mặt ngoài chi phối áp suất mao dẫn tại mặt phân chia pha lỏng - hơi.
  • Góc tiếp xúc và góc dính ướt: Xác định theo phương trình Young để đánh giá độ sụt áp tại mặt lõm chất lỏng dính ướt.
  • Bán kính mầm hơi tới hạn: Được tính toán theo độ quá nhiệt của bề mặt, trong đó khi độ quá nhiệt tăng từ 5°C lên 25°C thì bán kính mầm hơi giảm nhanh từ 6,7 µm xuống còn 1,3 µm.
  • Chu kỳ sản sinh bọt hơi: Xác định khoảng thời gian phát triển và tách ly bọt hơi tiêu chuẩn đạt 0,05 giây đối với nước ở áp suất khí quyển.

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

Nghiên cứu kết hợp chặt chẽ giữa tính toán mô phỏng số phần tử hữu hạn 3D và đo đạc thực nghiệm đối chứng:

Về nguồn dữ liệu và công cụ phân tích, tác giả ứng dụng phần mềm COMSOL Multiphysics phiên bản 5.3 với các module liên kết Đa vật lý (Conjugate Heat Transfer, Phase Change và Low Re k-ε). Dữ liệu thực nghiệm được thu thập trực tiếp qua hệ thống đo tự động gồm bộ thu nhận xử lý tín hiệu Yokogawa MX100, cân điện tử phân tích độ chính xác cao, các cảm biến nhiệt độ tiếp xúc và bộ nguồn gia nhiệt ổn định.

Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu chế tạo và khảo sát 2 bộ trao đổi nhiệt vi kênh bằng nhôm nguyên khối có kích thước tổng thể đồng nhất là 140 mm x 19,5 mm, chiều dày tấm đế 1,2 mm và được ốp kín bằng tấm cách nhiệt PMMA kích thước 152 mm x 31,5 mm:

  • Mẫu vi kênh chữ nhật: Gồm 25 kênh thẳng song song, kích thước tiết diện mỗi kênh là chiều rộng 0,3 mm x chiều sâu 0,5 mm, khoảng cách giữa các kênh là 0,5 mm, diện tích mặt cắt ngang đạt 0,15 mm².
  • Mẫu vi kênh chữ nhật bo tròn: Gồm 20 kênh song song, kích thước tiết diện là chiều rộng 0,07 mm x chiều sâu 0,5 mm kết hợp cung tròn bo đáy bán kính R = 0,25 mm, diện tích mặt cắt ngang được thiết kế tương đương 0,15 mm².

Lý do lựa chọn phương pháp phân tích này là nhằm cô lập hoàn toàn yếu tố diện tích tiết diện dòng chảy (cùng A = 0,15 mm²), giúp đánh giá chính xác tác động thuần túy của hình dáng hình học biên dạng kênh đến tốc độ bay hơi trong toàn bộ chu trình quan sát kéo dài 150 giây, đảm bảo sai số giữa lý thuyết và thực tế được kiểm soát dưới 9%.

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 và thực nghiệm đã làm sáng tỏ các phát hiện kỹ thuật mang tính bước ngoặt:

Thứ nhất, biên dạng bo tròn giúp kích hoạt quá trình chuyển pha sớm hơn đáng kể. Tại kênh hình chữ nhật được bo tròn, những bong bóng hơi đầu tiên xuất hiện và bắt đầu quá trình sôi tại giây thứ 17, nhanh hơn 5,56% so với giây thứ 18 ở kênh hình chữ nhật tiêu chuẩn.

Thứ hai, độ khô dòng hơi đầu ra ở kênh bo tròn vượt trội hơn. Tại thời điểm quan sát chuẩn ở giây thứ 150, chất lượng hơi và độ khô của kênh bo tròn đạt giá trị cao hơn từ 8% đến 12% so với kênh chữ nhật phẳng, cho thấy tốc độ bay hơi và khả năng chuyển hóa pha lỏng thành pha hơi diễn ra triệt để hơn.

Thứ ba, sự phân bố nhiệt độ và biên dạng nhiệt đồng đều hơn. Kênh hình chữ nhật bo tròn duy trì nhiệt độ đầu ra ổn định hơn, giảm hiện tượng quá nhiệt cục bộ trên bề mặt vách nhôm và hạn chế tối đa sự mất ổn định dòng chảy do tắc nghẽn bọt khí dạng khối dài (slug flow).

Thứ tư, độ chính xác và tương thích cao giữa mô phỏng và thực nghiệm. Dữ liệu ghi nhận từ hệ thống đo Yokogawa MX100 cho thấy sự trùng khớp rất cao với kết quả tính toán từ COMSOL Multiphysics 5.3, với sai số nhiệt độ trung bình đầu ra luôn duy trì dưới 6,5% trong suốt dải công suất gia nhiệt 132 W.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự vượt trội của kênh bo tròn nằm ở sự thay đổi phân bố lực căng bề mặt và động lực học mầm bọt. Ở kênh chữ nhật vuông góc, chất lỏng có xu hướng bị bẫy lại tại các góc nhọn do lực mao dẫn dính ướt, làm cản trở quá trình đối lưu màng và làm chậm tốc độ tách ly của bọt hơi. Ngược lại, bán kính cong R = 0,25 mm ở đáy kênh bo tròn tạo điều kiện cho các mầm hơi dễ dàng dịch chuyển lên phía mặt thoáng với vận tốc nâng bọt tăng từ 0,17 m/s lên 0,35 m/s, giúp bề mặt gia nhiệt liên tục được làm ướt lại và kích hoạt các chu kỳ sôi tiếp theo nhanh chóng.

So sánh với các nghiên cứu tiền nhiệm, kết quả này củng cố phát hiện của Sui về việc tăng cường truyền nhiệt qua các bề mặt vi kênh biến dạng, đồng thời bổ sung bằng chứng thực nghiệm rõ nét hơn so với nghiên cứu của Balasubramanian vốn chỉ khảo sát kênh thẳng mở rộng với độ giảm áp 30%. Độ chênh lệch nhiệt độ và lưu lượng trong đề tài này cũng hoàn toàn tương thích với mức sai số 9% đã được Dang công bố trên các hệ thống trao đổi nhiệt vi mô dòng ngược chiều.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua các đồ thị phân bố nhiệt độ theo thời gian từ giây thứ 0 đến giây thứ 150, biểu đồ so sánh độ khô đầu ra tại các mốc thời gian chuyển pha, và bảng tổng hợp đối chiếu thông số áp suất - nhiệt độ giữa hai mô hình hình học. Bản đồ trường nhiệt 3D trích xuất từ phần mềm mô phỏng cho thấy rõ ranh giới chuyển pha lỏng - hơi di chuyển nhanh hơn về phía cửa vào của kênh bo tròn, minh chứng rõ ràng cho hiệu suất gia nhiệt vượt bậc.

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

Nhằm chuyển giao và ứng dụng hiệu quả các kết quả nghiên cứu vào thực tiễn sản xuất, bốn nhóm giải pháp hành động cụ thể được đề xuất như sau:

Ứng dụng công nghệ phay vi cơ khí CNC hoặc gia công cắt dây tia lửa điện (Wire-EDM) độ chính xác cao để tạo hình góc lượn đáy R = 0,25 mm cho các tấm tản nhiệt vi kênh nhôm. Giải pháp này hướng tới mục tiêu nâng cao 15% hiệu suất giải nhiệt cho các khối làm mát CPU/GPU công suất lớn trong giai đoạn 2024 - 2025, do các kỹ sư R&D phần cứng điện tử chủ trì thực hiện.

Chuẩn hóa quy trình thiết kế mô phỏng đa vật lý trên nền tảng COMSOL Multiphysics kết hợp module dòng chảy Low Re k-ε tại các phòng thí nghiệm nhiệt lạnh. Mục tiêu là cắt giảm 30% thời gian thử nghiệm và tiết kiệm 40% chi phí chế tạo mẫu thử nghiệm cơ khí trong vòng 6 tháng triển khai tại các viện nghiên cứu chuyên ngành.

Tích hợp hệ thống kiểm soát lưu lượng dòng cấp siêu nhỏ tự động với độ chính xác đạt 0,01 g/s trên các dây chuyền bơm nhiệt vi mô. Giải pháp giúp duy trì độ khô đầu ra của thiết bị bay hơi đạt trên 90% một cách ổn định, hướng đến áp dụng cho các doanh nghiệp sản xuất điều hòa không khí chính xác trước năm 2026.

Mở rộng nghiên cứu thực nghiệm với các môi chất lạnh thân thiện với môi trường như khí CO2 tự nhiên hoặc dung dịch hydrofluoroether HFE-7100 ở dải mật độ dòng nhiệt từ 50 kW/m² đến 150 kW/m². Lộ trình thực hiện dự kiến kéo dài 12 đến 18 tháng, do các nhóm nghiên cứu khoa học cơ khí động lực phối hợp cùng doanh nghiệp năng lượng triển khai.

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

Nội dung và dữ liệu thực nghiệm trong luận văn mang lại giá trị học thuật và ứng dụng thiết thực cho bốn nhóm đối tượng chính:

Kỹ sư thiết kế nhiệt và hệ thống làm mát điện tử: Nắm bắt các thông số hình học chuẩn xác của kênh vi mô (kích thước 0,3 mm x 0,5 mm và 0,07 mm x 0,5 mm bo tròn) để ứng dụng trực tiếp vào việc tối ưu hóa khối tản nhiệt nước cho trung tâm dữ liệu, bộ biến tần năng lượng mặt trời và pin xe điện thế hệ mới.

Học viên cao học, nghiên cứu sinh và giảng viên ngành Kỹ thuật Nhiệt - Năng lượng: Sử dụng luận văn như một tài liệu tham khảo bài bản về phương pháp luận kết hợp giữa giải tích số học trên COMSOL Multiphysics 5.3 và kỹ thuật đo lường thực nghiệm dòng hai pha với thiết bị thu thập Yokogawa MX100.

Doanh nghiệp sản xuất thiết bị trao đổi nhiệt công nghiệp: Khai thác các kết quả so sánh độ khô và biên dạng nhiệt để cải tiến cấu trúc dàn bay hơi vi kênh, giúp giảm 20% lượng môi chất nạp vào hệ thống và thu gọn đáng kể diện tích chiếm chỗ của cụm tản nhiệt.

Chuyên gia phát triển phần mềm mô phỏng dòng chảy động lực học (CFD): Tiếp cận bộ dữ liệu đo đạc thực nghiệm độ chính xác cao về thời điểm hình thành bọt hơi (giây thứ 17 và 18) để hiệu chuẩn, tinh chỉnh các mô hình toán học dự đoán quá trình sôi vi mô và chuyển pha lỏng - hơi.

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

Tại sao kênh vi mô bo tròn lại bắt đầu quá trình sôi và bay hơi nhanh hơn kênh chữ nhật tiêu chuẩn? Biên dạng đáy bo tròn bán kính R = 0,25 mm giúp loại bỏ các điểm ứ đọng chất lỏng tại góc nhọn, phân bố đều lực căng bề mặt và giảm lực cản mao dẫn. Nhờ đó, các mầm hơi hình thành thuận lợi và bứt phá khỏi vách truyền nhiệt nhôm ngay tại giây thứ 17 (nhanh hơn kênh chữ nhật 1 giây), tạo chu kỳ sản sinh bọt hơi liên tục chỉ mất khoảng 0,05 giây.

Mức độ tin cậy của phần mềm COMSOL Multiphysics 5.3 trong bài toán mô phỏng dòng hai pha này như thế nào? COMSOL Multiphysics 5.3 thể hiện độ tin cậy vượt trội khi giải đồng thời các phương trình liên tục, Navier-Stokes và truyền nhiệt đa pha. Sai số nhiệt độ đầu ra giữa mô phỏng và dữ liệu đo thực tế từ hệ thống Yokogawa MX100 luôn duy trì ở mức dưới 6,5%, hoàn toàn đáp ứng các tiêu chuẩn khắt khe trong nghiên cứu học thuật quốc tế.

Tiêu chí nào được sử dụng để lựa chọn kích thước hình học giữa hai mẫu kênh vi mô khảo sát? Cả hai mẫu đều được thiết kế để duy trì diện tích mặt cắt ngang dòng chảy hoàn toàn bằng nhau ở mức 0,15 mm² (mẫu chữ nhật 0,3 mm x 0,5 mm với 25 kênh và mẫu bo tròn 0,07 mm x 0,5 mm với 20 kênh). Tiêu chí này giúp đồng nhất vận tốc và lưu lượng dòng nạp, cô lập tuyệt đối biến số hình dáng hình học để đánh giá chính xác tác động của góc bo.

Công suất nhiệt 132 W và lưu lượng 0,3 g/s phản ánh điều kiện làm việc thực tế nào? Bộ thông số này mô phỏng chính xác chế độ tải nhiệt của các vi xử lý bán dẫn hiệu năng cao hoặc mô-đun điện tử công suất có mật độ dòng nhiệt cục bộ từ 50 W/cm² đến 120 W/cm². Lưu lượng cấp 0,3 g/s đảm bảo dòng đối lưu hai pha hoạt động ổn định trong dải nhiệt độ nước vào 60°C mà không gây ra hiện tượng cháy khô bề mặt.

Kết quả của đề tài có thể mở rộng cho các loại môi chất làm mát khác ngoài nước cất không? Hoàn toàn có thể mở rộng. Khung lý thuyết và mô hình mô phỏng dòng hai pha trong luận văn là nền tảng chuẩn để áp dụng cho các môi chất lạnh như R134a, CO2 hoặc chất lỏng điện môi HFE. Khi chuyển đổi môi chất, kỹ sư chỉ cần cập nhật các thông số đặc trưng gồm ẩn nhiệt hóa hơi, sức căng mặt ngoài và độ nhớt động lực học tương ứng.

Kết luận

Những đóng góp cốt lõi và định hướng phát triển từ công trình nghiên cứu bao gồm:

  • Chứng minh rõ nét kênh vi mô bo tròn đáy (R = 0,25 mm) vượt trội hơn kênh chữ nhật tiêu chuẩn (0,3 mm x 0,5 mm) về tốc độ bắt đầu chuyển pha (sớm hơn tại giây thứ 17) và độ khô hơi đầu ra.
  • Xây dựng thành công quy trình mô phỏng 3D dòng hai pha chuyển pha trên phần mềm COMSOL Multiphysics 5.3 với sai số đối chiếu thực nghiệm được kiểm soát dưới 6,5%.
  • Xác lập bộ thông số vận hành nhiệt - thủy lực tối ưu trong không gian kênh vi mô với công suất 132 W, lưu lượng 0,3 g/s, nhiệt độ nước cấp 60°C trên nền vật liệu nhôm dẫn nhiệt 160 W/mK.
  • Cung cấp cơ sở lý thuyết chi tiết về động học mầm hơi vi mô với đường kính tới hạn từ 1,3 µm đến 6,7 µm và cơ chế tác động của lực mao dẫn đến quá trình sôi bọt.
  • Đưa ra giải pháp thực tiễn giúp thu nhỏ từ 15% đến 20% kích thước các bộ trao đổi nhiệt nhỏ gọn ứng dụng trong tản nhiệt vi điện tử và công nghệ điều hòa không khí.

Trong kế hoạch 12 tháng tới, nhóm nghiên cứu định hướng mở rộng thử nghiệm trên môi chất CO2 siêu tới hạn ở dải áp suất làm việc trên 80 bar. Hãy liên hệ bộ phận học thuật hoặc truy cập thư viện trường để khai thác trọn vẹn tập dữ liệu thiết kế vi kênh tối ưu cho các dự án nhiệt lạnh chuyên sâu của bạn.