Tổng quan về luận án

Biến đổi khí hậu toàn cầu đang là thách thức sinh thái nghiêm trọng nhất mà nhân loại phải đối mặt trong thế kỷ 21. Báo cáo đánh giá khoa học lần thứ năm của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC 2014) chỉ rõ các khí nhà kính gốc Flo (khí F) chiếm khoảng 2% tổng nguyên nhân gây ra hiện tượng nóng lên toàn cầu. Đáng chú ý, nhóm Hydrofluorocarbon (HFCs) chiếm tới 91% tổng lượng phát thải khí F, trong đó lĩnh vực kỹ thuật lạnh và điều hòa không khí (ĐHKK) tiêu thụ đến 86% tổng sản lượng HFCs toàn cầu (tính theo đương lượng $CO_2$). Sự chuyển dịch từ các chất làm suy giảm tầng ô-dôn như Chlorofluorocarbon (CFCs) và Hydrochlorofluorocarbon (HCFCs) sang HFCs đã giải quyết được chỉ số làm suy giảm tầng ô-dôn ($\text{ODP} = 0$), nhưng lại làm trầm trọng thêm chỉ số tiềm năng làm nóng toàn cầu ($\text{GWP}$). Điển hình, môi chất R134a có $\text{GWP} = 1430$ và R125 (thành phần chính trong hỗn hợp R410A) có $\text{GWP} = 3500$. Trước sức ép từ Bản sửa đổi, bổ sung Kigali thuộc Nghị định thư Montreal và Quyết định số 1888/QĐ-BTNMT của Bộ Tài nguyên và Môi trường Việt Nam nhằm cắt giảm 35% mức tiêu thụ cơ sở các chất HCFC giai đoạn 2020–2024 (tương đương loại bỏ hoàn toàn 1.000 tấn R22), việc nghiên cứu môi chất lạnh thế hệ thứ tư mang tính bền vững là một yêu cầu khoa học cấp thiết.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: Dòng môi chất mới Hydrofluoro-olefin (HFO) sở hữu $\text{ODP} = 0$ và $\text{GWP} < 1$, nhưng cơ sở dữ liệu nhiệt động học liên tục (nhiệt động trạng thái $p-v-T$, entanpy $h$, entropy $s$) của nhiều đồng phân tiềm năng như R1234ze(Z) và R1243zf còn rất rời rạc, thiếu vắng hoặc có độ sai lệch cao khi áp dụng các phương trình trạng thái cổ điển. Đồng thời, việc đánh giá hiệu quả chuyển hóa năng lượng thực nghiệm trên hệ thống máy lạnh thực tế của các môi chất HFO so với môi chất truyền thống (R134a, R22) trong điều kiện khí hậu nhiệt đới tại Việt Nam chưa được kiểm chứng đầy đủ.

Luận án tập trung giải quyết ba câu hỏi nghiên cứu (Research Questions) và kiểm chứng ba giả thuyết khoa học (Hypotheses) tương ứng:

  • RQ1: Làm thế nào để xây dựng mô hình toán học chính xác xác định liên tục toàn bộ các thông số nhiệt động của môi chất HFO từ dữ liệu thực nghiệm giới hạn?
    H1: Phương trình trạng thái dựa trên cấu trúc phân tử (MBEOS) kết hợp mô hình năng lượng tự do Helmholtz có thể mô tả chính xác trạng thái nhiệt động của R1234ze(Z) và R1243zf với sai số trung bình tuyệt đối ($\text{AAD}$) dưới 1%.
  • RQ2: Khả năng biến đổi năng lượng của các môi chất HFO (R1234ze(E), R1234ze(Z), R1243zf, R1234yf) thay đổi như thế nào trong các chu trình ĐHKK lý thuyết so với R134a và R22?
    H2: Các môi chất HFO có thể đạt hệ số hiệu quả năng lượng ($\text{COP}$) tương đương hoặc vượt trội so với môi chất truyền thống khi tối ưu hóa thông số vận hành.
  • RQ3: Mức độ tương thích và hiệu quả thực tế của R1234ze(E) khi thay thế trực tiếp ("drop-in replacement") cho R134a trên hệ thống máy làm lạnh nước (water chiller) thực nghiệm là gì?
    H3: R1234ze(E) vận hành ổn định trên hệ thống hiện hữu và duy trì $\text{COP}$ cạnh tranh mà không đòi hỏi thay đổi cấu trúc cơ khí căn bản.

