MÔ PHỎNG SỐ QUÁ TRÌNH LÀM MÁT MÔI CHẤT CO2 BẰNG NƯỚC: GIẢI PHÁP TIÊN PHONG CHO CÔNG NGHỆ NHIỆT LẠNH XANH TẠI VIỆT NAM

Tóm tắt nghiên cứu (Research Summary)

Nghiên cứu "Mô phỏng số quá trình làm mát $CO_2$ bằng nước" do nhóm nghiên cứu sinh viên Khoa Cơ khí Động lực, Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh thực hiện dưới sự hướng dẫn của PGS.TS Đặng Thành Trung, tập trung giải quyết bài toán cốt lõi: Làm thế nào để nâng cao hiệu quả giải nhiệt của môi chất lạnh tự nhiên $CO_2$ (R744) bằng nguồn nước tuần hoàn trong điều kiện khí hậu nhiệt đới, và mức độ tin cậy của phương pháp mô phỏng số đa vật lý trong việc dự đoán đặc tính nhiệt động lực học của thiết bị trao đổi nhiệt?

Phương pháp luận của đề tài kết hợp chặt chẽ giữa tính toán nhiệt động lý thuyết chu trình nén hơi một cấp và kỹ thuật mô phỏng phần tử hữu hạn (FEM) trên nền tảng COMSOL Multiphysics 5.6. Mô hình dàn trao đổi nhiệt ống đồng làm mát bay hơi được thiết kế chi tiết với công suất giải nhiệt danh định $1700\text{ W}$, bao gồm 5 tép dàn với tổng cộng 125 ống trao đổi nhiệt vi mô ($\varnothing 4\text{ mm}$, dài $300\text{ mm}$). Quá trình trao đổi nhiệt liên hợp và trường thủy động lực học được mô phỏng chính xác bằng mô hình dòng chảy rối $k-\varepsilon$.

Kết quả then chốt chỉ ra rằng dòng $CO_2$ ở áp suất ngưng tụ $74\text{ bar}$ với nhiệt độ đầu vào $65^\circ\text{C}$ sau khi qua dàn giải nhiệt bằng nước ở $28^\circ\text{C}$ đã hạ xuống mức $31,3^\circ\text{C}$ tại ngõ ra, với độ tổn thất áp suất chỉ $0,228\text{ bar}$. Độ sai lệch nhiệt độ giữa tính toán lý thuyết ($30^\circ\text{C}$) và mô phỏng số ($31,3^\circ\text{C}$) chỉ là $1,3^\circ\text{C}$ (sai số tương đối dưới 4,3%), chứng minh độ tin cậy vượt trội của mô hình. Nghiên cứu mang lại ý nghĩa học thuật và ứng dụng to lớn, khẳng định tính khả thi của việc ứng dụng $CO_2$ làm môi chất lạnh xanh hiệu suất cao tại Việt Nam, mở ra hướng đi đột phá thay thế hoàn toàn các chất làm lạnh truyền thống gây suy giảm tầng ôzôn và biến đổi khí hậu.

Bối cảnh và tầm quan trọng (Background & Significance)

Ngành công nghiệp nhiệt lạnh và điều hòa không khí (HVAC) toàn cầu đang đứng trước một bước ngoặt lịch sử khi các hiệp định quốc tế như Nghị định thư Montreal và Bản sửa đổi, bổ sung Kigali bắt buộc lộ trình loại bỏ dần các chất làm lạnh nhân tạo nhóm CFC, HCFC và HFC. Dù sở hữu đặc tính nhiệt động thuận lợi, các dòng môi chất này lại mang tiềm năng làm nóng lên toàn cầu (GWP) gấp hàng nghìn lần khí $CO_2$ và đe dọa trực tiếp đến tầng ôzôn. Thống kê tại các quốc gia phát triển cho thấy các hệ thống HVAC trong các tòa nhà chiếm tới gần 39% tổng năng lượng tiêu thụ cuối cùng và chịu trách nhiệm cho phần lớn lượng khí thải nhà kính trong lĩnh vực dân dụng.

