Tổng quan nghiên cứu

Trong bối cảnh công nghệ vi xử lý phát triển vượt bậc, các bộ vi xử lý máy tính cao cấp khi vận hành ở chế độ hiệu năng cao hoặc ép xung thường tỏa ra lượng nhiệt khổng lồ, khiến nhiệt độ bề mặt chip chạm ngưỡng giới hạn 85°C. Đồng thời, trong các ngành công nghiệp nặng như tua bin khí với nhiệt độ buồng đốt lên tới 1200°C đến 1600°C, các giải pháp tản nhiệt truyền thống bằng quạt khí đối lưu tự nhiên không còn đáp ứng được yêu cầu giải nhiệt nhanh trên diện tích bề mặt hẹp. Dòng phun tia tác động trực giao (Jet Impingement) nổi lên như một giải pháp làm mát có hiệu năng truyền nhiệt vượt trội nhờ khả năng loại bỏ thông lượng nhiệt cục bộ cực lớn.

Tuy nhiên, cơ chế thủy khí động lực học của dòng tia va chạm vuông góc là một bài toán phức tạp, bao gồm sự kết hợp đồng thời của dòng chảy trong ống, vùng phun tự do, vùng va đập điểm dừng và dòng tia bám tường. Tại Việt Nam, các nghiên cứu chuyên sâu về cấu trúc trường vận tốc và phân bố truyền nhiệt của dòng phun ngập chu trình kín còn khá khiêm tốn do hạn chế về trang thiết bị thí nghiệm quang học hiện đại.

Đề tài nghiên cứu hướng đến mục tiêu xây dựng cơ sở lý thuyết toàn diện và thực hiện mô phỏng số động lực học chất lưu trên phần mềm chuyên dụng ANSYS Fluent. Phạm vi nghiên cứu tập trung khảo sát cấu hình dòng phun ngập tác động vuông góc lên bề mặt phẳng ở hai chế độ chảy rối với hệ số Reynolds Re = 23.000 và Re = 70.000, tương ứng với hai khoảng cách vòi phun không thứ nguyên H/D = 2 và H/D = 6. Kết quả nghiên cứu có ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng giúp cải thiện hệ số truyền nhiệt đối lưu từ 25% đến 40% so với phương pháp thổi khí thông thường, mở ra giải pháp tối ưu hóa thiết kế bộ tản nhiệt khối nước cho thiết bị điện tử công suất lớn.

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 dựa trên hệ thống các định luật bảo toàn cơ bản của cơ học chất lưu và truyền nhiệt đối lưu:

  • Hệ phương trình vi phân liên tục và phương trình Navier-Stokes: Mô tả quá trình bảo toàn khối lượng và động lượng của dòng lưu chất Newton không nén được trong không gian ba chiều.
  • Phương trình năng lượng nhiệt đối lưu: Kết hợp với định luật Newton về làm mát bề mặt, xác định mối liên hệ giữa thông lượng nhiệt và độ chênh lệch nhiệt độ giữa tường chắn và dòng tia.
  • Các mô hình mô tả chảy rối: Luận văn tổng quan các mô hình rối hai phương trình phổ biến gồm mô hình k-epsilon tiêu chuẩn, mô hình k-epsilon Realizable, mô hình RNG k-epsilon và mô hình k-omega của Wilcox. Trong đó, các biến số chính bao gồm động năng rối k, tốc độ tiêu tán năng lượng rối epsilon và tần số tiêu tán rối riêng phần omega.

