Tổng quan nghiên cứu

Nhu cầu sử dụng năng lượng cho các hệ thống điều hòa không khí và làm lạnh công nghiệp hiện chiếm khoảng 40% tổng năng lượng tiêu thụ trong các công trình kiến trúc thương mại và dân dụng trên phạm vi toàn cầu. Trong hệ thống này, thiết bị trao đổi nhiệt đóng vai trò then chốt, quyết định trực tiếp đến hiệu suất năng lượng và kích thước vật lý của toàn bộ chu trình nhiệt động. Việc nâng cao hiệu quả truyền nhiệt khi sôi bên trong ống thông qua việc sử dụng ống có bề mặt tăng cường như ống vi cánh (micro-fin) đã được chứng minh giúp gia tăng hệ số truyền nhiệt từ 150% đến 250% so với ống trơn truyền thống. Tuy nhiên, sự xuất hiện của cấu trúc vi cánh bên trong lòng ống làm cho cơ chế động học dòng chảy hai pha và hiện tượng truyền nhiệt sôi trở nên vô cùng phức tạp, khiến các công thức bán thực nghiệm cổ điển áp dụng cho ống trơn thường gặp sai số lớn từ 25% đến 40%.

Mục tiêu cốt lõi của đề tài là xây dựng và phát triển phương trình tương quan dự đoán chính xác hệ số trao đổi nhiệt khi sôi dòng bên trong ống micro-fin nằm ngang, áp dụng cho hai môi chất lạnh phổ biến là R410A và R32 dựa trên mô hình cấu trúc dòng chảy thông qua phương pháp mô phỏng số. Công trình nghiên cứu được triển khai tại Trường Cơ khí thuộc Đại học Bách Khoa Hà Nội và hoàn thành vào tháng 5 năm 2022, tập trung khảo sát dải lưu lượng khối từ 200 đến 320 kg/(m²·s) cùng tỷ lượng thể tích biến thiên từ 0,2 đến 0,85. Kết quả nghiên cứu mang lại ý nghĩa kinh tế và kỹ thuật sâu sắc, hỗ trợ giảm từ 15% đến 20% khối lượng kim loại màu trong gia công ống trao đổi nhiệt, đồng thời góp phần nâng cao hệ số hiệu quả năng lượng (COP) của các thiết bị nhiệt lạnh lên khoảng 8% đến 12%.

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

Khung lý thuyết áp dụng

Khung lý thuyết của luận văn dựa trên nền tảng của nhiệt động lực học kỹ thuật và cơ học chất lưu hai pha, kết hợp giữa hai cơ chế truyền nhiệt sôi cơ bản: sôi mầm (nucleate boiling) và sôi đối lưu cưỡng bức (forced convective boiling). Trong ống micro-fin, hai cơ chế này diễn ra đồng thời nhưng có tỷ trọng đóng góp thay đổi tùy thuộc vào điều kiện vận hành. Sôi mầm chịu sự chi phối mạnh mẽ bởi mật độ dòng nhiệt và áp suất bão hòa, trong khi sôi đối lưu phụ thuộc chủ yếu vào lưu lượng khối của dòng môi chất và tỷ lượng thể tích của pha hơi.

Nghiên cứu kế thừa và phát triển mô hình cộng hưởng của Chen kết hợp với mô hình phân loại chế độ dòng chảy hiện đại của Wojtan, Ursenbacher và Thome. Khung phân tích tập trung vào các khái niệm trọng tâm bao gồm: hệ số trao đổi nhiệt hai pha toàn phần, hệ số ức chế sôi mầm, hệ số tăng cường trao đổi nhiệt đối lưu, tỷ lượng thể tích pha hơi và đường kính thủy lực tương đương của ống vi cánh. Bên cạnh đó, lý thuyết động lực học chất lưu tính toán nhiều pha được thiết lập dựa trên hệ phương trình liên tục, phương trình động lượng Navier-Stokes kết hợp mô hình lực căng bề mặt màng liên tục và mô hình dòng rối tiêu chuẩn nhằm mô tả chính xác sự phân bố của màng chất lỏng trên bề mặt vi cánh.

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

Nguồn dữ liệu nhiệt vật lý của các môi chất lạnh R410A và R32 được trích xuất đồng bộ từ phần mềm cơ sở dữ liệu nhiệt động chuyên dụng NIST REFPROP 9. Cỡ mẫu nghiên cứu bao gồm 36 trường hợp mô phỏng số 3D độc lập, được thiết lập dựa trên ma trận thực nghiệm kết hợp giữa hai loại môi chất lạnh, ba cấp lưu lượng khối gồm 200, 260 và 320 kg/(m²·s), cùng sáu giá trị tỷ lượng thể tích pha hơi trải đều từ 0,2; 0,4; 0,6; 0,7; 0,8 đến 0,85.