Trước áp lực môi sinh khẩn cấp, $CO_2$ (R744) nổi lên như một ứng cử viên thay thế hoàn hảo:

  • Chỉ số phá hủy tầng ôzôn bằng 0 ($\text{ODP} = 0$).
  • Tiềm năng làm nóng toàn cầu tối thiểu ($\text{GWP} = 1$).
  • Không cháy nổ, không độc hại, thân thiện tuyệt đối với môi trường.
  • Trữ lượng tự nhiên vô tận với chi phí sản xuất thấp.

Tuy nhiên, rào cản kỹ thuật lớn nhất của môi chất $CO_2$ là nhiệt độ tới hạn rất thấp ($31,1^\circ\text{C}$) ở áp suất tới hạn $73,8\text{ bar}$. Khi nhiệt độ môi trường xung quanh vượt quá ngưỡng này—đặc biệt tại các khu vực nhiệt đới gió mùa như TP. Hồ Chí Minh—hệ thống bắt buộc phải chuyển sang chu trình siêu tới hạn (transcritical cycle) với áp suất vận hành rất cao, dẫn đến suy giảm nghiêm trọng hệ số hiệu quả năng lượng (COP) nếu chỉ giải nhiệt bằng không khí thông thường.

Tại Việt Nam, các công bố khoa học chuyên sâu về hệ thống lạnh $CO_2$ còn rất hạn chế, phần lớn công nghệ cốt lõi vẫn phụ thuộc vào chuyển giao từ nước ngoài. Khoảng trống tri thức này đòi hỏi những công trình nghiên cứu bản địa hóa, đánh giá chính xác khả năng giải nhiệt cho $CO_2$ bằng môi trường nước (với nhiệt dung riêng và hệ số tỏa nhiệt cao hơn không khí nhiều lần). Nghiên cứu của nhóm tác giả ĐH Sư phạm Kỹ thuật TP.HCM xuất hiện đúng thời điểm, giải quyết trực tiếp thách thức kiểm soát nhiệt độ đầu ra của $CO_2$ tiệm cận điểm tới hạn bằng phương pháp giải nhiệt nước, tạo nền tảng vững chắc cho quá trình chuyển dịch công nghệ lạnh bền vững tại nước ta.

Methodology và Phương pháp tiếp cận (Methodology & Approach)

Đề tài áp dụng quy trình nghiên cứu tích hợp đa giai đoạn, kết hợp giải tích lý thuyết nhiệt động lực học, dựng hình tham số không gian ba chiều và mô phỏng số động lực học lưu chất (CFD):

1. Thiết kế và tính toán nhiệt động lý thuyết

  • Xác lập chu trình nén hơi một cấp: Dựa trên các thông số làm việc tiêu chuẩn: công suất giải nhiệt $Q_k = 1700\text{ W}$, nhiệt độ ngưng tụ $t_k = 31^\circ\text{C}$ ($p_k = 74\text{ bar}$), nhiệt độ bay hơi $t_0 = 10^\circ\text{C}$ ($p_0 = 45\text{ bar}$), nhiệt độ quá nhiệt $t_{qn} = 22^\circ\text{C}$, nhiệt độ quá lạnh $t_{ql} = 1^\circ\text{C}$. Tỉ số nén tính toán đạt $\pi = 1,64 < 9$, chứng minh tính tối ưu tuyệt đối của chu trình 1 cấp.
  • Xác định diện tích trao đổi nhiệt: Sử dụng phương pháp hiệu số nhiệt độ trung bình logarit ($\Delta t_{tb} = 12^\circ\text{C}$) và hệ số truyền nhiệt tổng $k = 300\text{ W/m}^2\text{K}$, tính toán được diện tích bề mặt cần thiết $F = 0,47\text{ m}^2$.
  • Cấu hình hình học: Lựa chọn ống đồng chuyên dụng chịu áp suất cao có đường kính ngoài $\varnothing 4\text{ mm}$, chiều dài mỗi ống $L = 30\text{ cm}$. Số lượng ống tính toán lý thuyết là 124,6 ống, được quy chuẩn thành 125 ống phân bố đều trên 5 tép dàn (mỗi tép gồm 25 ống song song) đấu nối qua ống góp $\varnothing 10\text{ mm}$ nhằm tối ưu hóa phân phối lưu lượng và giảm thiểu trở lực ma sát.

