Tổng quan về luận án

Ngành kỹ thuật nhiệt và điều hòa không khí (ĐHKK) toàn cầu đang đối mặt với thách thức kép: vừa phải đáp ứng nhu cầu tiện nghi vi khí hậu và làm mát công nghệ ngày càng gia tăng, vừa phải triệt tiêu các tác nhân gây suy giảm tầng ô-zôn cùng phát thải khí nhà kính. Theo báo cáo đánh giá khoa học của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC, 2014), các khí nhà kính Flo hóa (khí F) chiếm khoảng 2% nguyên nhân gây nóng lên toàn cầu; trong đó Hydrofluorocarbon (HFCs) chiếm tới 91%, Perfluorocarbon (PFCs) chiếm 6% và SF6 chiếm 3%. Đáng chú ý, có tới "86% lượng HFC được sử dụng làm môi chất lạnh trong ngành lạnh và ĐHKK (với tỷ lệ phân bổ 35% cho máy lạnh và 65% cho ĐHKK tính theo GWP tương đương)". Trước các chế tài pháp lý quốc tế nghiêm ngặt như Nghị định thư Montreal, Bản sửa đổi bổ sung Kigali (2016) và Quy định F-gas 517/2014 của Liên minh Châu Âu (cấm các môi chất có GWP ≥ 150 trên ô tô từ năm 2017 và ĐHKK thương mại từ 2020–2025), cùng Quyết định số 1888/QĐ-BTNMT của Việt Nam (loại bỏ 1.000 tấn HCFC-22 trong giai đoạn 2020–2024), việc chuyển dịch sang môi chất lạnh (MCL) thế hệ thứ tư dòng Hydrofluoro-olefin (HFO) là đòi hỏi cấp thiết.

Khoảng trống nghiên cứu then chốt (Research Gap) nằm ở sự thiếu hụt nghiêm trọng bộ cơ sở dữ liệu nhiệt động liên tục, toàn diện và có độ chính xác cao cho các đồng phân HFO tiềm năng, đặc biệt là R1234ze(Z) và R1243zf. Dữ liệu thực nghiệm công bố trên thế giới còn rời rạc, chủ yếu chỉ giới hạn ở một số điểm áp suất bão hòa hoặc khối lượng riêng cục bộ, không thể đo trực tiếp entropy và entanpy trong toàn bộ dải thông số nhiệt động. Hơn nữa, các phương trình trạng thái dạng khối bậc ba truyền thống (Cubic EOS) như Peng-Robinson hay Patel-Teja bộc lộ sai số rất lớn khi áp dụng cho phân tử phân cực HFO (sai lệch nhiệt dung riêng đẳng áp $C_p$ lên tới 20%). Ngoài ra, chưa có một công trình nào tích hợp đồng bộ từ việc mô hình hóa năng lượng tương tác phân tử, thiết lập bộ dữ liệu nhiệt động hoàn chỉnh, mô phỏng chu trình lý thuyết đa chế độ đến đánh giá thực nghiệm trên hệ thống ĐHKK thực tế trong điều kiện khí hậu nhiệt đới tại Việt Nam.

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập chính xác các thông số nhiệt động của môi chất HFO cấu trúc phức tạp khi số liệu thực nghiệm bị hạn chế? Giả thuyết 1 (H1): Mô hình phương trình trạng thái dựa trên năng lượng tự do Helmholtz kết hợp lực tương tác phân tử (MBEOS) và thuật toán tối ưu hóa đơn hình phi tuyến Nelder-Mead có thể nội suy và ngoại suy chính xác toàn bộ trường dữ liệu nhiệt động của R1234ze(Z) và R1243zf với sai số tuyệt đối trung bình (AAD) < 0,1% so với thực nghiệm.
  2. Câu hỏi nghiên cứu 2 (RQ2): Mức độ tương thích và hiệu quả nhiệt động của các đồng phân HFO (R1234yf, R1234ze(E), R1243zf, R1234ze(Z)) khi thay thế trực tiếp (drop-in replacement) cho R22 và R134a ra sao? Giả thuyết 2 (H2): Môi chất R1243zf và R1234ze(E) có đường bão hòa, áp suất làm việc và tỷ số nén tương đương R134a/R22, cho phép duy trì hệ số hiệu quả năng lượng (COP) suy giảm không quá 2,5% trên chu trình lý thuyết.
  3. Câu hỏi nghiên cứu 3 (RQ3): Hiệu suất nén không thuận nghịch và tổn thất áp suất thủy lực trên đường ống ảnh hưởng như thế nào đến COP thực tế của hệ thống chiller giải nhiệt gió sử dụng HFO? Giả thuyết 3 (H3): Tổn thất áp suất và sự suy giảm hiệu suất nén đẳng entropy thực nghiệm từ 100% xuống 53%–63% sẽ làm giảm COP của chu trình thực tế từ 30% đến 45% so với chu trình lý tưởng.
  4. Câu hỏi nghiên cứu 4 (RQ4): Nhiệt độ nước lạnh làm mát ảnh hưởng thế nào đến sự đáp ứng tải nhiệt và hiệu quả biến đổi năng lượng của R1234ze(E) so với R134a? Giả thuyết 4 (H4): Tăng nhiệt độ nước lạnh cung cấp từ 8°C lên 12°C sẽ cải thiện COP của hệ thống sử dụng R1234ze(E) vượt trội hơn (19%) so với môi chất nền R134a (15%).

