Tổng quan nghiên cứu

Trong bối cảnh các nguồn nhiên liệu hóa thạch truyền thống như than đá và dầu mỏ đang dần cạn kiệt, việc thu hồi các nguồn nhiệt thải công nghiệp có nhiệt thế thấp từ 80°C đến 150°C đóng vai trò đặc biệt quan trọng trong chiến lược tiết kiệm năng lượng toàn cầu. Tại Việt Nam, sự gia tăng nhanh chóng của các khu công nghiệp luyện kim, hóa chất, nhà máy nhiệt điện cùng các phụ phẩm nông nghiệp như trấu, bã mía đang thải ra môi trường một lượng nhiệt năng dư thừa rất lớn. Tuy nhiên, việc ứng dụng máy lạnh hấp thụ sử dụng nguồn nhiệt thải này còn gặp nhiều rào cản kỹ thuật do công tác tính toán nhiệt động truyền thống bằng phương pháp tra bảng và đồ thị entanpi - nồng độ (đồ thị i-c) tiêu tốn nhiều thời gian và có độ chuẩn xác chưa cao.

Mục tiêu trọng tâm của nghiên cứu là ứng dụng phần mềm Matlab để xây dựng chương trình số hóa tự động tính toán chu trình nhiệt động thực của hệ thống máy lạnh hấp thụ amoniac - nước (NH3 - H2O) một cấp hồi nhiệt. Đề tài tập trung giải quyết bài toán thiết kế với công suất lạnh định mức 175,8 kW (tương đương 50 TR), khảo sát dải nhiệt độ nguồn nhiệt cấp từ 120°C đến 165°C, nhiệt độ bay hơi từ 0°C đến -20°C và nhiệt độ nước giải nhiệt từ 22°C đến 45°C. Ý nghĩa của nghiên cứu thể hiện qua việc tối ưu hóa hệ số sử dụng nhiệt COP (đạt từ 0,22 đến 0,36) và kiểm soát bội số tuần hoàn dung dịch f (giảm xuống mức 2,21 đến 3,65 ở điều kiện tiêu chuẩn), giúp giảm hơn 80% lượng điện tiêu thụ so với máy lạnh nén hơi, đồng thời loại trừ hoàn toàn việc phát thải các chất làm suy giảm tầng ozone và gây hiệu ứng nhà kính.

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 3 trụ cột lý thuyết nhiệt động lực học và truyền nhiệt chuyên sâu:

  • Lý thuyết chu trình máy lạnh hấp thụ một cấp hồi nhiệt: Mô hình hóa toàn bộ quá trình biến đổi pha và trao đổi năng lượng giữa các thiết bị chính bao gồm bình phát sinh, bình hấp thụ, bình ngưng tụ, dàn bay hơi, bình ngưng hồi lưu và thiết bị trao đổi nhiệt dung dịch.
  • Phương trình trạng thái và nhiệt vật lý dung dịch hai cấu tử NH3 - H2O: Ứng dụng hệ thống công thức giải tích giải mã mối quan hệ tương quan giữa áp suất bão hòa p, nhiệt độ sôi T, nồng độ dung dịch c, nhiệt dung riêng cp và hệ số dẫn nhiệt lambda.
  • Lý thuyết cân bằng nhiệt và bảo toàn vật chất: Xác định cân bằng năng lượng từng thiết bị với hệ số hiệu suất trao đổi nhiệt đạt 95% và hệ số dự phòng tổn thất lạnh 5%.

Bên cạnh đó, nghiên cứu làm rõ 4 khái niệm nhiệt kỹ thuật cốt lõi:

  • Hệ số sử dụng nhiệt COP: Tỷ số giữa năng suất lạnh nhận được tại dàn bay hơi và tổng năng lượng cấp cho bình phát sinh cùng điện năng tiêu hao của bơm dung dịch.
  • Bội số tuần hoàn f: Tỷ lệ khối lượng dung dịch đậm đặc tuần hoàn trên một đơn vị khối lượng môi chất hơi amoniac sinh ra.
  • Nồng độ dung dịch đậm đặc cs và nồng độ dung dịch loãng cw: Đại lượng đặc trưng cho khả năng hấp thụ và phát sinh khí NH3 trong chu trình.
  • Quá trình tinh luyện và hồi lưu: Giải pháp công nghệ tách nước để nâng cao nồng độ hơi NH3 trước khi vào bình ngưng đạt độ tinh khiết tối thiểu 99,8%.

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