2. Mô hình hóa số học với COMSOL Multiphysics 5.6

Nghiên cứu sử dụng phương pháp phần tử hữu hạn để giải quyết bài toán ghép đôi đa vật lý phức tạp (Multiphysics Coupling):

  • Truyền nhiệt liên hợp (Conjugate Heat Transfer): Giải đồng thời phương trình truyền nhiệt trong chất rắn (vách ống đồng) và truyền nhiệt đối lưu trong dòng lưu chất (khối nước giải nhiệt và dòng gas $CO_2$).
  • Mô hình thủy động lực học dòng chảy rối: Ứng dụng mô hình $k-\varepsilon$ chuẩn với xử lý biên thành (wall functions), giải quyết hệ phương trình bảo toàn khối lượng (phương trình liên tục), bảo toàn động lượng (Navier-Stokes) và bảo toàn năng lượng: $$\rho (\mathbf{u} \cdot \nabla)\mathbf{u} = \nabla \cdot [-p\mathbf{I} + (\mu + \mu_T)(\nabla\mathbf{u} + (\nabla\mathbf{u})^T)] + \mathbf{F}$$ $$\rho C_p \mathbf{u} \cdot \nabla T = \nabla \cdot (k \nabla T) + Q$$
  • Thiết lập điều kiện biên chính xác:
    • Đầu vào dòng $CO_2$ (Inflow): Lưu lượng khối lượng $G_0 = 0,009\text{ kg/s}$, nhiệt độ vào $T_{in} = 65^\circ\text{C}$ (khảo sát mở rộng tại $68^\circ\text{C}$).
    • Khối chất lỏng giải nhiệt: Nước duy trì ổn định ở $T_{water} = 28^\circ\text{C}$ (và $26^\circ\text{C}$).
    • Áp suất đầu ra (Pressure Outlet): $P_k = 74\text{ bar}$.
    • Điều kiện biên vách ống: Điều kiện không trượt (no-slip condition) $u = 0$.

Quy trình giải bài toán dừng (Stationary Study) kết hợp chia lưới tam giác tự do (Free Tetrahedral Mesh) đa tỷ lệ, tinh chỉnh mật độ lưới tại lớp biên sát vách ống đồng nhằm đảm bảo độ hội tụ và tính chính xác cao nhất cho gradient vận tốc và nhiệt độ.

Các phát hiện chính (Key Findings)

Từ các kết quả mô phỏng trường nhiệt độ, phân bố áp suất và đối chiếu lý thuyết, nghiên cứu ghi nhận 4 phát hiện khoa học mang tính bước ngoặt:

Thông số đánh giá Giá trị tính toán lý thuyết Kết quả mô phỏng số Sai lệch tuyệt đối Đánh giá kỹ thuật
Nhiệt độ $CO_2$ đầu ra ($T_{out}$) $30,0^\circ\text{C}$ $31,3^\circ\text{C}$ $+1,3^\circ\text{C}$ Độ chính xác đạt $95,7%$, nằm trong dung sai kỹ thuật
Nhiệt độ nước giải nhiệt $28,0^\circ\text{C}$ $28,0^\circ\text{C}$ $0,0^\circ\text{C}$ Duy trì ổn định theo điều kiện biên
Áp suất làm việc ngõ vào $74,000\text{ bar}$ $74,228\text{ bar}$ $+0,228\text{ bar}$ Xuất hiện tổn thất áp suất động lực học
Áp suất làm việc ngõ ra $74,000\text{ bar}$ $74,000\text{ bar}$ $0,000\text{ bar}$ Khớp hoàn toàn thông số chu trình
Tổn thất áp suất ($\Delta P$) $0\text{ bar}$ (lý tưởng) $0,228\text{ bar}$ $0,228\text{ bar}$ Trở lực thủy lực cực thấp, bảo toàn năng lượng bơm