Khung khái niệm trọng tâm bao gồm: hệ số Reynolds đặc trưng cho chế độ chảy rối của lưu chất; hệ số Nusselt đại diện cho tỷ số giữa truyền nhiệt đối lưu và dẫn nhiệt thuần túy; hệ số Prandtl mô tả tỷ số giữa độ khuếch tán động lượng và độ khuếch tán nhiệt độ; vùng lõi tiềm năng nơi vận tốc trục dòng tia chưa bị suy giảm; và điểm dừng nơi vận tốc dòng triệt tiêu về 0 tạo áp suất tĩnh cực đại.

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

  • Nguồn dữ liệu kiểm chứng: Dữ liệu thực nghiệm chuẩn quốc tế từ công trình của Cooper và các cộng sự công bố năm 1993 với đường kính vòi phun D = 101,6 mm, bề dày thành ống s = 2 mm, kết hợp dữ liệu truyền nhiệt của Baughn và Shimizu đối với bài toán dòng va chạm trực giao.
  • Phương pháp phân tích: Sử dụng phương pháp thể tích hữu hạn (Finite Volume Method - FVM) tích hợp trong bộ giải ANSYS Fluent. Lý do lựa chọn FVM là vì phương pháp này đảm bảo tính bảo toàn khối lượng, động lượng và năng lượng một cách tuyệt đối trên từng ô thể tích điều khiển rời rạc, có khả năng thích ứng linh hoạt với các gradient áp suất cao và hiện tượng phân tách dòng phức tạp tại bề mặt tường chắn.
  • Kỹ thuật chọn mẫu và chia lưới: Bài toán được quy về mô hình hình học 2D đối xứng trục nhằm tiết kiệm tài nguyên tính toán nhưng vẫn bảo toàn đầy đủ bản chất vật lý 3D. Miền tính toán được rời rạc hóa với mật độ hơn 200.000 phần tử lưới có cấu trúc. Lưới được làm mịn đặc biệt tại vùng sát vách với giá trị y+ nhỏ hơn 1 để nắm bắt chính xác cấu trúc phân lớp nhớt và lớp biên nhiệt. Chiều dài ống cấp được kéo dài theo công thức chuẩn Le/D = 4,4 nhân với Re lũy thừa 1/6 nhằm đảm bảo dòng chảy đã phát triển hoàn toàn trước khi thoát ra khỏi miệng vòi. Toàn bộ chu trình mô phỏng và phân tích dữ liệu được thực hiện hoàn chỉnh trong năm 2014.

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

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

  • Phân bố trường vận tốc và cấu trúc dòng: Dòng chảy khi thoát ra khỏi miệng vòi duy trì vùng lõi tiềm năng dạng hình nón thu hẹp. Vận tốc dọc trục đạt cực đại tại tâm vòi x/D = 0 và giảm nhanh về 0 khi chạm điểm dừng trên tường chắn. Sau cú va đập, dòng đổi hướng 90 độ, tăng tốc cục bộ tạo thành dòng tia bám tường chảy tràn theo phương bán kính r/D.
  • Ảnh hưởng của hệ số Reynolds: Khi tăng hệ số Reynolds từ Re = 23.000 lên Re = 70.000 (tăng hơn 204%), động năng rối k tại vùng biên cắt tiếp giáp môi trường xung quanh tăng gấp hơn 3,2 lần. Sự gia tăng xung lượng này thúc đẩy hệ số truyền nhiệt đối lưu trung bình trên toàn bộ bề mặt tấm phẳng tăng từ 65% đến 82%.
  • Ảnh hưởng của khoảng cách vòi phun H/D: Tại khoảng cách H/D = 6, hệ số Nusselt tại điểm dừng đạt giá trị tuyệt đối cao nhất do chiều dài lõi tiềm năng kết thúc ngay trước bề mặt va đập, cho phép các xoáy rối quy mô lớn xâm nhập trực tiếp và phá vỡ lớp biên nhiệt.
  • Hiện tượng cực đại thứ cấp của hệ số Nusselt: Ở khoảng cách hẹp H/D = 2, đường cong phân bố hệ số Nusselt xuất hiện một đỉnh cực đại thứ hai tại vị trí bán kính r/D xấp xỉ 1,8 đến 2,0 do sự chuyển tiếp đột ngột từ chế độ chảy tầng sang chảy rối trong lớp biên bám tường.

Thảo luận kết quả

Cơ chế vật lý của hiện tượng truyền nhiệt được sáng tỏ thông qua sự tương tác giữa động học dòng xoáy và lớp biên nhiệt. Khi H/D = 6, các cấu trúc xoáy hình khuyên hình thành từ tầng trượt của dòng phun tự do phát triển mạnh mẽ, làm tăng cường độ xáo trộn giữa lưu chất mát và bề mặt nóng. Ngược lại, khi H/D = 2, lớp biên trên vách ban đầu phát triển ở chế độ tầng rất mỏng, sau đó mất ổn định do gradient áp suất nghịch và chuyển pha sang chảy rối, tạo ra đỉnh thứ cấp của hệ số Nusselt.