Khung lý thuyết của nghiên cứu được neo trên Lý thuyết Thống kê Nhiệt động học Cổ điển và Hiện đại (Statistical Associating Fluid Theory - SAFT), Phương trình Trạng thái cơ bản Helmholtz ($F = F^0 + F^r$) và Nhiệt động học Chu trình Lạnh nén hơi ngược chiều. Quy mô nghiên cứu bao quát khảo sát 19 đồng phân propen flo hóa, mô hình hóa chuyên sâu 2 môi chất thiếu hụt dữ liệu (R1234ze(Z), R1243zf), mô phỏng nhiệt động 9 môi chất (4 HFO và 5 môi chất tham chiếu: R134a, R22, R32, R152a, R245fa) trên dải nhiệt độ bay hơi từ -5°C đến 15°C, cùng hệ thống thực nghiệm chiller công suất chuẩn thử nghiệm trên 140 chế độ vận hành ổn định.

Literature Review và Positioning

Lịch sử phát triển môi chất lạnh trải qua bốn thế hệ mang tính quy luật:

  • Thế hệ thứ nhất (1830–1930): Sử dụng các hợp chất tự nhiên sẵn có ($SO_2, CO_2, NH_3$, Ether) với nhược điểm độc hại, dễ cháy nổ và áp suất làm việc khắc nghiệt.
  • Thế hệ thứ hai (1931–1990): Sự thống trị của các hợp chất chlorofluorocarbon (CFC) và hydrochlorofluorocarbon (HCFC) như R11, R12, R22, R114. Tuy nhiên, phát hiện chấn động năm 1974 của hai nhà khoa học đoạt giải Nobel là Rowland và Molina đã chỉ rõ nguyên tử Clo trong CFC/HCFC là tác nhân xúc tác phá hủy tầng ô-zôn bình lưu.
  • Thế hệ thứ ba (1990–2010): Sự chuyển dịch sang hydrofluorocarbon (HFC) với $ODP = 0$ điển hình là R134a, R404A, R410A, R407C, R507A. Mặc dù bảo vệ tầng ô-zôn, nhóm HFC lại có chỉ số tiềm năng làm nóng toàn cầu rất cao (GWP của R134a là 1.430; R410A là 2.088; R125 là 3.500), thúc đẩy phát thải khí nhà kính gián tiếp và trực tiếp.
  • Thế hệ thứ tư (2010–nay): Kỷ nguyên của môi chất lạnh thân thiện môi trường ($ODP = 0$, $GWP < 10$), trọng tâm là dòng Hydrofluoro-olefin (HFO) và các chất tự nhiên cải tiến.
TIẾN TRÌNH PHÁT TRIỂN CÁC THẾ HỆ MÔI CHẤT LẠNH

Trong dòng văn hiến quốc tế, tồn tại hai luồng quan điểm tranh luận khoa học sâu sắc về phương pháp xác định tính chất nhiệt động của môi chất mới:

  1. Luồng quan điểm thực nghiệm thuần túy và phương trình bán thực nghiệm bậc ba: Tiêu biểu là các nghiên cứu của Katsuyuki (2010) sử dụng phương pháp chiết tách (extraction method) để đo áp suất bão hòa của R1234ze(E) (526–3.462 kPa) và R1234ze(Z) (185–2.309 kPa); Brown (2010) công bố 83 điểm dữ liệu pvT của R1243zf (233,80 K đến 372,86 K). Tuy nhiên, nhóm nghiên cứu của Jacobsen (1997) và Akasaka (2010) đã chỉ ra rằng khi áp dụng phương trình Peng-Robinson hay Patel-Teja cho R1234yf, sai lệch khối lượng riêng pha lỏng lên tới 2,5% và sai số nhiệt dung riêng đẳng áp $C_p$ vượt quá 20%, không thể phản ánh chính xác tương tác dị hướng và mômen đa cực của phân tử phân cực.
  2. Luồng quan điểm mô hình hóa cơ sở phân tử (Molecular-based EOS): Tiêu biểu là nghiên cứu của Lemmon (2009, NIST) phát triển phương trình đa tham số cho Propane và Dimethyl Carbonate, và Akasaka (2011) cho R1234ze(Z). Dù đạt độ chính xác cực cao ($AAD < 0,02%$), các phương trình này đòi hỏi hàng trăm điểm đo thực nghiệm phân bố đều ở mọi pha – một điều bất khả thi đối với các môi chất mới tổng hợp.