Khung lý thuyết của công trình tích hợp phương trình trạng thái dựa trên lực tương tác phân tử (MBEOS), nhiệt động học thống kê lượng tử và mô hình chu trình nhiệt động kỹ thuật. Tác động của nghiên cứu thể hiện qua việc công bố bộ thông số nhiệt động học toàn diện cho 2 môi chất HFO, đánh giá chu trình lý thuyết trên dải nhiệt độ bay hơi từ $-5^\circ\text{C}$ đến $15^\circ\text{C}$, nhiệt độ ngưng tụ từ $30^\circ\text{C}$ đến $55^\circ\text{C}$, và kiểm chứng thực nghiệm hoàn chỉnh trên mô hình máy làm lạnh nước với hệ thống thu thập dữ liệu tự động đạt độ không đảm bảo đo $U_{95}$ theo tiêu chuẩn quốc gia.

Literature Review và Positioning

Lịch sử phát triển môi chất lạnh trải qua 4 thế hệ: thế hệ 1 (các chất tự nhiên như $CO_2, NH_3, SO_2$), thế hệ 2 (CFCs như R11, R12 với tính năng nhiệt động ưu việt nhưng phá hủy tầng ô-dôn cực mạnh), thế hệ 3 (HCFCs như R22 và HFCs như R134a, R404A, R410A giải quyết vấn đề ô-dôn nhưng gây hiệu ứng nhà kính nặng nề), và thế hệ 4 (HFOs và chất tự nhiên cải tiến). Dòng HFO là các hydrocacbon không no (olefin) được flo hóa, nổi bật là các đồng phân của propen gồm R1225 (4 đồng phân), R1234 (7 đồng phân), và R1243 (8 đồng phân).

Trong bức tranh tổng quan học thuật, tồn tại hai luồng quan điểm và phương pháp tiếp cận chính để xác định tính chất nhiệt động của môi chất:

                  ┌──────────────────────────────────────────────────────────┐
                  │   CÁC PHƯƠNG PHÁP XÁC ĐỊNH TÍNH CHẤT NHIỆT ĐỘNG HFO      │
                  └────────────────────────────┬─────────────────────────────┘
                                               │
               ┌───────────────────────────────┴───────────────────────────────┐
               ▼                                                               ▼
┌─────────────────────────────┐                                 ┌─────────────────────────────┐
│  PHƯƠNG TRÌNH DẠNG BẬC BA   │                                 │    PHƯƠNG TRÌNH MBEOS /     │
│   (Peng-Robinson, Patel-    │                                 │    PC-SAFT NĂNG LƯỢNG TỰ DO │
│          Teja)              │                                 │           HELMHOLTZ         │
├─────────────────────────────┤                                 ├─────────────────────────────┤
│ • Đơn giản, ít tham số.     │                                 │ • Dựa trên cơ học thống kê. │
│ • Sai số vùng lỏng: 5-10%.  │                                 │ • Tách biệt lực phân tử:    │
│ • Sai số nhiệt dung Cp:     │                                 │   F = F⁰ + Fᴴ + Fᴬ + Fᑫ     │
│   lên tới 20% (Akasaka).    │                                 │ • Sai số cực nhỏ: < 0.5-1%. │
│ • Không mô tả tốt cấu trúc  │                                 │ • Độ tin cậy cao trên toàn  │
│   bất đối xứng phân tử.     │                                 │   dải trạng thái lỏng - hơi.│
└─────────────────────────────┘                                 └─────────────────────────────┘
  1. Nhóm phương trình trạng thái dạng bậc ba cải tiến (Cubic EOS): Điển hình là nghiên cứu của Jacobsen (1992) trên Peng-Robinson cho R134a với sai số trung bình tuyệt đối ($\text{AAD}$) trong vùng lỏng từ 5% đến 10%, vùng hơi từ 2% đến 5%. Brown (2009, 2010) áp dụng phương trình Peng-Robinson để ước lượng tính chất nhiệt động của 8 đồng phân HFO, nhưng chỉ ra sai lệch áp suất pha hơi của R1234yf so với thực nghiệm luôn lớn hơn 1% và lên tới 4% ở dải nhiệt độ thấp (240–260 K). Nghiên cứu của Akasaka (2010) sử dụng phương trình Patel-Teja cho R1234yf cho thấy sai lệch khối lượng riêng pha lỏng là 2.5%, và sai lệch nhiệt dung riêng đẳng áp ($C_p$) lên tới 20%, cho thấy giới hạn rõ rệt của phương trình bậc ba trong việc mô tả cấu trúc phân tử bất đối xứng.
  2. Nhóm phương trình trạng thái đa tham số và phương trình tương tác phân tử (MBEOS / PC-SAFT): Các công trình của Wagner (2002), Lemmon (2009 cho Propan; 2012 cho Dimethyl Carbonate) đạt độ chính xác rất cao nhưng đòi hỏi hàng trăm điểm thực nghiệm phân bố đều trên mọi vùng trạng thái. Khi dữ liệu thực nghiệm khan hiếm, lý thuyết tương tác phân tử dựa trên năng lượng tự do Helmholtz của Chapman, Gross, Sadowski (PC-SAFT) và Muller (xử lý mômen đa cực) là hướng đi đột phá.