Chi tiết các phát hiện:

  1. Hiệu quả hạ nhiệt ấn tượng qua cấu trúc 5 tép dàn: Dòng môi chất $CO_2$ với nhiệt độ đầu vào lên tới $65^\circ\text{C}$ được hạ nhiệt nhanh chóng và đồng đều qua từng tầng dàn ống. Mặt cắt nhiệt độ qua các tép dàn 1 đến 5 minh chứng sự truyền nhiệt cưỡng bức diễn ra vô cùng mãnh liệt ở 2 tép đầu tiên nhờ chênh lệch nhiệt độ lớn ($\Delta T_{max} = 37^\circ\text{C}$), sau đó ổn định dần và đạt $31,3^\circ\text{C}$ ở đầu ra tép thứ 5. Mức nhiệt này xấp xỉ nhiệt độ ngưng tụ thiết kế, bảo đảm môi chất hóa lỏng hoàn toàn trước khi đi vào van tiết lưu.
  2. Kiểm soát tổn thất áp suất tối ưu: Một trong những lo ngại lớn nhất khi sử dụng dàn ống đường kính nhỏ ($\varnothing 4\text{ mm}$) là hiện tượng sụt áp ma sát cao dẫn đến lãng phí công nén. Kết quả mô phỏng cho thấy tổn thất áp suất tổng dọc chiều dài toàn bộ hệ thống chỉ là $0,228\text{ bar}$ (tương đương mức sụt giảm chỉ $0,3%$ áp suất vận hành). Điều này xác nhận rằng việc chia dòng $CO_2$ thành 25 nhánh song song trên mỗi tép dàn đã triệt tiêu hiệu quả ma sát dòng chảy rối mà vẫn giữ được diện tích tiếp xúc nhiệt tối đa.
  3. Sự tương quan xuất sắc giữa mô hình số và giải tích: Độ chênh lệch nhiệt độ đầu ra giữa mô phỏng ($31,3^\circ\text{C}$) và tính toán lý thuyết ($30,0^\circ\text{C}$) chỉ là $1,3^\circ\text{C}$. Độ sai lệch nhỏ này phản ánh trung thực thực tế cơ học chất lỏng: tính toán lý thuyết thường giả định hệ số truyền nhiệt không đổi và bỏ qua điện trở nhiệt tiếp xúc cục bộ, trong khi mô phỏng số xét đến sự suy giảm lớp biên nhiệt và biến thiên độ nhớt động học dọc theo đường ống.
  4. Tính khả quan của giải pháp giải nhiệt bay hơi kết hợp nước: Nghiên cứu chứng minh rằng chỉ cần nguồn nước làm mát thông thường ở mức $28^\circ\text{C}$ (hoàn toàn khả thi với tháp giải nhiệt cooling tower tại miền Nam Việt Nam), hệ thống đã có thể đưa môi chất $CO_2$ về trạng thái ngưng tụ ổn định mà không cần vận hành chu trình siêu tới hạn phức tạp và tốn kém điện năng.

Đóng góp khoa học và giá trị ứng dụng (Scientific Contributions)

1. Đóng góp về mặt lý thuyết và học thuật

  • Xây dựng bộ thông số nhiệt động hoàn chỉnh cho chu trình lạnh $CO_2$ một cấp ứng dụng phương pháp giải nhiệt nước trong điều kiện vận hành áp suất cao ($74\text{ bar}$).
  • Xác thực tính chuẩn xác của mô hình chảy rối $k-\varepsilon$ khi giải bài toán truyền nhiệt liên hợp chất lỏng - chất rắn cho môi chất có mật độ cao và độ nhớt thấp như $CO_2$ ở trạng thái bão hòa tiệm cận điểm tới hạn.