Nghiên cứu thiết lập tập dữ liệu mô phỏng gồm 48 điểm tính toán nhiệt động cho bài toán biến thiên nhiệt độ nước làm mát và 30 kịch bản mô phỏng cho bài toán nhiệt độ bay hơi. Phương pháp chọn mẫu là chọn mẫu phân tầng có chủ đích (purposive stratified sampling) dựa trên các khoảng biến thiên thông số thực tế của công nghiệp lạnh (nhiệt độ nước làm mát 24°C đến 45°C, nhiệt độ nguồn nhiệt cấp 125°C đến 140°C).

Lý do lựa chọn giải thuật số trị trên nền tảng Matlab là khả năng giải quyết các hệ phương trình đại số phi tuyến bậc cao mô tả hàm trạng thái nhiệt động dung dịch NH3 - H2O một cách tự động, khắc phục triệt để sai số trực quan từ 5% đến 10% khi tra cứu thủ công trên đồ thị i-c và t-c. Toàn bộ nguồn dữ liệu kiểm chứng được đối soát trực tiếp với các số liệu chuẩn của ngành kỹ thuật nhiệt lạnh quốc tế, hoàn thiện quy trình phân tích và thử nghiệm mô phỏng trong chu kỳ 12 tháng.

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

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

Nghiên cứu đã lượng hóa chi tiết các mối quan hệ nhiệt động then chốt trong hệ thống máy lạnh hấp thụ NH3 - H2O:

  • Ảnh hưởng của nhiệt độ nước làm mát đến bội số tuần hoàn: Tại mức nhiệt độ nguồn nhiệt cấp 134°C, khi nhiệt độ nước làm mát tăng từ 24°C lên 42°C, bội số tuần hoàn f tăng vọt từ 2,21 lên 20,14 (tương đương mức tăng hơn 811%). Ở nhiệt độ nước làm mát 44°C và nguồn nhiệt 130°C, giá trị f đạt đỉnh 59,83, làm gia tăng đột biến công tiêu hao của bơm tuần hoàn.
  • Ảnh hưởng của nhiệt độ bay hơi đến hệ số sử dụng nhiệt COP: Khi nhiệt độ bay hơi te giảm từ -1°C xuống -10°C ở nhiệt độ nguồn cấp 140°C, hệ số COP giảm liên tục từ 0,36 xuống 0,06 (suy giảm 83,3%). Với nguồn nhiệt cấp 134°C, hệ số COP duy trì ổn định ở mức 0,30 tại te = -2°C và 0,32 tại te = -1°C.
  • Giới hạn tối ưu của nguồn nhiệt cấp: Khi tăng nhiệt độ nguồn cấp từ 125°C lên 140°C, hệ số COP được cải thiện từ 15% đến 25% trong toàn dải nhiệt độ bay hơi. Tuy nhiên, nếu nhiệt độ nguồn cấp vượt ngưỡng 150°C, lượng hơi nước bốc theo NH3 tăng mạnh, làm quá tải tháp tinh luyện và giảm 8% năng suất dàn bay hơi.
  • Độ chính xác của phần mềm mô phỏng: Kết quả so sánh giữa chương trình Matlab và phương pháp tra đồ thị i-c cho thấy entanpi lỏng ngưng tụ tại áp suất ngưng 20 bar đạt 585 kJ/kg và entanpi hơi sau bốc hơi tại áp suất bay hơi 4 bar đạt 1605 kJ/kg, sai số tương đối giữa hai phương pháp chỉ dưới 1,5%.

Thảo luận kết quả

Nguyên nhân vật lý của sự thay đổi hiệu suất xuất phát từ hiện tượng dịch chuyển cân bằng pha lỏng - hơi của dung dịch NH3 - H2O. Khi nhiệt độ nước làm mát tăng, áp suất ngưng tụ trong bình ngưng tăng lên mức 20 bar, làm nhiệt độ bình hấp thụ tăng theo và làm suy giảm khả năng hòa tan khí amoniac vào dung dịch loãng. Hệ quả là hiệu số nồng độ (cs - cw) bị thu hẹp đáng kể, buộc hệ thống phải tăng lưu lượng dung dịch tuần hoàn f để đáp ứng cùng một công suất lạnh 175,8 kW, từ đó làm giảm hiệu quả năng lượng toàn hệ thống.