Về mặt thực nghiệm, các công bố quốc tế chủ yếu tập trung vào R1234yf và R1234ze(E). Katsuyuki (2010) đo khối lượng riêng bão hòa của R1234ze(E) và R1234ze(Z) bằng phương pháp chiết tách (extraction method). An (2011) công bố số liệu $p-v-T$ của R1234ze(E) pha khí bằng phương pháp Burnett với độ không đảm bảo đo nhiệt độ 2.3 mK và áp suất 600 Pa. Higashi (2010) xác định thông số tới hạn của R1234ze(Z) ($T_c = 423.27\text{ K}$, $p_c = 3.533\text{ MPa}$, $\rho_c = 473\text{ kg/m}^3$). Về chu trình ứng dụng, Sad Jarall (2012) thử nghiệm R1234yf trên máy nén rô-to lăn; Claudio Zilio (2011) thử nghiệm trên hệ thống ĐHKK ô tô; Barbara Minor (2010) khảo sát máy làm lạnh đồ uống. Tại Việt Nam, các tác giả như Võ Chí Chính, Lê Xuân Hòa, Nguyễn Đức Lợi chủ yếu khảo sát chu trình lý thuyết trên các môi chất HCFC/HFC hoặc hỗn hợp hydrocacbon (R290/R600a), chưa có công trình nào xây dựng hệ thống dữ liệu nhiệt động hoàn chỉnh kết hợp thử nghiệm máy làm lạnh nước chuyên dụng cho dòng HFO.

Luận án định vị chính xác vào giao điểm học thuật: Ứng dụng mô hình tương tác phân tử năng lượng tự do Helmholtz để thiết lập bộ dữ liệu nhiệt động đầy đủ cho R1234ze(Z) và R1243zf, đồng thời thiết lập chu trình thực nghiệm máy làm lạnh nước kiểm chứng khả năng thay thế của R1234ze(E) trong điều kiện vận hành tại Việt 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 thêm lý thuyết nhiệt động học thống kê vi mô áp dụng cho các chất lưu hữu cơ có cấu trúc phân tử bất đối xứng cao:

  1. Hoàn thiện mô hình năng lượng tự do Helmholtz mở rộng: Tách biệt chính xác các thành phần tương tác năng lượng vi mô của phân tử HFO theo biểu thức: $$F = F^0 + F^r = F^0 + (F^H + F^A + F^Q)$$ Trong đó:
    • $F^0$: Năng lượng tự do Helmholtz của khí lý tưởng.
    • $F^H$: Năng lượng do tương tác đẩy của phần lõi cứng (Hard-body contribution).
    • $F^A$: Năng lượng do lực hấp dẫn phân tán (Attractive dispersion force contribution) dựa trên thế Lennard-Jones cải biên.
    • $F^Q$: Năng lượng do tương tác định hướng đa cực (Quadrupolar contribution) phản ánh sự phân bố điện tích phi đối xứng của các liên kết C-F và C=C.
  2. Xác lập tập tham số đặc trưng phân tử vi mô: Xác định chính xác các đại lượng vật lý đặc trưng gồm đường kính hạt rắn tương đương ($\sigma$), chiều sâu giếng thế năng lượng tương tác ($\epsilon/k_B$), mômen tứ cực không thứ nguyên ($Q^{*2}$), và hệ số bất đối xứng phân tử ($\omega$) cho hai chất R1234ze(Z) và R1243zf.
  3. Chứng minh quy luật chuyển pha và quan hệ trạng thái: Thiết lập mô hình giải tích cho phép dẫn xuất chính xác các hàm nhiệt động bậc một và bậc hai ($h, s, C_p, C_v$, tốc độ âm thanh $w$) mà không cần dựa vào các giả thiết bán thực nghiệm thiếu cơ sở vật lý.

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

Khung phân tích của luận án tích hợp ba trụ cột phương pháp luận:

  • Cơ sở tương tác phân tử MBEOS: Mô hình hóa toán học các lực tương tác nội tại giữa các mắt xích phân tử.
  • Giải thuật tối ưu hóa đơn hình (Nelder-Mead Simplex Algorithm): Cực tiểu hóa hàm mục tiêu sai số tuyệt đối giữa giá trị tính toán từ mô hình và các điểm dữ liệu thực nghiệm rời rạc đã công bố: $$\text{AAD} = \frac{1}{N}\sum_{i=1}^N \left|\frac{X_{\text{exp},i} - X_{\text{cal},i}}{X_{\text{exp},i}}\right| \times 100% \to \min$$
  • Mô hình chu trình nhiệt động kỹ thuật hệ thống lạnh: Đánh giá hiệu quả biến đổi năng lượng thông qua hệ số $\text{COP}$, năng suất lạnh thể tích ($q_v$), tỷ số nén ($\Pi$), và hiệu suất nén không thuận nghịch ($\eta_i$).
                      ┌────────────────────────────────────────────────────────┐
                      │    KHUNG PHÂN TÍCH TỔNG HỢP VÀ TỐI ƯU HÓA LUẬN ÁN      │
                      └───────────────────────────┬────────────────────────────┘
                                                  │
             ┌────────────────────────────────────┼────────────────────────────────────┐
             ▼                                    ▼                                    ▼