2. Đổi mới về phương pháp nghiên cứu (Methodological Innovations)

  • Chuẩn hóa quy trình thiết kế khép kín cho thiết bị trao đổi nhiệt áp lực cao: Từ giải tích toán học $\rightarrow$ Thiết kế 3D CAD tham số hóa $\rightarrow$ Phân tích số CFD đa vật lý $\rightarrow$ Đánh giá đối chuẩn.
  • Phương pháp này cho phép loại bỏ các thử nghiệm thực tế tốn kém, rút ngắn thời gian nghiên cứu và phát triển sản phẩm (R&D), giảm thiểu rủi ro nổ vỡ khi làm việc với môi chất áp suất cao.

3. Khả năng ứng dụng thực tiễn và chính sách

  • Công nghiệp chế tạo thiết bị lạnh nội địa: Cung cấp bản vẽ kỹ thuật và thông số chuẩn cho các doanh nghiệp sản xuất dàn ngưng tụ, tháp làm mát bằng nước cho hệ thống máy lạnh $CO_2$ thương mại và công nghiệp.
  • Đóng góp vào mục tiêu Net-Zero: Khẳng định tiềm năng hiện thực hóa các dòng máy lạnh "xanh", hỗ trợ đắc lực cho chiến lược quốc gia về biến đổi khí hậu, hướng tới mục tiêu phát thải ròng bằng 0 vào năm 2050.
  • Công bố khoa học uy tín: Kết quả nghiên cứu đã được phản biện độc lập và chấp nhận công bố chính thức tại Kỷ yếu Hội nghị Khoa học Quốc tế về Công nghệ Xanh và Phát triển Bền vững (GTSD 2022).

Đối tượng đón nhận nghiên cứu (Target Audience & Benefits)

  • Giới nghiên cứu học thuật và giảng viên: Nhận được nguồn tài liệu tham khảo thực nghiệm - mô phỏng có giá trị cao cho các bộ môn Kỹ thuật Nhiệt, Điện lạnh, Động lực học chất lưu tính toán (CFD). Cung cấp case study mẫu về khai thác phần mềm COMSOL Multiphysics trong phân tích thiết bị trao đổi nhiệt.
  • Kỹ sư thiết kế và R&D ngành HVAC-R: Nắm bắt được phương pháp tính toán kết cấu ống góp, số lượng ống chùm và dự báo chính xác trường nhiệt - áp của môi chất $CO_2$, từ đó tự tin chuyển đổi các hệ thống sử dụng gas freon (R22, R134a, R410A) sang gas $CO_2$.
  • Doanh nghiệp vận hành kho lạnh và siêu thị: Tiếp cận giải pháp công nghệ làm lạnh tự nhiên an toàn, tiết kiệm chi phí năng lượng nhờ tận dụng nước tuần hoàn, đồng thời đáp ứng các tiêu chuẩn chứng chỉ công trình xanh (LEED, LOTUS).
  • Cơ quan quản lý và hoạch định chính sách: Có thêm căn cứ khoa học vững chắc để ban hành các quy chuẩn, tiêu chuẩn kỹ thuật thúc đẩy việc ứng dụng chất làm lạnh tự nhiên tại Việt Nam.

Câu hỏi thường gặp (Frequently Asked Questions - FAQ)

1. Phát hiện quan trọng nhất của nghiên cứu này là gì?

Phát hiện cốt lõi là việc khẳng định tính khả thi vượt trội của phương pháp giải nhiệt $CO_2$ bằng nguồn nước $28^\circ\text{C}$. Hệ thống đã hạ nhiệt độ $CO_2$ từ $65^\circ\text{C}$ xuống $31,3^\circ\text{C}$ với mức sụt áp cực thấp ($0,228\text{ bar}$), chứng minh môi chất có thể hóa lỏng an toàn trong chu trình dưới tới hạn mà không cần cấu hình siêu tới hạn phức tạp.