Phương pháp chọn mẫu thông số vận hành được thực hiện theo nguyên tắc quy hoạch thực nghiệm có hệ thống, phản ánh chính xác các chế độ làm việc điển hình của thiết bị bay hơi từ tải cục bộ đến đầy tải trong thực tế kỹ thuật. Phương pháp phân tích chủ đạo là mô phỏng động lực học dòng chảy (CFD) sử dụng mô hình thể tích chất lưu (Volume of Fluid - VOF) kết hợp xử lý lưới phi cấu trúc mật độ cao tại vùng rãnh cánh và đỉnh cánh vi mô có chiều cao từ 0,15 đến 0,25 mm và góc xoắn cánh 18 độ. Lý do lựa chọn mô hình VOF trong không gian ba chiều là khả năng bắt bám chính xác ranh giới pha lỏng - hơi và tái hiện chuyển động xoáy của màng chất lỏng dọc theo rãnh cánh, điều mà các phương pháp đo đạc quang học truyền thống trong thực nghiệm rất khó ghi nhận đầy đủ. Toàn bộ chu trình nghiên cứu, xây dựng mô hình lưới, giải số và phân tích dữ liệu được thực hiện liên tục trong khoảng thời gian từ tháng 10 năm 2021 đến tháng 5 năm 2022.

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

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

Nghiên cứu đã ghi nhận những phát hiện quan trọng về động học dòng chảy và đặc tính truyền nhiệt của môi chất sôi trong ống micro-fin:

  • Quy luật chuyển dịch chế độ dòng chảy: Sự hiện diện của cánh vi mô tạo ra hiệu ứng dẫn dòng và lực mao dẫn mạnh mẽ, khiến cấu trúc dòng chảy chuyển đổi từ dòng phân tầng lượn sóng sang dòng màng hình khuyên sớm hơn khoảng 15% về giá trị tỷ lượng thể tích so với ống trơn cùng đường kính thủy lực. Khi lưu lượng khối tăng từ 200 lên 320 kg/(m²·s), vùng xuất hiện dòng màng hình khuyên mở rộng đáng kể, giúp bề mặt ống luôn duy trì được trạng thái làm ướt tốt.
  • Hiệu quả duy trì màng chất lỏng và hạn chế hiện tượng khô màng: Lực ly tâm sinh ra từ góc xoắn của cánh vi mô liên tục đẩy pha lỏng lên đỉnh ống, ngăn ngừa sự xuất hiện của các điểm khô cục bộ sớm. Môi chất R32 cho thấy khả năng duy trì màng lỏng ổn định hơn R410A, giúp hệ số trao đổi nhiệt của R32 cao hơn từ 18% đến 22% ở cùng điều kiện lưu lượng khối 260 kg/(m²·s) và cùng mật độ dòng nhiệt.
  • Quy luật biến thiên hệ số truyền nhiệt theo tỷ lượng thể tích: Trong dải tỷ lượng thể tích từ 0,2 đến 0,6, hệ số trao đổi nhiệt hai pha tăng mạnh với mức tăng trung bình từ 35% đến 50% do vận tốc dòng hơi tăng làm tăng cường đối lưu cưỡng bức. Hệ số này đạt giá trị cực đại trong khoảng tỷ lượng thể tích từ 0,75 đến 0,85 trước khi có xu hướng suy giảm nhẹ khi màng chất lỏng bắt đầu mỏng dần.
  • Độ chính xác của phương trình tương quan mới: Phương trình tương quan được phát triển dựa trên mô hình dòng chảy có khả năng dự đoán hệ số trao đổi nhiệt hai pha với độ lệch chuẩn trung bình chỉ khoảng 11,5%, trong khi sai số tối đa không vượt quá 18% đối với cả hai môi chất R410A và R32, vượt trội hoàn toàn so với các công thức truyền thống có sai lệch lên tới 35%.

Thảo luận kết quả

Sự vượt trội về đặc tính trao đổi nhiệt của R32 so với R410A xuất phát từ các đặc tính nhiệt vật lý thuận lợi, bao gồm nhiệt ẩn hóa hơi cao hơn khoảng 38% và độ nhớt pha lỏng thấp hơn, tạo điều kiện cho các bọt hơi hình thành nhanh chóng và thoát ly dễ dàng khỏi rãnh cánh. Khi so sánh với các nghiên cứu thực nghiệm nhiệt lạnh công bố gần đây, kết quả từ mô hình CFD hoàn toàn tương thích và phản ánh đúng bản chất vật lý của quá trình sôi màng hỗn hợp.