┌───────────────────────────┐        ┌───────────────────────────┐        ┌───────────────────────────┐
│   CẤU TRÚC PHÂN TỬ VI MÔ  │        │   GIẢI THUẬT NELDER-MEAD  │        │  CHU TRÌNH NHIỆT ĐỘNG HỌC │
├───────────────────────────┤        ├───────────────────────────┤        ├───────────────────────────┤
│ • Thế Lennard-Jones       │        │ • Hàm mục tiêu sai số     │        │ • Mô phỏng chu trình nén  │
│ • Tương tác đa cực Fᑫ     │───────►│   AAD -> min              │───────►│   hơi 1 cấp, hồi nhiệt    │
│ • Phân tách năng lượng    │        │ • Xác định bộ tham số     │        │ • Đánh giá COP, qᵥ, Π     │
│   Helmholtz: F = F⁰ + Fʳ  │        │   tối ưu: σ, ε/k_B, Q*²   │        │ • Kiểm chứng Water Chiller│
└───────────────────────────┘        └───────────────────────────┘        └───────────────────────────┘

Điều kiện biên của khung phân tích: Mô hình áp dụng hiệu quả cho môi chất HFO có mạch carbon ngắn ($C_3$), không liên kết hydro mạnh (non-associating fluids), trong dải áp suất từ áp suất chân không đến cận tới hạn ($p < p_c$) và nhiệt độ từ nhiệt độ ba trạng thái đến nhiệt độ quá nhiệt $100^\circ\text{C}$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ lập trường nhận thức luận thực chứng (Positivism) kết hợp mô hình hóa giải tích số trị và kiểm chứng thực nghiệm vật lý khách quan. Thiết kế nghiên cứu đa cấp độ (multi-level design) bao gồm:

  • Cấp độ vi mô: Tính toán cơ học phân tử và nhiệt động học thống kê.
  • Cấp độ trung gian (Meso-level): Xây dựng cơ sở dữ liệu số trị, lập bảng nhiệt động trạng thái bão hòa và vùng 1 pha (lỏng chưa sôi, hơi quá nhiệt), xây dựng đồ thị $\lg p - h$.
  • Cấp độ vĩ mô: Mô phỏng chu trình nhiệt động và thử nghiệm trên hệ thống thiết bị thực tế.

Quy trình nghiên cứu rigorous

┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                 QUY TRÌNH NGHIÊN CỨU 4 BƯỚC CHẶT CHẼ                                   │
└───────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                    │
 ┌──────────────────────────────────────────────────┴──────────────────────────────────────────────────┐
 │ 1. SÀNG LỌC ĐỒNG PHÂN HFO                                                                           │
 │    19 đồng phân -> Đánh giá ODP=0, GWP<1, dải áp suất phù hợp -> Chọn R1234ze(Z), R1243zf,          │
 │    R1234ze(E), R1234yf.                                                                             │
 └──────────────────────────────────────────────────┬──────────────────────────────────────────────────┘
                                                    │
 ┌──────────────────────────────────────────────────┴──────────────────────────────────────────────────┐
 │ 2. XÁC ĐỊNH THÔNG SỐ NHIỆT ĐỘNG MBEOS                                                               │
 │    Thu thập dữ liệu thực nghiệm cơ sở -> Tối ưu hóa Simplex -> Xác lập bảng dữ liệu và đồ thị lg p-h│
 │    cho R1234ze(Z) và R1243zf.                                                                       │
 └──────────────────────────────────────────────────┬──────────────────────────────────────────────────┘
                                                    │
 ┌──────────────────────────────────────────────────┴──────────────────────────────────────────────────┐
 │ 3. KHẢO SÁT CHU TRÌNH LÝ THUYẾT ĐA BIẾN                                                             │
 │    Biến thiên t₀ (-5 đến 15°C), tₖ (30 đến 55°C), ηᵢ (0.6 đến 1.0) -> Đánh giá so sánh với R134a,  │
 │    R22, R245fa.                                                                                     │
 └──────────────────────────────────────────────────┬──────────────────────────────────────────────────┘
                                                    │
 ┌──────────────────────────────────────────────────┴──────────────────────────────────────────────────┐
 │ 4. THỰC NGHIỆM ĐỐI CHỨNG WATER CHILLER                                                              │
 │    Chế tạo hệ thống thử nghiệm -> Hiệu chuẩn đầu đo cấp chính xác cao tại Viện Đo lường Việt Nam    │
 │    -> Vận hành đối chứng R1234ze(E) vs R134a.                                                       │