Về mặt biểu diễn trực quan, các tương quan này được thể hiện rõ ràng thông qua hệ thống biểu đồ đường đa biến (Line Chart) biểu diễn biến thiên của f và COP theo từng đường đẳng nhiệt nguồn cấp (125°C, 130°C, 134°C, 140°C), đi kèm bảng ma trận 10 điểm nút trạng thái nhiệt động (áp suất, nhiệt độ, nồng độ, entanpi). Kết quả nghiên cứu hoàn toàn tương thích với các báo cáo thực nghiệm của ngành nhiệt lạnh thế giới, đồng thời chứng minh phần mềm xây dựng giúp rút ngắn 90% thời gian tính toán chu trình thực so với phương pháp truyền thống.

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

Dựa trên kết quả nghiên cứu mô phỏng, 4 giải pháp cụ thể được khuyến nghị triển khai:

  • Phát triển giao diện người dùng trực quan cho phần mềm tính toán: Xây dựng bộ công cụ GUI trên nền tảng số nhằm tự động hóa hoàn toàn việc thiết lập thông số chu trình, hướng tới mục tiêu giảm 95% thời gian thiết kế sơ bộ máy lạnh hấp thụ, hoàn thành trong vòng 6 tháng bởi nhóm nghiên cứu và kỹ sư lập trình nhiệt lạnh.
  • Tối ưu hóa hệ thống tháp giải nhiệt công nghiệp: Lắp đặt cụm tháp giải nhiệt hiệu suất cao nhằm duy trì nhiệt độ nước làm mát tuần hoàn luôn dưới 32°C, đảm bảo bội số tuần hoàn f không vượt quá 6,5 và duy trì hệ số COP trên mức 0,30, triển khai trong 3 đến 6 tháng tại các đơn vị chế tạo thiết bị lạnh.
  • Cải tiến kết cấu tháp tinh luyện và bộ hồi lưu: Thiết kế tháp chưng cất đa tầng tích hợp bộ trao đổi nhiệt dung dịch nhằm nâng nồng độ hơi NH3 vào bình ngưng đạt ngưỡng 99,8%, giúp tiết kiệm 15% nhiệt lượng cấp cho bình phát sinh, áp dụng lộ trình 9 tháng bởi các kỹ sư công nghệ tại các nhà máy có nguồn thải nhiệt lớn.
  • Chế tạo mô hình thử nghiệm pilot quy mô 50 TR: Triển khai lắp đặt hệ thống máy lạnh hấp thụ NH3 - H2O thực nghiệm tại nhà máy công nghiệp để đối soát số liệu thực tế, kiểm soát sai số mô phỏng dưới 2%, thực hiện trong giai đoạn 12 đến 18 tháng bởi các viện nghiên cứu phối hợp cùng doanh nghiệp cơ khí chế tạo.

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

Nội dung luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế và tư vấn cơ điện lạnh: Tiếp cận thuật toán số hóa để tính toán kích thước bình ngưng, tháp tinh luyện và bình phát sinh; ứng dụng thiết kế hệ thống cấp đông bảo quản thủy hải sản công suất 175,8 kW sử dụng nhiệt khói thải lò hơi.
  • Giảng viên và sinh viên ngành Kỹ thuật Nhiệt - Năng lượng: Sử dụng cơ sở lý thuyết nhiệt động dung dịch hai cấu tử và phương pháp mô phỏng trên Matlab làm tài liệu giảng dạy, nghiên cứu chuyên đề hoặc đồ án tốt nghiệp chuyên ngành nhiệt lạnh.
  • Cán bộ quản lý năng lượng tại các nhà máy công nghiệp: Đánh giá tiềm năng kinh tế và kỹ thuật trong việc thu hồi nguồn nhiệt thải từ 80°C đến 150°C để tạo nước đá hoặc làm mát xưởng sản xuất; xây dựng báo cáo khả thi đầu tư dự án tiết kiệm năng lượng.
  • Doanh nghiệp sản xuất và tích hợp thiết bị công nghệ xanh: Khai thác công cụ phần mềm để tối ưu hóa vật liệu chế tạo thiết bị trao đổi nhiệt, phát triển dòng máy chiller hấp thụ thay thế máy nén lạnh Freon nhằm đạt các chứng chỉ công trình xanh.

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

Máy lạnh hấp thụ NH3 - H2O khác gì so với máy lạnh nén hơi truyền thống? Hệ thống hấp thụ thay thế máy nén tiêu tốn nhiều điện năng bằng cụm thiết bị gồm bình hấp thụ, bơm dung dịch và bình phát sinh sử dụng nhiệt thừa từ 80°C đến 150°C. Cơ cấu này chỉ sử dụng một động cơ bơm nhỏ, giúp giảm hơn 80% lượng điện tiêu thụ và vận hành êm ái, rất phù hợp cho những khu vực thiếu điện hoặc giá điện cao.