So sánh với các công bố thực nghiệm của Cooper và cộng sự, kết quả mô phỏng vận tốc dọc trục và vận tốc tiếp tuyến dọc tường chắn có độ sai lệch tương đối dưới 7%, khẳng định độ tin cậy của mô hình số. Trong báo cáo kỹ thuật, dữ liệu được trực quan hóa sinh động thông qua biểu đồ đường dòng thể hiện tâm xoáy thứ cấp, biểu đồ vector trường vận tốc 2D và đồ thị phân bố hệ số Nusselt theo khoảng cách bán kính không thứ nguyên r/D. Kết quả này chứng minh rằng việc bố trí khoảng cách vòi phun H/D = 6 mang lại hiệu quả tản nhiệt tổng thể đồng đều nhất cho các tấm làm mát CPU.

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

  • Tối ưu hóa khoảng cách vòi phun trong thiết kế tản nhiệt: Các kỹ sư nhiệt tại các doanh nghiệp chế tạo phần cứng máy tính nên thiết lập khoảng cách vòi phun cố định ở tỷ lệ H/D = 6 cho các hệ thống tản nhiệt nước chu trình kín, giúp nâng cao hiệu suất truyền nhiệt tổng thể thêm 15% đến 20%, triển khai áp dụng ngay trong quý 1 đến quý 2 năm tới.
  • Ứng dụng cấu hình vòi phun dạng côn thu: Bộ phận R&D sản xuất thiết bị làm mát công nghiệp cần thay thế miệng vòi phẳng truyền thống bằng vòi phun côn thu hẹp ở đầu ra nhằm gia tăng vận tốc tia và phân bố động năng rối đồng đều hơn, hướng tới mục tiêu giảm nhiệt độ hoạt động của chip xuống dưới 70°C.
  • Phát triển mô hình mảng đa vòi phun: Các nhóm nghiên cứu cơ khí cần mở rộng bài toán từ vòi phun đơn lẻ sang mảng vòi phun đa điểm với khoảng cách bước giữa các vòi S/D từ 4 đến 8, nhằm xử lý bài toán tản nhiệt diện rộng cho các trạm máy chủ và hệ thống pin xe điện trong lộ trình 6 đến 12 tháng.
  • Chuẩn hóa quy trình chia lưới thích ứng trong mô phỏng CFD: Các viện nghiên cứu và phòng thí nghiệm tính toán cần áp dụng tiêu chuẩn kiểm soát lớp biên nhớt y+ nhỏ hơn 1 kết hợp với mô hình rối k-omega SST, giúp giảm sai số dự báo nhiệt động học xuống dưới mức 5%.

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

  • Kỹ sư thiết kế nhiệt và phần cứng điện tử: Khai thác chi tiết các thông số hình học vòi phun, lưu lượng dòng chảy và quy luật phân bố hệ số Nusselt để thiết kế các khối block làm mát chất lỏng hiệu năng cao cho vi xử lý máy trạm và GPU.
  • Học viên cao học và nghiên cứu sinh ngành Cơ kỹ thuật - Kỹ thuật Nhiệt: Sử dụng luận văn như một tài liệu tham khảo chuẩn mực về phương pháp thiết lập bài toán 2D đối xứng trục, kỹ thuật chia lưới phi cấu trúc và phương pháp giải phương trình vi phân phi tuyến trên phần mềm ANSYS Fluent.
  • Chuyên gia công nghệ tua bin và hàng không vũ trụ: Vận dụng cơ chế dòng va đập trực giao để cải tiến hệ thống làm mát bên trong cánh tua bin khí hoạt động trong môi trường nhiệt độ khắc nghiệt từ 1200°C trở lên.
  • Kỹ sư vận hành trong ngành sấy dệt nhuộm và luyện kim: Ứng dụng nguyên lý dòng tia tốc độ cao để tối ưu hóa dây chuyền sấy khô vải, giấy tự động hoặc hệ thống phun sương làm nguội phôi thép cán nóng, giúp tiết kiệm 10% đến 15% năng lượng tiêu thụ.

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

  • Tại sao dòng phun va chạm trực giao lại đạt hiệu quả truyền nhiệt cao hơn dòng chảy song song? Dòng phun vuông góc ép các cột lưu chất vận tốc cao đập trực tiếp vào bề mặt nóng, làm mỏng tối đa lớp biên thủy động và lớp biên nhiệt tại vùng điểm dừng. Quá trình này tạo ra gradient nhiệt độ cực lớn tại sát vách, thúc đẩy hệ số truyền nhiệt đối lưu tăng gấp 2 đến 3 lần so với dòng chảy dọc theo bề mặt thông thường.

  • Vì sao hệ số Nusselt tại điểm dừng đạt giá trị cao nhất khi khoảng cách vòi phun đạt H/D = 6? Ở khoảng cách H/D = 6, chiều dài của vùng lõi tiềm năng vừa vặn kết thúc ngay sát bề mặt va đập. Khi đó, các xoáy rối quy mô lớn hình thành từ lớp biên cắt xung quanh tia phun thâm nhập trực tiếp vào dòng chính, tạo ra cường độ dao động mạch động vận tốc cực đại phá hủy lớp biên nhiệt.

  • Hiện tượng đỉnh cực đại thứ cấp của hệ số Nusselt ở H/D = 2 hình thành như thế nào? Khi vòi phun đặt rất sát tường ở khoảng cách H/D = 2, lớp biên bám tường ban đầu phát triển theo chế độ chảy tầng mỏng. Khi dòng tiếp tục lan ra ngoài đến vị trí bán kính r/D từ 1,8 đến 2,0, sự gia tăng bất ổn định dòng chảy kích hoạt quá trình chuyển tiếp sang chảy rối, làm tăng đột biến tốc độ truyền nhiệt cục bộ.

  • Việc kéo dài đường ống theo công thức Le/D = 4,4 nhân với Re lũy thừa 1/6 có ý nghĩa gì trong mô phỏng? Công thức này xác định chiều dài tối thiểu của đoạn ống dẫn để biên dạng vận tốc của dòng chảy bên trong đạt trạng thái phát triển hoàn toàn trước khi phun ra môi trường. Điều này loại bỏ sai số đầu vào của bài toán mô phỏng, đảm bảo dữ liệu đầu ra tương thích hoàn hảo với các điều kiện thực nghiệm tiêu chuẩn.

  • Mô hình chảy rối nào cho kết quả tính toán tối ưu nhất trong bài toán dòng va chạm? Mô hình k-epsilon Realizable và k-omega là hai lựa chọn tối ưu nhờ khả năng mô tả chính xác sự tiêu tán năng lượng rối và ứng suất nhớt tại vùng dòng va đập chịu gradient áp suất lớn, giúp rút ngắn thời gian hội tụ và kiểm soát sai số tính toán dưới 8%.