└────────────────────────────────────────────────────────────────────────────────────────────────────────┘

Công tác kiểm chuẩn thiết bị đo được thực hiện nghiêm ngặt tại Viện Đo lường Việt Nam:

  • 12 cảm biến nhiệt độ Pt100 bọc kim loại được hiệu chuẩn trong dải $-20^\circ\text{C}$ đến $100^\circ\text{C}$, đạt độ không đảm bảo đo mở rộng $U_{95} = 0.05\text{ K} - 0.08\text{ K}$.
  • Các cảm biến áp suất dải $0 - 30\text{ bar}$ đạt cấp chính xác 0.1, $U_{95} = 0.015\text{ bar}$.
  • Lưu lượng kế đo lưu lượng nước lạnh và lưu lượng môi chất đạt sai số tương đối $< 0.2%$.
  • Dữ liệu được thu thập đồng bộ qua module chuyển đổi tín hiệu chuyên dụng kết nối máy tính thời gian thực (chu kỳ quét 1 giây/lần, kiểm tra độ ổn định trạng thái xác lập trong 10–15 phút).

Data và phân tích

Hệ số tin cậy và sai số của mô hình được kiểm định qua phân tích thống kê so sánh với dữ liệu thực nghiệm của các phòng thí nghiệm hàng đầu thế giới:

Môi chất Đại lượng kiểm định Dữ liệu thực nghiệm đối chứng Sai số trung bình tuyệt đối ($\text{AAD}$) Độ lệch cực đại ($\text{Max Deviation}$)
R1234ze(Z) Áp suất bão hòa $p_s$ Yukihiro (2011), Higashi (2010) 0.32% 0.85%
R1234ze(Z) Khối lượng riêng lỏng $\rho'_s$ Katsuyuki (2010), Yukihiro (2011) 0.48% 1.12%
R1234ze(Z) Thể tích riêng pha hơi $v_g$ Yukihiro (2011) 0.65% 1.40%
R1243zf Áp suất bão hòa $p_s$ Brown (2010) 0.28% 0.72%
R1243zf Khối lượng riêng lỏng $\rho'_s$ Nicola (2012) 0.41% 0.95%
R1243zf Áp suất pha hơi $p_g$ Brown (2010) 0.53% 1.25%

Các kết quả kiểm định độ vững (robustness checks) khẳng định mô hình MBEOS đạt độ chuẩn xác vượt bậc so với các phương trình Peng-Robinson và Patel-Teja truyền thống trên toàn bộ dải thông số làm việc.

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

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

  1. Xác lập thành công bảng số liệu nhiệt động và đồ thị $\lg p - h$ chuẩn xác cho R1234ze(Z) và R1243zf: Lần đầu tiên cung cấp đầy đủ thông số entanpy ($h$), entropy ($s$), khối lượng riêng ($\rho$), thể tích riêng ($v$) trên đường bão hòa và vùng một pha (áp suất lên tới $6\text{ MPa}$, nhiệt độ tới $150^\circ\text{C}$) cho hai môi chất này, phục vụ trực tiếp cho công tác tính toán thiết kế kỹ thuật.
  2. Khám phá tiềm năng vượt trội của R1234ze(Z) cho ứng dụng nhiệt độ cao: R1234ze(Z) sở hữu nhiệt độ sôi tiêu chuẩn $t_b = 9.8^\circ\text{C}$ và thông số tới hạn cao ($T_c = 426.8\text{ K} \approx 153.65^\circ\text{C}$, $p_c = 3.533\text{ MPa}$). Kết quả phân tích chu trình khẳng định R1234ze(Z) là môi chất lý tưởng thay thế R245fa trong các hệ thống bơm nhiệt nhiệt độ cao (nhiệt độ cấp nước nóng lên tới $90–100^\circ\text{C}$) và chu trình Rankine hữu cơ (ORC) thu hồi nhiệt thải công nghiệp với hiệu suất biến đổi năng lượng cao hơn R245fa từ 4.5% đến 8.2%.
  3. Đặc tính tương đương hoàn hảo của R1243zf đối với R134a: R1243zf có đường bão hòa áp suất - nhiệt độ gần như trùng khít với R134a trong dải nhiệt độ ĐHKK ($t_0 = 0 - 10^\circ\text{C}, t_k = 40 - 50^\circ\text{C}$). Tỷ số nén $\Pi$ của R1243zf chỉ chênh lệch $< 1.5%$ so với R134a, năng suất lạnh riêng thể tích $q_v$ đạt xấp xỉ 95% R134a, và $\text{COP}$ lý thuyết tương đương (độ lệch $< 0.8%$), cho thấy đây là ứng viên thay thế "drop-in" không đòi hỏi điều chỉnh thiết kế máy nén.
  4. Hiệu quả thực nghiệm của R1234ze(E) trên hệ thống làm lạnh nước: Kết quả thử nghiệm thực tế trên mô hình máy làm lạnh nước khẳng định:
    • Ở chế độ nhiệt độ nước lạnh ra khỏi dàn bay hơi cao ($t_{w2} = 10 - 15^\circ\text{C}$), $\text{COP}$ thực nghiệm của R1234ze(E) đạt từ $3.25 - 4.10$, cao hơn R134a từ 1.2% đến 3.5%.
    • Ở chế độ nhiệt độ nước lạnh thấp ($t_{w2} = 5 - 7^\circ\text{C}$), $\text{COP}$ của R1234ze(E) thấp hơn R134a từ 2.0% đến 4.8%.
    • Nhiệt độ cuối tầm nén của R1234ze(E) luôn thấp hơn R134a từ $3 - 6\text{ K}$, giúp tăng tuổi thọ dầu bôi trơn và độ tin cậy của máy nén.
    • Năng suất lạnh riêng thể tích $q_v$ của R1234ze(E) thấp hơn R134a khoảng 18–22%, dẫn đến công suất lạnh của hệ thống giảm tương ứng nếu giữ nguyên thể tích quét của máy nén.
                      ┌────────────────────────────────────────────────────────┐
                      │    MA TRẬN ĐÁNH GIÁ TIỀM NĂNG CỦA CÁC MÔI CHẤT HFO     │
                      └───────────────────────────┬────────────────────────────┘
                                                  │
             ┌────────────────────────────────────┼────────────────────────────────────┐
             ▼                                    ▼                                    ▼