2. Phương pháp mô phỏng số trong đề tài có điểm gì đặc biệt?

Nghiên cứu sử dụng COMSOL Multiphysics 5.6 với module Conjugate Heat Transfer kết hợp mô hình dòng chảy rối $k-\varepsilon$. Điểm nổi bật là mô phỏng đồng thời cả 3 hiện tượng: dẫn nhiệt qua thành ống đồng dày chịu áp, dòng đối lưu cưỡng bức của $CO_2$ áp suất cao bên trong ống và truyền nhiệt bề mặt với khối nước bên ngoài.

3. Kết quả mô phỏng này có thể áp dụng cho các quy mô công suất khác không?

Có. Phương pháp tính toán số ống ($n = 125$), hệ số truyền nhiệt ($k = 300\text{ W/m}^2\text{K}$) và phương pháp NTU/Log-mean được thiết lập dưới dạng tham số hóa. Các kỹ sư có thể dễ dàng mở rộng quy mô (scale-up) cho các hệ thống công suất lớn hơn (hàng chục đến hàng trăm kW) bằng cách tăng số tép dàn hoặc chiều dài ống theo công thức đã đề xuất.

4. Hướng phát triển tiếp theo của nghiên cứu là gì?

Nhóm nghiên cứu đề xuất tiến hành chế tạo mô hình thực nghiệm vật lý tại phòng thí nghiệm để kiểm chứng chéo với dữ liệu mô phỏng, đồng thời nghiên cứu chuyên sâu về quá trình sôi dòng hai pha (two-phase flow boiling) và tối ưu hóa thêm hình dạng cánh tản nhiệt nhằm nâng cao hơn nữa hệ số tỏa nhiệt của dàn ống.

5. Khí $CO_2$ có áp suất làm việc lên tới 74 bar, liệu có an toàn trong thực tế?

Áp suất $74\text{ bar}$ là rất cao so với các hệ thống lạnh truyền thống ($15 - 25\text{ bar}$). Tuy nhiên, thiết kế dàn ống đồng đường kính nhỏ ($\varnothing 4\text{ mm}$) có độ bền chịu áp nổ rất lớn (theo công thức vỏ mỏng chịu áp lực), đảm bảo hệ số an toàn cao trong vận hành công nghiệp.

Kết luận (Conclusion)

Báo cáo nghiên cứu khoa học của nhóm tác giả Trường ĐH Sư phạm Kỹ thuật TP.HCM là một công trình tiêu biểu, kết hợp nhuần nhuyễn giữa nền tảng lý thuyết nhiệt - lạnh vững chắc và công nghệ mô phỏng số tiên tiến COMSOL Multiphysics. Bằng việc chứng minh khả năng làm mát dòng $CO_2$ đạt $31,3^\circ\text{C}$ với sai số dưới 4,3% so với lý thuyết, đề tài đã phá vỡ rào cản ứng dụng môi chất tự nhiên trong điều kiện khí hậu nóng ẩm tại Việt Nam.

Trong bối cảnh chuyển dịch năng lượng và bảo vệ môi trường đang trở thành mệnh lệnh sống còn, các kết quả từ công trình này chính là bước đệm then chốt. Nhóm nghiên cứu kêu gọi các viện nghiên cứu, trường đại học và cộng đồng doanh nghiệp tiếp tục đầu tư nguồn lực, đẩy mạnh thử nghiệm thực tế và sớm thương mại hóa các hệ thống điều hòa không khí và bơm nhiệt sử dụng gas sạch $CO_2$, đưa nền kỹ thuật nhiệt lạnh nước nhà bắt nhịp cùng xu thế phát triển xanh của thế giới.