Công trình này định vị chính xác tại giao điểm chiến lược: Ứng dụng lý thuyết tương tác phân tử Helmholtz (Müller, 1995; Gross & Sadowski, 2001 - PC-SAFT) cho chất lỏng phân tử nhỏ gọn có tính dị hướng và phân cực, giải quyết triệt để bài toán thiếu dữ liệu mà vẫn đảm bảo độ chính xác vượt trội so với các PTTT bậc ba kinh điển.

Khi đối chiếu với các công trình quốc tế cùng thời kỳ:

  • So sánh với Sad Jarall (2014): Tác giả Jarall khảo sát chu trình nén hơi có hồi nhiệt sử dụng R1234yf thay thế R134a trên máy nén rô-to lăn 550 W, kết luận R1234yf có COP thấp hơn R134a từ 4% đến 9%. Luận án mở rộng phạm vi vượt bậc khi không chỉ khảo sát R1234yf mà còn chứng minh R1243zf và R1234ze(E) có tiềm năng thay thế vượt trội hơn hẳn R1234yf (COP của R1243zf chỉ thấp hơn R134a 0,52%).
  • So sánh với Claudio Zilio (2011) và AB Pearson (2014): Zilio thử nghiệm R1234yf trên hệ thống ĐHKK ô tô với >140 chế độ, trong khi Pearson nghiên cứu R1234ze trên chiller ly tâm 250 kW ($T_0 = 5^\circ\text{C}, T_k = 50^\circ\text{C}$). Luận án đã tiến xa hơn khi phân tích định lượng đồng thời cả hai yếu tố suy giảm phi lý tưởng: hiệu suất nén cơ học thực tế ($\eta = 53%–63%$) và trở lực thủy lực đường ống ($\Delta P = 25\text{ kPa}$ và $63\text{ kPa}$) trên mô hình thực nghiệm hoàn chỉnh, mang lại giá trị tham chiếu thực tiễn cao cho điều kiện Đông Nam Á.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và làm sâu sắc Lý thuyết Tương tác Phân tử Helmholtz trong việc xác định các đặc tính nhiệt động liên tục của môi chất lạnh HFO phân cực. Cấu trúc năng lượng tự do Helmholtz ($F$) được phân rã tường minh thành hai thành phần cơ bản: phần khí lý tưởng ($F^0$) và phần thực dư thừa ($F^r$):

$$F(T, \rho) = F^0(T, \rho) + F^r(T, \rho)$$

Trong đó, phần thực dư thừa $F^r$ được tích hợp từ ba dạng năng lượng tương tác phân tử vi mô độc lập:

  1. Năng lượng tương tác lõi rắn đẩy (Hard-body interaction, $F^H$): Biểu diễn lực đẩy tĩnh điện cản trở giữa các đám mây electron của các phần tử phân tử khi tiếp cận khoảng cách nhỏ hơn đường kính tương đương $\sigma$.
  2. Năng lượng do lực hấp dẫn phân tán (Attractive dispersion force, $F^A$): Biểu diễn lực hút Van der Waals giữa các phân tử theo thế năng Lennard-Jones cải tiến, phụ thuộc vào chiều sâu giếng thế năng $\epsilon/k_B$ và hệ số dị hướng phân tử $\alpha$.
  3. Năng lượng do tương tác đa cực điện (Quadrupolar interaction, $F^Q$): Đóng góp đột phá của mô hình Müller & Gubbins nhằm mô tả chính xác sự phân bố điện tích bất đối xứng của liên kết đôi $C=C$ và các nguyên tử Flo có độ âm điện cực cao trong cấu trúc HFO thông qua mômen tứ cực không thứ nguyên $Q^{*2}$.

$$F^r = F^H + F^A + F^Q$$

KHUNG CẤU TRÚC PHÂN TÍCH NĂNG LƯỢNG TỰ DO HELMHOLTZ

Mô hình lý thuyết thiết lập mối quan hệ giải tích chặt chẽ giữa hàm thế vi mô và các thông số trạng thái vĩ mô: áp suất $p = \rho^2 (\partial F/\partial \rho)T$, nội năng $u = F - T(\partial F/\partial T)\rho$, entanpy $h = u + p/\rho$, entropy $s = -(\partial F/\partial T)\rho$, và nhiệt dung riêng đẳng tích $C_v = -T(\partial^2 F/\partial T^2)\rho$.