┌───────────────────────────┐        ┌───────────────────────────┐        ┌───────────────────────────┐
│       R1234ze(Z)          │        │         R1243zf           │        │       R1234ze(E)          │
├───────────────────────────┤        ├───────────────────────────┤        ├───────────────────────────┤
│ • Tc cao (~153.65°C)      │        │ • Áp suất bão hòa trùng   │        │ • Đã thử nghiệm Water     │
│ • Thay thế R245fa         │        │   khít R134a              │        │   Chiller thực tế         │
│ • Ứng dụng: Bơm nhiệt     │        │ • Tỷ số nén lệch < 1.5%   │        │ • COP cao hơn R134a khi   │
│   nhiệt độ cao, chu trình │        │ • Ứng dụng: Thay thế      │        │   nước lạnh ra > 10°C     │
│   ORC thu hồi nhiệt thải. │        │   "drop-in" cho R134a.    │        │ • Nhiệt độ nén thấp hơn.  │
└───────────────────────────┘        └───────────────────────────┘        └───────────────────────────┘

Implications đa chiều

  • Ý nghĩa khoa học: Cung cấp công cụ toán học và hệ thống tham số tương tác phân tử có độ chính xác cao, mở ra phương pháp luận chuẩn mực để tính toán số liệu nhiệt động cho các môi chất hữu cơ mới mà không cần tốn kém chi phí chế tạo thiết bị đo thực nghiệm phức tạp cho từng điểm trạng thái.
  • Ý nghĩa thực tiễn & công nghiệp: Cung cấp cơ sở dữ liệu tính toán thiết kế cho các nhà sản xuất máy lạnh, điều hòa không khí và bơm nhiệt; đưa ra khuyến nghị kỹ thuật cụ thể: Khi thay thế R1234ze(E) cho R134a trong hệ thống chiller, cần tăng thể tích hút của máy nén khoảng 20% hoặc sử dụng máy nén biến tần để bù đắp năng suất lạnh thể tích, đồng thời tận dụng hiệu quả biến đổi năng lượng vượt trội ở phụ tải bán phần và nhiệt độ bay hơi cao.
  • Ý nghĩa chính sách: Đóng góp luận cứ khoa học thực chứng để các cơ quan quản lý nhà nước (Bộ Tài nguyên và Môi trường, Bộ Công thương) xây dựng tiêu chuẩn kỹ thuật quốc gia và lộ trình thay thế môi chất lạnh theo cam kết quốc tế.

Limitations và Future Research

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

  1. Giới hạn về đặc tính truyền nhiệt và động học dòng chảy: Luận án tập trung vào tính chất nhiệt động học trạng thái ($p-v-T-h-s$) và chu trình nén hơi, chưa đi sâu khảo sát các hệ số vận chuyển (transport properties) như độ dẫn nhiệt ($\lambda$), độ nhớt động lực học ($\mu$), và hệ số truyền nhiệt sôi/ngưng tụ cục bộ bên trong ống trao đổi nhiệt.
  2. Tính bắt cháy nhẹ (A2L Safety Classification): Các môi chất HFO nghiên cứu đều thuộc nhóm an toàn A2L (cháy nhẹ), đòi hỏi các giải pháp thông gió, cảm biến rò rỉ và tiêu chuẩn an toàn phòng chống cháy nổ trong không gian kín chưa được khảo sát chi tiết trong mô hình thực nghiệm.
  3. Quy mô thực nghiệm: Thử nghiệm thực tế mới tiến hành trên hệ thống máy làm lạnh nước (Water Chiller) giải nhiệt nước sử dụng máy nén pittông nửa kín, chưa khảo sát trên máy nén trục vít, máy nén ly tâm công suất lớn hoặc hệ thống ĐHKK trực tiếp (VRV/VRF).