Về mặt biểu diễn dữ liệu, toàn bộ quy luật biến thiên có thể được thể hiện trực quan thông qua các hình thức sau:

  • Bản đồ chế độ dòng chảy dạng đồ thị hai chiều với trục hoành là tỷ lượng thể tích và trục tung là lưu lượng khối, phân định rõ ràng các ranh giới giữa dòng nút, dòng phân tầng lượn sóng và dòng màng hình khuyên.
  • Biểu đồ đường phân bố hệ số truyền nhiệt theo tỷ lượng thể tích, cho phép so sánh trực quan các đường cong đặc tính giữa ống micro-fin và ống trơn ở ba cấp lưu lượng khối 200, 260 và 320 kg/(m²·s).
  • Bảng đối chiếu sai số tương đối giữa số liệu dự đoán từ phương trình tương quan mới và số liệu thực nghiệm, minh chứng cho độ tin cậy cao của mô hình tính toán.

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

Dựa trên các kết quả đạt được, luận văn đưa ra bốn nhóm khuyến nghị có tính ứng dụng cao cho ngành kỹ thuật nhiệt:

  • Chuẩn hóa thông số hình học của ống vi cánh trong sản xuất: Các kỹ sư thiết kế tại các nhà máy chế tạo thiết bị nhiệt lạnh cần tối ưu hóa góc xoắn cánh trong khoảng từ 15 độ đến 18 độ và duy trì chiều cao cánh từ 0,18 đến 0,22 mm. Giải pháp này giúp tăng hệ số trao đổi nhiệt thêm từ 10% đến 15% trong khi chỉ làm tăng tổn thất áp suất dưới mức 8%, dự kiến hoàn thành trong lộ trình 6 tháng.
  • Tích hợp phương trình tương quan vào phần mềm tính toán thiết kế: Các phòng nghiên cứu và phát triển (R&D) cần nhúng trực tiếp thuật toán tương quan mới vào các bộ công cụ tính toán nhiệt của doanh nghiệp nhằm giảm thiểu 40% thời gian thử nghiệm nguyên mẫu thực tế, triển khai trong vòng 9 tháng tới.
  • Thúc đẩy chuyển đổi môi chất lạnh sang R32: Các doanh nghiệp sản xuất thiết bị điều hòa không khí cần đẩy nhanh lộ trình thay thế môi chất R410A bằng R32 cho các dòng sản phẩm thương mại nhằm giảm 67% chỉ số tiềm năng gây nóng lên toàn cầu (GWP) và nâng cao hiệu suất chu trình thêm khoảng 5% đến 7% trong giai đoạn từ 1 đến 2 năm.
  • Mở rộng cơ sở dữ liệu mô phỏng cho các môi chất lạnh thế hệ mới: Các viện nghiên cứu và trường đại học chuyên ngành kỹ thuật nhiệt cần tiếp tục ứng dụng mô hình số để khảo sát các môi chất hydrofluoroolefin (HFO) thân thiện môi trường trên dải lưu lượng mở rộng từ 100 đến 500 kg/(m²·s), hoàn thành trong thời gian 12 tháng.

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