Vì sao nghiên cứu lựa chọn cặp môi chất NH3 - H2O thay vì H2O - LiBr? Cặp dung dịch H2O - LiBr sử dụng nước làm môi chất bay hơi nên chỉ có thể tạo ra nhiệt độ làm lạnh trên 0°C phục vụ điều hòa không khí. Ngược lại, cặp NH3 - H2O cho phép nhiệt độ bay hơi hạ sâu xuống từ -2°C đến -20°C, là lựa chọn bắt buộc cho các hệ thống cấp đông thủy hải sản và sản xuất nước đá công nghiệp.

Tại sao phải lắp đặt tháp chưng cất và bộ ngưng tụ hồi lưu trong hệ thống? Khi dung dịch NH3 - H2O sôi trong bình phát sinh, một lượng nhỏ hơi nước sẽ bị cuốn theo amoniac. Nếu không loại bỏ hoàn toàn, lượng nước này sẽ đóng băng tại van tiết lưu và làm tắc nghẽn đường ống. Tháp chưng cất và bộ hồi lưu giúp duy trì nồng độ hơi amoniac tinh khiết đạt 99,8%, bảo vệ hệ thống vận hành liên tục.

Nhiệt độ nước làm mát ảnh hưởng như thế nào đến hiệu suất hệ thống? Khi nhiệt độ nước làm mát tăng từ 24°C lên 42°C, áp suất ngưng tụ tăng cao làm dung dịch loãng hấp thụ amoniac kém hơn, dẫn đến bội số tuần hoàn f tăng vọt từ 2,21 lên 20,14 ở nguồn nhiệt 134°C. Điều này làm tăng điện năng tiêu thụ của bơm và làm suy giảm nghiêm trọng hệ số COP của chu trình.

Việc ứng dụng phần mềm Matlab mang lại ưu thế gì so với phương pháp tra đồ thị cổ điển? Phương pháp tra đồ thị entanpi - nồng độ (i-c) thủ công thường gặp sai số từ 5% đến 10% do thao tác gióng đường kẻ trực quan và tốn nhiều thời gian khi thay đổi thông số đầu vào. Phần mềm Matlab giải chính xác các hệ phương trình trạng thái phi tuyến trong vài giây, hỗ trợ tối ưu hóa đa biến nhanh chóng và chuẩn xác.

Kết luận

  • Đã xây dựng hoàn chỉnh cơ sở lý thuyết nhiệt động và hệ thống phương trình giải tích xác định chính xác các thông số trạng thái của máy lạnh hấp thụ NH3 - H2O một cấp hồi nhiệt.
  • Ứng dụng thành công ngôn ngữ lập trình Matlab để phát triển phần mềm tính toán chu trình nhiệt động với sai số kiểm nghiệm so với đồ thị i-c dưới mức 1,5%.
  • Xác lập định lượng quy luật biến thiên của hệ số sử dụng nhiệt COP (từ 0,06 đến 0,36) và bội số tuần hoàn f (từ 2,21 đến 59,83) dưới tác động của nhiệt độ nguồn nhiệt cấp, nhiệt độ bay hơi và nhiệt độ nước giải nhiệt.
  • Chứng minh tính khả thi vượt trội của giải pháp tận dụng nguồn nhiệt thải có nhiệt thế thấp từ 80°C đến 150°C để cấp đông và làm lạnh công nghiệp, bảo vệ tầng ozone và giảm phát thải khí nhà kính.
  • Đặt nền móng kỹ thuật vững chắc cho việc thiết kế, chế tạo và nội địa hóa thiết bị máy lạnh hấp thụ công nghiệp tại Việt Nam trong giai đoạn 12 đến 24 tháng tới.

Đóng góp chính của luận văn là đã số hóa thành công quy trình tính toán nhiệt động phức tạp thành công cụ phần mềm tiện ích, tháo gỡ điểm nghẽn về mặt công nghệ thiết kế máy lạnh hấp thụ tại Việt Nam.

Các đơn vị sản xuất công nghiệp và doanh nghiệp cơ điện lạnh nên chủ động ứng dụng ngay mô hình tính toán này trong vòng 6 đến 12 tháng tới để tối ưu hóa việc thu hồi nhiệt thải, cắt giảm chi phí vận hành và đón đầu xu thế công nghệ xanh bền vững.