Định hướng nghiên cứu tiếp theo:

  • Mở rộng xác định tính chất vận chuyển ($\lambda, \mu$) và cơ chế truyền nhiệt biến đổi pha của R1234ze(Z) và R1243zf trong ống vi rãnh (microchannel).
  • Nghiên cứu các hệ hỗn hợp môi chất (blends) giữa HFO và HFC (ví dụ R1234yf/R32 hoặc R1234ze(E)/R32) nhằm trung hòa tính cháy, tối ưu hóa năng suất lạnh thể tích $q_v$ mà vẫn duy trì chỉ số $\text{GWP} < 150$.
  • Thử nghiệm độ bền hóa học, độ hòa tan và tương thích lâu dài của môi chất HFO với các loại dầu bôi trơn tổng hợp (POE, PAG, PVE) và vật liệu làm kín (elastomers) trong điều kiện vận hành khắc nghiệt.

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

  • Tác động học thuật: Kết quả luận án đóng góp vào kho tàng nhiệt động học kỹ thuật nguồn dữ liệu số trị có độ tin cậy cao, dự kiến thu hút trích dẫn khoa học lớn trong các nghiên cứu về môi chất thế hệ mới, chu trình ORC và bơm nhiệt công nghiệp.
  • Chuyển đổi công nghiệp: Định hướng cho các tập đoàn sản xuất thiết bị HVAC&R tại Việt Nam và khu vực chuyển đổi công nghệ sản xuất chiller, ĐHKK thương mại và bơm nhiệt sấy công nghiệp thân thiện môi trường.
  • Đóng góp bảo vệ môi trường toàn cầu: Ứng dụng các môi chất HFO với $\text{GWP} < 1$ thay thế R134a ($\text{GWP} = 1430$) và R22 ($\text{ODP} = 0.055, \text{GWP} = 1810$) trực tiếp loại bỏ nguy cơ phá hủy tầng ô-dôn và giảm phát thải hàng triệu tấn $CO_2$ tương đương mỗi năm từ rò rỉ môi chất và tiêu thụ điện năng.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận mẫu mực về ứng dụng MBEOS và giải thuật tối ưu hóa tham số vi mô phân tử; khai thác bảng dữ liệu nhiệt động chi tiết để phát triển các mô hình mô phỏng nâng cao.
  • Kỹ sư R&D ngành Lạnh & ĐHKK: Sở hữu bộ công cụ tra cứu $p-h$ chuẩn xác để tính toán cân bằng nhiệt, lựa chọn máy nén, thiết kế thiết bị trao đổi nhiệt dàn bay hơi/ngưng tụ tối ưu cho dòng HFO.
  • Doanh nghiệp vận hành hệ thống năng lượng: Có cơ sở khoa học để xây dựng phương án chuyển đổi môi chất ("retrofit" hoặc "drop-in") tiết kiệm năng lượng và giảm thiểu chi phí phát thải carbon.
  • Nhà hoạch định chính sách: Có dữ liệu thực chứng khách quan phục vụ việc ban hành quy chuẩn kỹ thuật và kiểm soát nhập khẩu môi chất lạnh theo hạn ngạch quốc gia.

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 độc đáo nhất là việc tích hợp thành công mô hình tương tác phân tử MBEOS dựa trên thế năng Helmholtz mở rộng ($F = F^0 + F^H + F^A + F^Q$) để giải mã chính xác các đặc tính nhiệt động của các phân tử HFO bất đối xứng cao (R1234ze(Z) và R1243zf). Nghiên cứu đã mở rộng lý thuyết cơ học thống kê lượng tử từ các phân tử đối xứng đơn giản sang các hợp chất hydrocacbon không no đa cực chứa liên kết đôi và nhiều nguyên tử Flo, khắc phục triệt để sai số lớn của các phương trình trạng thái dạng bậc ba (Cubic EOS).

2. Tính đổi mới trong phương pháp nghiên cứu so với các công trình quốc tế trước đây?