Bằng cách chuyển dịch từ mô hình hóa hồi quy thực nghiệm sang mô hình hóa cơ sở tương tác hạt vi mô, luận án đã tạo ra bước tiến quan trọng: giải thích được bản chất vật lý của sự biến thiên tính chất nhiệt động vùng bão hòa và quá nhiệt thông qua 4 đại lượng đặc trưng phân tử $(T_0, \rho_0, \alpha, Q^{*2})$.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên hoàn 3 lý thuyết chuyên sâu:

  1. Lý thuyết cấu trúc phân tử và đóng góp nhóm (Joback-Reid cải tiến): Dự đoán nhiệt dung riêng khí lý tưởng $C_p^0(T) = a_i + b_i T + c_i T^2 + d_i T^3$ và ước lượng sơ bộ các thông số tới hạn $(T_c, P_c, \rho_c)$.
  2. Lý thuyết tương tác đa cực Helmholtz (MBEOS): Giải bài toán tối ưu hóa phi tuyến đa biến bằng thuật toán đơn hình Nelder-Mead để trích xuất bộ thông số phân tử tối ưu từ tập dữ liệu thực nghiệm giới hạn.
  3. Lý thuyết chu trình nhiệt động học thực tế (Exergy-Enthalpy Analysis): Đánh giá hiệu quả biến đổi năng lượng COP, năng suất lạnh riêng thể tích $q_v$, tỷ số nén $\Pi$, nhiệt độ đầu đẩy máy nén $T_d$ có xét đến suy hao áp suất đường ống và hiệu suất nén đẳng entropy thực tế.

Điều kiện biên xác lập cho khung phân tích: Áp dụng cho các phân tử đơn chất HFO mạch hở có liên kết đôi, khối lượng mol từ 90 đến 140 g/mol, dải nhiệt độ làm việc từ $-30^\circ\text{C}$ đến $150^\circ\text{C}$, áp suất từ 0,05 MPa đến 5,0 MPa, không áp dụng cho vùng siêu tới hạn cực vi lân cận điểm kỳ dị critical point ($|T - T_c| < 0,5\text{ K}$).

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng định lượng (Positivism/Post-positivism Paradigm), kết hợp chặt chẽ giữa giải tích mô hình hóa toán học số trị và thực nghiệm vật lý kiểm chứng đa cấp độ (Multi-level validation). Thiết kế bao gồm 3 cấp độ:

  • Level 1 (Molecular Level): Tối ưu hóa hàm thế phân tử trên nền tảng phương trình trạng thái Helmholtz.
  • Level 2 (Thermodynamic Cycle Level): Mô phỏng chu trình nén hơi lý thuyết và bán lý tưởng trên 9 môi chất lạnh với các biến số quét liên tục ($T_0 \in [-5^\circ\text{C}, 15^\circ\text{C}]$, $T_k \in [35^\circ\text{C}, 55^\circ\text{C}]$, $\eta \in [0,5; 1,0]$).
  • Level 3 (System/Empirical Level): Thử nghiệm trên trạm thực nghiệm Chiller làm lạnh nước giải nhiệt gió hoàn chỉnh tại phòng thí nghiệm Viện Khoa học & Công nghệ Nhiệt - Lạnh, Đại học Bách khoa Hà Nội.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ các bước nghiêm ngặt:

  1. Thuật toán tối ưu hóa tham số đơn hình Nelder-Mead: Xây dựng hàm mục tiêu sai số tổng hợp phi tuyến:

$$\Phi(T_0, \rho_0, \alpha, Q^{*2}) = \sum_{j=1}^{N_p} \left[ \frac{p_{s,\exp}(T_j) - p_{s,\text{cal}}(T_j)}{p_{s,\exp}(T_j)} \right]^2 + \sum_{k=1}^{N_\rho} \left[ \frac{\rho'{\exp}(T_k) - \rho'{\text{cal}}(T_k)}{\rho'_{\exp}(T_k)} \right]^2 \to \min$$

Thuật toán thực hiện co, dãn, phản xạ và thu hẹp đa diện đơn hình trong không gian 4 chiều cho đến khi chuẩn hội tụ sai số $\epsilon < 10^{-6}$.

QUY TRÌNH THUẬT TOÁN ĐƠN HÌNH TỐI ƯU HÓA THAM SỐ PHÂN TỬ
                 ĐÚNG                               KHÔNG
  1. Giao thức thu thập và hiệu chuẩn thực nghiệm: Hệ thống Chiller được trang bị cụm cảm biến áp suất dải 0–2,5 MPa và cảm biến nhiệt độ Pt100 chuẩn 4 dây.
  2. Đánh giá độ không đảm bảo đo (Measurement Uncertainty): Hiệu chuẩn các kênh đo theo tiêu chuẩn ISO/IEC Guide 98-3 (GUM). Độ không đảm bảo đo mở rộng với hệ số phủ $k=2$ (độ tin cậy 95%) đạt:
    • Cảm biến nhiệt độ: $U_{95}(T) = \pm 0,05\text{ K}$.
    • Cảm biến áp suất: $U_{95}(p) = \pm 0,15\text{ kPa}$.
    • Lưu lượng kế khối lượng: $U_{95}(\dot{m}) = \pm 0,2%$.

Data và phân tích

Kết quả tối ưu hóa đã xác lập thành công bộ 4 thông số đặc trưng phân tử độc bản:

  • Môi chất R1234ze(Z): $T_0 = 413,83\text{ K}$; $\rho_0 = 4,057\text{ mol/l}$; $\alpha = 1,44$; $Q^{*2} = 2,756$.
  • Môi chất R1243zf: $T_0 = 326,83\text{ K}$; $\rho_0 = 4,375\text{ mol/l}$; $\alpha = 1,39$; $Q^{*2} = 1,136$.

Kiểm định độ chuẩn xác (Robustness Check) đối chiếu với các bộ dữ liệu thực nghiệm độc lập của Higashi (2015) và Yukihiro (2011) cho thấy:

  • Sai lệch áp suất bão hòa tính toán ($p_{s,\text{cal}}$) so với thực nghiệm ($p_{s,\exp}$) của R1234ze(Z) đạt độ lệch tuyệt đối trung bình $AAD = 0,08%$, sai số lớn nhất không vượt quá 0,35% trên toàn bộ dải từ 280 K đến 420 K.
  • Sai lệch khối lượng riêng lỏng bão hòa ($\rho'_{s}$) của R1243zf có $AAD = 0,12%$ so với số liệu của Nicola (2014).
  • Bộ bảng dữ liệu nhiệt động vùng bão hòa, vùng lỏng chưa sôi và hơi quá nhiệt cùng đồ thị logarit áp suất - entanpy ($\lg p - h$) được thành lập với độ mịn bước nhiệt độ $\Delta T = 1\text{ K}$, bước áp suất $\Delta p = 0,01\text{ bar}$.

Phát hiện đột phá và implications

Những phát hiện then chốt

HIỆU QUẢ NĂNG LƯỢNG VÀ THAY ĐỔI NHIỆT ĐỘNG CỦA CÁC MCL THAY THẾ
  1. Phát hiện 1 (Khẳng định R1243zf là ứng viên thay thế toàn diện nhất): Trong dải nhiệt độ bay hơi từ $-5^\circ\text{C}$ đến $15^\circ\text{C}$, R1243zf chứng minh khả năng thay thế trực tiếp tối ưu cho cả R134a và R22. So với R134a, COP của R1243zf chỉ thấp hơn $0,52%$, nhưng có ưu điểm vượt trội: nhiệt độ đầu đẩy máy nén ($T_d$) thấp hơn $1,33%$ và tỷ số nén ($\Pi$) thấp hơn $5,82%$, giúp giảm tải cơ học và tăng tuổi thọ máy nén. Khi so sánh với R22, COP của R1243zf chỉ thấp hơn $2,02%$, $T_d$ thấp hơn $5,92%$.
  2. Phát hiện 2 (Khả năng ứng dụng thực tế xuất sắc của R1234ze(E)): Môi chất R1234ze(E) có COP lý thuyết chỉ thấp hơn R134a $0,83%$, tỷ số nén cao hơn $2,48%$, nhiệt độ đầu đẩy thấp hơn $2,06%$. Trong thử nghiệm hệ thống Chiller thực tế, COP trung bình của R1234ze(E) chỉ sai khác dưới $7,8%$ so với R134a mà hoàn toàn không cần thay đổi kết cấu máy nén hay dàn trao đổi nhiệt.
  3. Phát hiện 3 (Tác động phi tuyến mang tính quyết định của hiệu suất nén $\eta$): Khi hiệu suất nén đẳng entropy lý tưởng $\eta = 100%$, COP chu trình đạt từ 5,65 đến 6,02. Tuy nhiên, khi đưa giá trị $\eta$ đo đạc thực tế từ máy nén thực nghiệm (dao động từ $53%$ đến $63%$), COP sụt giảm mạnh xuống dải $3,01$ đến $3,80$ (giảm tương ứng 43,2% đến 46,7%). Điều này chứng minh việc bỏ qua hiệu suất nén thực nghiệm trong các nghiên cứu lý thuyết trước đây dẫn đến sai số dự báo hiệu quả năng lượng rất lớn.
  4. Phát hiện 4 (Tổn thất áp suất thủy lực gây suy giảm COP nghiêm trọng): Luận án định lượng chính xác: độ suy giảm áp suất $\Delta p = 25\text{ kPa}$ trên đường đẩy (từ máy nén đến trước van tiết lưu) và $\Delta p = 63\text{ kPa}$ trên đường hút (từ thiết bị bay hơi về đầu hút máy nén) làm COP thực tế của R1234ze(E) giảm từ dải $3,01–3,80$ xuống chỉ còn $2,33–3,01$ (mức sụt giảm thêm 20,8% đến 22,6%).
  5. Phát hiện 5 (Độ nhạy vượt trội của R1234ze(E) theo nhiệt độ nước lạnh): "Khi nhiệt độ nước lạnh tăng từ 8°C đến 12°C, mức tăng COP tương ứng của chu trình sử dụng môi chất R134a và R1234ze(E) lần lượt là 15% và 19%". R1234ze(E) có độ nhạy nhiệt động cao hơn, mang lại hiệu quả tiết kiệm điện vượt bậc khi vận hành ở dải nhiệt độ nước lạnh cao (chế độ điều hòa tiện nghi).
BIẾN THIÊN COP THEO HIỆU SUẤT NÉN VÀ TỔN THẤT ÁP SUẤT CỦA R1234ze(E)

Implications đa chiều

  • Học thuật & Lý thuyết: Cung cấp phương pháp luận chuẩn hóa kết hợp tối ưu hóa vi mô đa biến Nelder-Mead với phương trình trạng thái Helmholtz đa cực, mở đường cho việc xây dựng cơ sở dữ liệu số cho hàng trăm dẫn xuất fluor hóa chưa từng được đo đạc.
  • Phương pháp luận: Khẳng định quy trình đánh giá MCL bắt buộc phải thực hiện phép thử 3 bước: Mô hình hóa phân tử $\to$ Tính toán chu trình nhiệt động $\to$ Hiệu chỉnh thực nghiệm tổn thất áp suất và $\eta$, chấm dứt việc kết luận vội vã chỉ dựa trên COP lý thuyết thuần túy.
  • Thực tiễn kỹ thuật: Cung cấp hướng dẫn kỹ thuật thay thế trực tiếp R1234ze(E) và R1243zf vào các hệ thống Chiller và ĐHKK đang sử dụng R134a/R22 mà không cần đầu tư tái thiết kế cơ khí, chỉ cần tối ưu lại lượng nạp môi chất và cài đặt van tiết lưu.
  • Chính sách môi trường: Cung cấp luận cứ khoa học thực chứng để Bộ Tài nguyên và Môi trường, Bộ Công Thương xây dựng lộ trình chi tiết thực hiện Tiêu chuẩn TCVN về môi chất lạnh xanh, đẩy nhanh việc thực hiện Bản sửa đổi Kigali tại Việt Nam.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu:

  1. Đặc tính dễ cháy nhẹ (A2L): Luận án tập trung vào tính chất nhiệt động và hiệu quả năng lượng, chưa đi sâu thực nghiệm đặc tính cháy nổ, vận tốc lan truyền ngọn lửa và giới hạn bắt cháy dưới (LFL) của R1243zf và R1234ze(E) trong không gian kín.
  2. Độ bền hóa học và tính tương thích dầu bôi trơn lâu dài: Nghiên cứu thực nghiệm chu trình tiến hành trong các chu kỳ đo ổn định ngắn hạn, chưa khảo sát độ thoái hóa của dầu tổng hợp POE (Polyolester) và PAG (Polyalkylene Glycol) sau 5.000–10.000 giờ vận hành liên tục cùng HFO.
  3. Hiện tượng trượt nhiệt độ và hệ số truyền nhiệt đối lưu: Chưa nghiên cứu sâu sự biến thiên cục bộ của hệ số truyền nhiệt sôi dòng chảy (flow boiling heat transfer coefficient) bên trong ống trao đổi nhiệt có cánh micro-fin.

Chương trình nghiên cứu tương lai (Future Agenda) gồm 5 hướng cụ thể:

  1. Khảo sát thực nghiệm tính chất hòa trộn và độ nhớt động lực học của hỗn hợp HFO/dầu POE ở dải nhiệt độ từ $-20^\circ\text{C}$ đến $120^\circ\text{C}$.
  2. Thiết lập mô hình mô phỏng truyền nhiệt hai pha sôi và ngưng tụ của R1234ze(Z) trong các chu trình Rankine hữu cơ (ORC) thu hồi nhiệt thải công nghiệp nhiệt độ cao ($>90^\circ\text{C}$).
  3. Nghiên cứu hỗn hợp nhị phân và tam phân giữa HFO và HFC (ví dụ R1234yf/R32 hoặc R1234ze(E)/R152a) nhằm triệt tiêu tính cháy mà vẫn duy trì GWP < 150.
  4. Mở rộng trạm thực nghiệm đa năng để kiểm thử khả năng thay thế drop-in của R1243zf trên hệ thống điều hòa ô tô và kho lạnh bảo quản thực phẩm.
  5. Ứng dụng kỹ thuật trí tuệ nhân tạo (Machine Learning/Deep Neural Networks) để dự đoán trực tiếp các tham số tương tác phân tử Helmholtz từ cấu trúc không gian 3D của phân tử.

Tác động và ảnh hưởng

  • Tác động học thuật: Bộ dữ liệu nhiệt động của R1243zf được công bố trên tạp chí khoa học quốc tế uy tín trong danh mục ISI, trở thành tài liệu tham khảo chuẩn cho cộng đồng kỹ thuật nhiệt quốc tế. Ước tính công trình sẽ thu hút hàng chục trích dẫn từ các nhóm nghiên cứu về thermodynamic property modeling và green refrigeration.
  • Tác động ngành công nghiệp HVAC&R: Tạo tiền đề kỹ thuật trực tiếp cho các nhà sản xuất máy lạnh, điều hòa không khí tại Việt Nam (như Nagakawa, Funiki, Reetech) và các doanh nghiệp dịch vụ bảo dưỡng chuyển đổi môi chất lạnh, giảm phụ thuộc vào các HFC đang bị hạn chế nhập khẩu.
  • Tác động môi trường & xã hội: Nếu chuyển đổi thành công 1.000 tấn HCFC-22 và hàng ngàn tấn R134a sang R1234ze(E) và R1243zf, lượng phát thải khí nhà kính trực tiếp sẽ giảm hàng triệu tấn $CO_2$ tương đương mỗi năm (do GWP giảm từ 1.430–1.810 xuống dưới 10), đóng góp trực tiếp vào cam kết Net Zero carbon của Việt Nam tại COP26.

Đối tượng hưởng lợi

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận được phương pháp luận mô hình hóa phân tử Helmholtz và thuật toán tối ưu hóa phi tuyến ứng dụng trong kỹ thuật nhiệt.
  • Giảng viên & Chuyên gia Kỹ thuật Nhiệt: Sở hữu bộ giáo trình, bảng số liệu nhiệt động tra cứu chuẩn xác ($\lg p - h$, $T - s$) của R1234ze(Z) và R1243zf phục vụ công tác đào tạo đại học và sau đại học.
  • Kỹ sư R&D và Thiết kế hệ thống HVAC: Nắm vững các công thức hiệu chỉnh hiệu suất nén $\eta$ và tổn thất áp suất thực tế để tính toán thiết kế kích thước dàn bay hơi, dàn ngưng tụ và chọn công suất máy nén chính xác.
  • Cơ quan quản lý nhà nước (Cục Biến đổi khí hậu, Cục Quản lý Môi trường): Có cơ sở dữ liệu kỹ thuật tin cậy để thẩm định các dự án chuyển đổi công nghệ làm lạnh xanh được tài trợ bởi Quỹ Đa phương thi hành Nghị định thư Montreal.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập thành công mô hình phương trình trạng thái dựa trên năng lượng tự do Helmholtz đa cực (MBEOS) kết hợp thuật toán tối ưu hóa đa biến Nelder-Mead để giải mã chính xác 4 thông số tương tác phân tử vi mô $(T_0, \rho_0, \alpha, Q^{*2})$ của R1234ze(Z) và R1243zf. Nghiên cứu đã mở rộng Lý thuyết Trường Thế lực Thống kê Phân tử (Statistical Mechanics Field Theory) của Müller và Gubbins sang nhóm chất hydrofluoro-olefin có liên kết đôi bất đối xứng, chứng minh rằng tương tác tứ cực $Q^{*2}$ đóng vai trò quyết định trong việc hiệu chỉnh độ cong của đường bão hòa mà các phương trình Cubic EOS kinh điển không thể mô tả được.

2. Đột phá phương pháp luận so với ít nhất 2 nghiên cứu quốc tế trước đây?

So với nghiên cứu của Akasaka (2010) (sử dụng phương trình Patel-Teja gây sai lệch nhiệt dung riêng lên tới 20%) và Sad Jarall (2014) (chỉ khảo sát chu trình lý thuyết đơn thuần của R1234yf), phương pháp luận của luận án mang tính đột phá toàn diện: (1) Giảm sai số xác định thông số bão hòa xuống $AAD < 0,08%$ thông qua mô hình hóa vi mô; (2) Thiết lập quy trình đánh giá liên hoàn tích hợp đồng thời cả tổn thất trở lực dòng chảy ($\Delta P$) và sự suy giảm hiệu suất nén đẳng entropy thực nghiệm ($\eta = 53%–63%$), loại bỏ hoàn toàn các kết luận thiên lệch từ mô hình chu trình lý tưởng.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu hỗ trợ?

Phát hiện bất ngờ nhất là môi chất R1243zf – vốn ít được chú ý hơn so với R1234yf trên thị trường quốc tế – lại thể hiện tính năng nhiệt động vượt trội toàn diện khi thay thế cho cả R134a và R22. Dữ liệu chứng minh: COP của R1243zf chỉ thấp hơn $0,52%$ so với R134a, nhưng có nhiệt độ đầu đẩy thấp hơn $1,33%$ và tỷ số nén thấp hơn $5,82%$. Đối với R22, R1243zf có COP chỉ kém $2,02%$ trong khi nhiệt độ khí xả thấp hơn tới $5,92%$, giải quyết triệt để nguy cơ quá nhiệt dầu bôi trơn máy nén thường gặp ở vùng nhiệt đới.

4. Giao thức lặp lại nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp giao thức lặp lại chi tiết và minh bạch 100%: (1) Toàn bộ hệ phương trình toán học giải tích xác định $F^0, F^H, F^A, F^Q$; (2) Mã nguồn thuật toán tối ưu hóa đơn hình với các điều kiện dừng sai số $\Phi < 10^{-6}$; (3) Bảng thông số chi tiết hệ thống thực nghiệm Chiller, bao gồm kích thước dàn ống trao đổi nhiệt, đặc tính máy nén, sơ đồ bố trí cảm biến đo lường và bảng phân tích độ không đảm bảo đo $U_{95}$ cho từng kênh thu thập dữ liệu.

5. Kế hoạch nghiên cứu 10 năm (10-year Research Agenda) được phác thảo ra sao?

Kế hoạch nghiên cứu 10 năm được định hình theo 3 giai đoạn:

  • Giai đoạn 1 (Năm 1–3): Khảo sát toàn diện hệ số truyền nhiệt đối lưu hai pha sôi/ngưng tụ và độ hòa tan của R1243zf với các gốc dầu bôi trơn tổng hợp mới.
  • Giai đoạn 2 (Năm 4–6): Nghiên cứu ứng dụng R1234ze(Z) cho chu trình bơm nhiệt công nghiệp nhiệt độ cao ($>100^\circ\text{C}$) và chu trình nhiệt điện Rankine hữu cơ (ORC) tận dụng nhiệt thải mặt trời và địa nhiệt.
  • Giai đoạn 3 (Năm 7–10): Phát triển phần mềm thương mại tính toán nhiệt động học tự động cho các hỗn hợp HFO đa cấu tử dựa trên nền tảng trí tuệ nhân tạo và vật lý số lượng tử.

Kết luận

  1. Xây dựng thành công phương pháp và quy trình xác định thông số nhiệt động của môi chất lạnh theo lý thuyết năng lượng tự do Helmholtz tương tác phân tử, giải quyết triệt để bài toán thiếu hụt dữ liệu nhiệt động của môi chất thế hệ mới.
  2. Trích xuất chính xác bộ 4 đại lượng đặc trưng phân tử độc bản cho hai môi chất tiềm năng: R1234ze(Z) ($T_0 = 413,83\text{ K}; \rho_0 = 4,057\text{ mol/l}; \alpha = 1,44; Q^{*2} = 2,756$) và R1243zf ($T_0 = 326,83\text{ K}; \rho_0 = 4,375\text{ mol/l}; \alpha = 1,39; Q^{*2} = 1,136$).
  3. Thiết lập bộ dữ liệu nhiệt động hoàn chỉnh, có độ chính xác cao ($AAD < 0,08%$) trên toàn bộ các vùng lỏng chưa sôi, đường bão hòa, hơi quá nhiệt và đồ thị $\lg p - h$ cho R1234ze(Z) và R1243zf, được công nhận quốc tế qua công bố ISI.
  4. Chứng minh khoa học R1243zf và R1234ze(E) là hai ứng viên thay thế trực tiếp hoàn hảo nhất cho R134a và R22 trong lĩnh vực ĐHKK với độ suy giảm COP lý thuyết cực nhỏ (dưới $0,83%$ với R134a và dưới $2,02%$ với R22) cùng ưu điểm giảm nhiệt độ đầu đẩy máy nén.
  5. Chế tạo thành công hệ thống thực nghiệm Chiller giải nhiệt gió, định lượng chính xác ảnh hưởng của hiệu suất nén thực tế ($\eta = 53%–63%$) và tổn thất áp suất ($\Delta p = 25\text{ kPa}$ và $63\text{ kPa}$) làm COP thực tế dao động từ 2,33 đến 3,01; đồng thời xác nhận COP của R1234ze(E) chỉ sai khác dưới $7,8%$ so với R134a và có tốc độ tăng COP vượt trội (19% so với 15%) khi nhiệt độ nước lạnh nâng từ 8°C lên 12°C.
  6. Mở ra 3 hướng nghiên cứu mới mang tính liên ngành: Ứng dụng HFO trong bơm nhiệt thu hồi nhiệt độ cao, tối ưu hóa hỗn hợp môi chất lạnh giảm phát thải carbon, và phát triển vật liệu trao đổi nhiệt tăng cường dành riêng cho môi chất thế hệ thứ tư.