Công trình nghiên cứu mang lại giá trị học thuật và thực tiễn phong phú cho các nhóm độc giả chuyên sâu:

  • Kỹ sư nghiên cứu và phát triển thiết bị nhiệt lạnh: Cung cấp công cụ tính toán chính xác để tối ưu hóa kích thước bề mặt dàn bay hơi, giúp giảm từ 15% đến 20% không gian lắp đặt và hạ giá thành vật tư sản xuất.
  • Nhà nghiên cứu và học viên cao học chuyên ngành Kỹ thuật nhiệt: Cung cấp tài liệu tham khảo chuẩn mực về phương pháp số mô phỏng dòng chảy hai pha, kỹ thuật xây dựng lưới vi kênh và quy trình xử lý mô hình thể tích chất lưu (VOF).
  • Doanh nghiệp sản xuất và gia công ống đồng trao đổi nhiệt: Tham chiếu các thông số hình học vi mô tối ưu để cải tiến khuôn kéo ống và dây chuyền gia công cơ khí chính xác.
  • Giảng viên các trường đại học khối kỹ thuật: Bổ sung nguồn học liệu giá trị cho các học phần Thiết bị trao đổi nhiệt, Truyền nhiệt nâng cao và Ứng dụng CFD trong nhiệt lạnh.

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

  • Ống micro-fin giúp tăng cường hiệu quả trao đổi nhiệt dựa trên cơ chế nào? Cánh vi mô bên trong ống làm tăng diện tích tiếp xúc truyền nhiệt thực tế lên từ 1,3 đến 1,8 lần so với ống trơn. Đồng thời, cấu trúc rãnh xoắn tạo ra lực xoáy ly tâm đẩy màng chất lỏng lên làm ướt chu vi ống và tạo nhiều điểm mầm sôi nhân tạo, giúp hệ số truyền nhiệt tổng thể tăng từ 150% đến 250%.

  • Tại sao môi chất lạnh R32 lại cho hiệu quả trao đổi nhiệt cao hơn R410A? R32 sở hữu các đặc tính nhiệt động vượt trội với nhiệt ẩn hóa hơi cao hơn khoảng 38% và hệ số dẫn nhiệt pha lỏng lớn hơn R410A. Điều này giúp tăng cường độ sôi mầm bên trong các khe cánh và tạo vận tốc dòng hơi cao hơn, thúc đẩy quá trình đối lưu cưỡng bức diễn ra mạnh mẽ hơn.

  • Mô hình thể tích chất lưu (VOF) đem lại lợi thế gì trong bài toán mô phỏng dòng hai pha? Mô hình VOF giải phương trình liên tục cho tỷ lượng thể tích của từng pha, cho phép theo dõi chính xác sự biến dạng liên tục của giao diện phân cách giữa pha lỏng và pha hơi. Nhờ đó, mô hình tái hiện chân thực sự phân bố màng chất lỏng dưới tác động của trọng lực và lực căng bề mặt.

  • Phương trình tương quan trong luận văn có thể áp dụng cho các điều kiện nào? Phương trình được xây dựng tối ưu cho các môi chất R410A và R32 chảy sôi bên trong ống micro-fin nằm ngang, hoạt động trong dải lưu lượng khối từ 200 đến 320 kg/(m²·s), tỷ lượng thể tích từ 0,2 đến 0,85 và dải nhiệt độ bão hòa từ 5 đến 15 độ C với độ chính xác đạt trên 88%.

  • Việc gia tăng mật độ cánh micro-fin có luôn mang lại hiệu quả kinh tế cao hơn không? Không hoàn toàn, bởi vì việc tăng quá nhiều cánh vi mô sẽ làm tăng sức cản ma sát thủy lực, dẫn đến tổn thất áp suất tăng vọt từ 30% đến 60%. Điều này buộc máy nén phải tiêu thụ nhiều điện năng hơn, do đó cần lựa chọn mật độ cánh tối ưu để cân bằng giữa truyền nhiệt và tổn thất áp suất.

Kết luận

  • Luận văn đã mô phỏng thành công động học dòng chảy hai pha của môi chất R410A và R32 sôi bên trong ống micro-fin nằm ngang bằng phương pháp CFD với mô hình VOF 3D.
  • Xác lập được bản đồ chế độ dòng chảy chi tiết, chỉ ra sự chuyển dịch sớm sang chế độ dòng màng hình khuyên ở tỷ lượng thể tích trên 0,4 nhờ tác động dẫn hướng của cánh vi mô.
  • Phát triển thành công phương trình tương quan cải tiến dự đoán hệ số trao đổi nhiệt hai pha với độ chính xác vượt trội, giảm độ lệch chuẩn trung bình xuống dưới 12%.
  • Đóng góp luận cứ khoa học quan trọng cho việc chuyển đổi môi chất lạnh thân thiện môi trường R32 và tối ưu hóa kết cấu thiết bị bay hơi trong ngành công nghiệp lạnh.
  • Trong giai đoạn từ 6 đến 12 tháng tới, hướng phát triển tiếp theo cần mở rộng mô phỏng cho các môi chất hỗn hợp mới có chỉ số GWP cực thấp và tích hợp dữ liệu vào các tiêu chuẩn thiết kế công nghiệp.

Quý độc giả, kỹ sư và nhà nghiên cứu quan tâm đến mô hình tính toán nhiệt và giải pháp tối ưu hóa thiết bị trao đổi nhiệt có thể liên hệ trực tiếp với nhóm tác giả tại Khoa Năng lượng Nhiệt, Trường Cơ khí - Đại học Bách Khoa Hà Nội để cùng hợp tác, trao đổi chuyên môn và ứng dụng thực tiễn công nghệ.