So với Brown (2010) sử dụng phương trình Peng-Robinson (sai số áp suất hơi tới 4%) và Akasaka (2010) dùng phương trình Patel-Teja (sai số nhiệt dung riêng $C_p$ tới 20%), luận án đã ứng dụng mô hình tương tác phân tử kết hợp thuật toán tối ưu hóa đa biến Nelder-Mead Simplex, giúp hạ thấp sai số áp suất bão hòa xuống 0.28% – 0.32% và sai số khối lượng riêng lỏng xuống 0.41% – 0.48%. Đồng thời, luận án thiết lập chu trình kiểm chứng thực nghiệm hoàn chỉnh có kiểm chuẩn sai số cảm biến đạt độ không đảm bảo đo $U_{95}$ khắt khe tại Viện Đo lường Quốc gia.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và lời giải thích nhiệt động học?

Phát hiện thực nghiệm bất ngờ nhất là môi chất R1234ze(E) có hệ số $\text{COP}$ thực tế cao hơn R134a từ 1.2% đến 3.5% khi nhiệt độ nước lạnh ra dàn bay hơi ở mức cao ($10 - 15^\circ\text{C}$), mặc dù năng suất lạnh thể tích $q_v$ của nó thấp hơn R134a khoảng 20%. Giải thích nhiệt động học: R1234ze(E) có tỷ số nén $\Pi$ thấp hơn và tổn thất tiết lưu nhỏ hơn ở dải nhiệt độ bay hơi cao; đồng thời, nhiệt độ cuối tầm nén thấp hơn giúp giảm tổn thất truyền nhiệt không thuận nghịch tại xi lanh máy nén, nâng cao hiệu suất chỉ thị ($\eta_i$) của chu trình thực tế.

4. Luận án có cung cấp quy trình tái lặp nghiên cứu (replication protocol) không?

Luận án cung cấp quy trình tái lặp đầy đủ và minh bạch: Toàn bộ hệ phương trình toán học xác định các thành phần năng lượng $F^0, F^H, F^A, F^Q$; bảng tham số vi mô phân tử tối ưu ($\sigma, \epsilon/k_B, Q^{*2}, \omega$); lưu đồ giải thuật chi tiết của thuật toán Simplex; sơ đồ nguyên lý hệ thống thực nghiệm, quy cách cảm biến, vị trí lắp đặt đầu đo và quy trình kiểm chuẩn đo lường.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm tập trung vào 3 giai đoạn:

  • Giai đoạn 1 (1–3 năm): Khảo sát hệ số dẫn nhiệt, độ nhớt và mô hình hóa truyền nhiệt sôi/ngưng tụ của R1234ze(Z) và R1243zf trong ống trao đổi nhiệt hiệu năng cao.
  • Giai đoạn 2 (4–6 năm): Nghiên cứu tối ưu hóa hỗn hợp đa cấu tử HFO/HFC/Hydrocarbon hướng tới môi chất có $\text{GWP} < 10$, triệt tiêu tính cháy (nhóm A1) và tối ưu hóa $q_v$.
  • Giai đoạn 3 (7–10 năm): Thiết kế, chế tạo thử nghiệm công nghiệp các dòng máy bơm nhiệt nhiệt độ cao ($120–150^\circ\text{C}$) và hệ thống ORC phát điện mini tận dụng nhiệt thải sử dụng R1234ze(Z).

Kết luận

  1. Đã giải quyết thành công bài toán xác định tính chất nhiệt động học của hai môi chất HFO tiềm năng là R1234ze(Z) và R1243zf bằng phương pháp mô hình hóa tương tác phân tử MBEOS kết hợp giải thuật tối ưu hóa đơn hình Simplex.
  2. Công bố bộ dữ liệu nhiệt động học liên tục, toàn diện và đồ thị $\lg p - h$ đạt độ chính xác cao với sai số trung bình tuyệt đối ($\text{AAD}$) dưới 0.5% đối với áp suất bão hòa và khối lượng riêng so với thực nghiệm quốc tế.
  3. Chứng minh trên phương diện lý thuyết nhiệt động học: R1243zf là chất thay thế "drop-in" tối ưu cho R134a trong điều hòa không khí; R1234ze(Z) là môi chất vượt trội thay thế R245fa trong bơm nhiệt nhiệt độ cao và chu trình ORC.
  4. Chế tạo thành công mô hình thực nghiệm máy làm lạnh nước (Water Chiller) và chứng minh R1234ze(E) có thể thay thế an toàn, hiệu quả cho R134a, đạt $\text{COP}$ vượt trội ở chế độ nhiệt độ bay hơi cao và giảm nhiệt độ cuối tầm nén từ 3–6 K.
  5. Đóng góp luận cứ khoa học thực tiễn vững chắc phục vụ lộ trình quốc gia về giảm phát thải khí nhà kính và loại bỏ môi chất lạnh phá hủy tầng ô-dôn theo Nghị định thư Montreal và Thỏa thuận Kigali.
  6. Mở ra các hướng nghiên cứu liên ngành mới về hỗn hợp môi chất HFO tiên tiến, kỹ thuật trao đổi nhiệt vi kênh và phát triển các hệ thống thiết bị nhiệt lạnh siêu hiệu suất cho nền kinh tế carbon thấp.