Kết luận

  • Luận văn đã hệ thống hóa hoàn chỉnh cơ sở lý thuyết nhiệt - thủy khí động lực học của dòng phun tia ngập tác động trực giao lên bề mặt phẳng trong hệ thống làm mát chu trình kín.
  • Mô hình số FVM trên ANSYS Fluent mô phỏng thành công trường vận tốc và phân bố truyền nhiệt tại hai chế độ Re = 23.000 và Re = 70.000 với độ chính xác cao.
  • Nghiên cứu đã xác lập quy luật biến thiên của hệ số Nusselt, làm rõ cơ chế xuất hiện đỉnh truyền nhiệt thứ cấp tại H/D = 2 và cực đại điểm dừng tại H/D = 6.
  • Đóng góp chính của đề tài là cung cấp giải pháp thiết kế định lượng cho các hệ thống tản nhiệt chất lỏng vi mạch, khống chế nhiệt độ chip dưới ngưỡng 85°C một cách an toàn.
  • Hướng phát triển trong 12 đến 24 tháng tới là mở rộng mô hình sang bài toán dòng phun 3D đa pha kết hợp chất lưu nano. Các kỹ sư và nhà nghiên cứu quan tâm nên ứng dụng ngay quy trình mô phỏng này để tối ưu hóa sản phẩm tản nhiệt công nghiệp.