Tổng quan về luận án
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 chiếm tỷ trọng từ 30% đến 50% tổng điện năng tiêu thụ trong khu vực dân dụng và thương mại. Trước sức ép kép từ cuộc khủng hoảng năng lượng hóa thạch và cam kết giảm phát thải khí nhà kính toàn cầu theo Nghị định thư Montreal và Thỏa thuận Paris, việc đổi mới công nghệ làm lạnh theo hướng nâng cao hiệu suất nhiệt động và tận dụng các nguồn năng lượng tái tạo đã trở thành một mệnh lệnh nghiên cứu cấp thiết. Trong các giải pháp tiềm năng, chu trình làm lạnh bằng bơm phun nhiệt (ejector refrigeration system) vận hành bằng nguồn nhiệt có nhiệt thế thấp (low-grade thermal energy) trong dải nhiệt độ từ 60°C đến 100°C—như năng lượng bức xạ mặt trời hoặc nhiệt thải từ các lò hơi công nghiệp—nổi lên như một công nghệ triển vọng nhờ cấu tạo không chi tiết chuyển động cơ học, độ tin cậy cao và chi phí bảo trì thấp. Tuy nhiên, rào cản kỹ thuật lớn nhất khiến chu trình ejector độc lập chưa thể thương mại hóa rộng rãi là hệ số hiệu quả năng lượng nhiệt ($COP_n$) rất thấp, thường chỉ dao động từ 0,1 đến 0,3.
Khoảng trống nghiên cứu (research gap) cốt lõi nằm ở sự thiếu hụt các mô hình toán học giải tích kết hợp mô phỏng khí động lực học 3D hoàn chỉnh nhằm tối ưu hóa đa thông số hình học ejector khi tích hợp vào chu trình máy lạnh có máy nén hơi truyền thống (Vapor Compression Refrigeration System - VCRS) thông qua cơ chế ghép tầng (cascade). Các nghiên cứu trước đây (như Dorantes & Lallemand, Chen et al.) chủ yếu tập trung vào kết nối trực tiếp máy nén vào ejector, dẫn đến các nhược điểm nghiêm trọng như hiện tượng ejector hút lẫn dầu bôi trơn máy nén, hơi quá nhiệt làm tăng đột biến tải nhiệt bình ngưng và gây mất ổn định thủy lực. Luận án của nghiên cứu sinh Nguyễn Trung Kiên dưới sự hướng dẫn của GS. TS Lê Chí Hiệp tại Trường Đại học Bách Khoa – ĐHQG-HCM đã giải quyết triệt để vấn đề này bằng cách thiết lập cấu trúc ghép tầng gián tiếp qua bình trao đổi nhiệt trung gian (cascade intercooler).
Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác định chặt chẽ:
- Câu hỏi nghiên cứu 1 (RQ1): Mối quan hệ định lượng giữa kích thước hình học tối ưu của ejector siêu âm và các thông số vận hành động học trong chu trình tích hợp ghép tầng là gì?
- Giả thuyết 1 (H1): Tồn tại hàm tương quan hồi quy tuyến tính đa biến cho phép xác định chính xác đường kính cổ ống phun, đường kính cửa ra và đường kính buồng hòa trộn với độ tin cậy $R^2 > 80%$.
- Câu hỏi nghiên cứu 2 (RQ2): Vị trí tổn thất exergy cực đại tập trung ở thành phần nào và điều kiện vận hành nào cực tiểu hóa tính bất thuận nghịch của hệ thống?
- Giả thuyết 2 (H2): Ejector là nguồn phát sinh tổn thất exergy lớn nhất do ma sát nhớt, sóng xung kích và hòa trộn không đẳng áp; tồn tại một nhiệt độ trung gian tối ưu tại đó tổn thất exergy tổng thể đạt cực tiểu.
- Câu hỏi nghiên cứu 3 (RQ3): Mức độ cải thiện hiệu năng năng lượng ($COP_c$, $COP_n$) và khả năng tiết kiệm điện năng của chu trình tích hợp so với các chu trình đơn lẻ tương ứng trong điều kiện khí hậu nhiệt đới Việt Nam đạt được bao nhiêu?
- Giả thuyết 3 (H3): Chu trình ghép tầng kép cho phép giảm công nén trên 35% và tiết kiệm từ 20% đến 30% điện năng tiêu thụ thực tế.
Khung lý thuyết của nghiên cứu dựa trên sự giao thoa giữa Định luật I và II của Nhiệt động lực học (Thermodynamic First and Second Laws), Lý thuyết dòng khí nén được siêu âm (Compressible Gas Dynamics), và Lý thuyết phân tích Exergy (Gouy-Stodola Theorem). Đóng góp đột phá của luận án đã được định lượng hóa rõ ràng: "Trong điều kiện thí nghiệm, COPc tăng thêm 30,8% khi so sánh với COP của chu trình máy lạnh có máy nén hơi truyền thống và COPn tăng 122% so với COP của chu trình ejector truyền thống. Ngoài ra, công nén giảm 39,01% so với chu trình máy nén hơi thông thường." Phạm vi nghiên cứu áp dụng cho dải năng suất lạnh công nghiệp nhỏ và trung bình từ 2 kW đến 15 kW, với mô hình thực nghiệm kiểm chứng chế tạo ở công suất 2,8 kW, sử dụng cặp môi chất tối ưu R134a (tiểu chu trình ejector) và R410A (tiểu chu trình máy nén hơi).
Literature Review và Positioning
Các công trình nghiên cứu trên thế giới về công nghệ ejector nhiệt phân hóa thành hai trường phái chính. Trường phái thứ nhất tập trung vào việc hoàn thiện chu trình ejector đơn thuần bằng cách khảo sát động học dòng môi chất trong vòi phun siêu âm và buồng hòa trộn (Keenan et al., 1950; Eames et al., 1995; Huang et al., 1999). Các tác giả này đã xây dựng lý thuyết sóng áp suất 1D (constant-pressure and constant-area mixing models) nhằm dự đoán tỷ lệ lôi cuốn (entrainment ratio $\omega$). Tuy nhiên, hạn chế lớn của các mô hình này là giả định dòng một chiều lý tưởng, bỏ qua tổn thất xoáy cục bộ và sự dịch chuyển vị trí sóng xung kích (shock wave position) khi áp suất ngược thay đổi. Trường phái thứ hai khảo sát chu trình tích hợp ejector - máy nén hơi (Sokolov & Hershgal, 1990; Sun, 1997; Alexis, 2004) nhằm mục đích giải tỏa gánh nặng công nén cho máy nén cơ học.
Tranh luận học thuật cốt lõi diễn ra xung quanh phương thức ghép nối giữa chu trình ejector và chu trình máy nén hơi:
- Quan điểm ghép nối trực tiếp (Direct-coupling configuration): Các nghiên cứu của Arbel & Shitzer (2003) và Yan et al. (2015) đề xuất đưa trực tiếp dòng môi chất từ máy nén vào buồng hút của ejector. Tuy nhiên, phương án này vấp phải sự phản đối gay gắt từ các nhà thực nghiệm (như Pridasawas & Chaiyaporn, 2007) do nguy cơ dầu bôi trơn máy nén bị cuốn vào ejector gây tắc nghẽn vòi phun và làm sai lệch nghiêm trọng trường dòng chảy siêu âm, đồng thời khí xả quá nhiệt làm tăng đột biến phụ tải nhiệt tại bình ngưng tụ.
- Quan điểm ghép tầng gián tiếp (Cascade intercooler configuration): Phương án sử dụng bình trao đổi nhiệt trung gian để liên kết hai vòng tuần hoàn độc lập được đề xuất nhằm cô lập hoàn toàn dầu bôi trơn máy nén, cho phép lựa chọn độc lập hai môi chất lạnh chuyên biệt tối ưu cho từng cấp nhiệt độ.
Luận án định vị chính xác khoảng trống khoa học bằng cách lựa chọn cấu trúc ghép tầng gián tiếp, đồng thời giải quyết bài toán thiếu hụt công cụ thiết kế định lượng hình học ejector chuyên dụng cho dải nhiệt độ nguồn phát sinh thế thấp ($T_{ps} = 60^\circ\text{C} - 100^\circ\text{C}$). Để xác lập độ tin cậy khoa học, nghiên cứu đã tiến hành đối chuẩn (benchmarking) phương pháp mô phỏng số với nghiên cứu kinh điển của Scott et al. (2008) [25] sử dụng môi chất R141b. Kết quả cho thấy mô hình động lực học chất lưu tính toán 3D (3D CFD) do tác giả thiết lập tái hiện chính xác cấu trúc trường áp suất tĩnh, phân bố số Mach và hiện tượng nghẽn dòng cấp hai (secondary choking), vượt trội hoàn toàn so với các mô hình tính toán 1D truyền thống vốn đánh giá sai lệch tỷ lệ lôi cuốn từ 15% đến 30%. Luận án tiếp tục so sánh với các công bố của Chen et al. (2014) và Shen et al. (2017) để khẳng định tính ưu việt của cặp môi chất R134a-R410A trong điều kiện nhiệt độ môi trường nhiệt đới ($T_{nt} = 30^\circ\text{C} - 40^\circ\text{C}$).
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Nghiên cứu tạo ra bước tiến quan trọng trong việc mở rộng Lý thuyết Nhiệt động học Ghép tầng (Cascade Refrigeration Theory) và Lý thuyết Dòng khí nén Siêu âm (Compressible Aerodynamic Theory) thông qua các đóng góp cụ thể:
- Thiết lập mô hình toán học giải tích nhiệt động học toàn diện: Xây dựng hệ phương trình bảo toàn liên tục, động lượng và năng lượng khép kín cho chu trình ghép tầng, mô tả chính xác sự tương tác nhiệt - thủy lực giữa tiểu chu trình ejector (R134a) và tiểu chu trình máy nén hơi (R410A) qua bình trung gian.
- Lượng hóa bản chất tổn thất năng lượng thông qua phân tích Exergy: Xác lập quy luật phân bố tổn thất exergy cục bộ. Kết quả chỉ ra rằng tổn thất exergy trong ejector luôn chiếm tỷ trọng lớn nhất toàn hệ thống do ba cơ chế bất thuận nghịch nội tại: ma sát biên tại lớp biên của ống phun hội tụ - phân kỳ Laval, gradient vận tốc lớn trong quá trình hòa trộn dòng sơ cấp - thứ cấp, và hệ thống sóng xung kích (shock train) không đẳng entropy xảy ra tại buồng hòa trộn và ống khuếch tán. "Tổn thất exergy tổng nhỏ nhất đạt được là 1,5 kW tại nhiệt độ trung gian là 26,5°C trong điều kiện thí nghiệm."
- Chuyển đổi hình mẫu (Paradigm Shift) trong thiết kế hình học Ejector: Khắc phục hạn chế của phương pháp thiết kế thử - sai (trial-and-error) bằng cách phát triển bộ phương trình hồi quy tuyến tính đa biến có tính khái quát hóa cao, cho phép dự báo trực tiếp kích thước hình học tới hạn dựa trên các điều kiện biên nhiệt động học.
+-------------------------------------------------------------------------+
| KHUNG PHÂN TÍCH NHIỆT ĐỘNG VÀ KHÍ ĐỘNG |
+-------------------------------------------------------------------------+
|
+---------------------------+---------------------------+
| |
v v
+-------------------------------+ +-------------------------------+
| TIỂU CHU TRÌNH EJECTOR | | TIỂU CHU TRÌNH VCRS |
| (Môi chất nạp R134a) | | (Môi chất nạp R410A) |
| Nhiệt phát sinh: 60 - 100°C | | Nhiệt bay hơi: 0 - 10°C |
+-------------------------------+ +-------------------------------+
| |
| +-----------------------------------+ |
+-------->| BÌNH TRAO ĐỔI NHIỆT TRUNG GIAN |<--------+
| Nhiệt độ tối ưu: Ttg = 26,5°C |
| Tổn thất Exergy min: 1,5 kW |
+-----------------------------------+
|
v
+-------------------------------------------------------------------------+
| BỘ BA PHƯƠNG TRÌNH HỒI QUY TUYẾN TÍNH ĐA BIẾN |
| dt, d1r, dht = f(Tps, Tnt, Tbh, Ttg, Qbh) | R² = 84,98% - 89,71% |
+-------------------------------------------------------------------------+
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp chặt chẽ ba lý thuyết nền tảng:
- Lý thuyết dòng môi chất thực nén được (Real Gas Compressible Flow Theory): Áp dụng phương trình trạng thái nhiệt động Helmholtz thông qua liên kết thư viện số REFPROP với EES.
- Lý thuyết mô hình rối 2 phương trình (Launder & Spalding $k-\varepsilon$ Standard Model): Mô tả chính xác ứng suất Reynolds và hiện tượng khuếch tán rối trong lớp hòa trộn tự do siêu âm.
- Lý thuyết phân tích hiệu năng đa mục tiêu: Kết hợp đồng thời hai chỉ số hiệu quả năng lượng là hệ số COP điện năng ($COP_c = Q_{bh}/N_{mn}$) và hệ số COP nhiệt năng ($COP_n = Q_{bh}/(Q_{ps} + N_{mn} + N_b)$).
Phương trình hồi quy tuyến tính đa biến xác định kích thước hình học ejector được xác lập dưới dạng tổng quát:
$$D = f(T_{ps}, T_{nt}, T_{bh}, T_{tg}, Q_{bh}) = a_0 + a_1 T_{ps} + a_2 T_{nt} + a_3 T_{bh} + a_4 T_{tg} + a_5 Q_{bh}$$
Bộ hệ số hồi quy được tối ưu hóa với độ tin cậy thống kê cao:
- Đường kính cổ ống phun ($d_t$): $a_0 = 2,172316 \times 10^{-3}$, $a_1 = -6,051767 \times 10^{-5}$, $a_2 = 1,437332 \times 10^{-4}$, $a_3 = -7,965921 \times 10^{-6}$, $a_4 = -8,541482 \times 10^{-5}$, $a_5 = 5,542732 \times 10^{-4}$ với hệ số xác định $R^2 = 84,98%$.
- Đường kính cửa ra ống phun ($d_{1r}$): $a_0 = 3,609567 \times 10^{-3}$, $a_1 = -1,045933 \times 10^{-4}$, $a_2 = 2,882167 \times 10^{-4}$, $a_3 = -7,376264 \times 10^{-5}$, $a_4 = -1,886310 \times 10^{-4}$, $a_5 = 1,125799 \times 10^{-3}$ với hệ số xác định $R^2 = 87,00%$.
- Đường kính buồng hòa trộn ($d_{ht}$): $a_0 = 4,651666 \times 10^{-3}$, $a_1 = -9,910666 \times 10^{-5}$, $a_2 = 2,431008 \times 10^{-4}$, $a_3 = -1,151184 \times 10^{-4}$, $a_4 = -1,328218 \times 10^{-4}$, $a_5 = 1,385460 \times 10^{-4}$ với hệ số xác định $R^2 = 89,71%$.
Điều kiện biên áp dụng cho khung phân tích bao gồm dải nhiệt độ phát sinh $T_{ps} \in [60, 100]^\circ\text{C}$, nhiệt độ ngưng tụ $T_{nt} \in [30, 40]^\circ\text{C}$, nhiệt độ trung gian $T_{tg} \in [10, 20]^\circ\text{C}$, nhiệt độ bay hơi $T_{bh} \in [0, 10]^\circ\text{C}$, độ chênh nhiệt độ bình trung gian $\Delta T_{tg} = 5^\circ\text{C}$, và dải năng suất lạnh $Q_{bh} \in [0,5, 5]\text{ kW}$.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án áp dụng lập trường nhận thức luận thực chứng (Positivism) kết hợp với chủ nghĩa hiện thực phản biện (Critical Realism). Thiết kế nghiên cứu tuân thủ quy trình đa tầng tam giác hóa (Tri-level Triangulation Design) bao gồm ba cấu phần tương hỗ:
- Tầng 1 - Mô hình hóa lý thuyết nhiệt động 1D: Lập trình giải hệ phương trình bảo toàn phi tuyến trên phần mềm Engineering Equation Solver (EES) để xác định các trạng thái nhiệt động điểm nút và phân tích năng lượng - exergy.
- Tầng 2 - Mô phỏng số khí động lực học vi mô 3D CFD: Ứng dụng phần mềm thương mại ANSYS-FLUENT giải hệ phương trình Reynolds-Averaged Navier-Stokes (RANS) nhằm tối ưu hóa các thông số hình học vi mô không thể tính toán bằng mô hình 1D (tỷ lệ diện tích $TLDT$, chiều dài buồng hòa trộn $L_{ht}$, khoảng cách vòi phun $NXP$).
- Tầng 3 - Kiểm chứng thực nghiệm trên băng thử vật lý: Chế tạo hệ thống thử nghiệm thực tế công suất 2,8 kW với 4 cấu hình ejector chế tạo cơ khí chính xác để đo đạc và đối soát dữ liệu với mô phỏng.
Quy trình nghiên cứu rigorous
Quy trình mô phỏng số và thực nghiệm được thiết lập với tính chuẩn mực học thuật cao:
- Xây dựng hệ thống lưới và độc lập lưới (Grid Independence Study): Miền tính toán 2D đối xứng trục và 3D của ejector được rời rạc hóa bằng hệ thống lưới có cấu trúc (structured mesh) kết hợp làm mịn cục bộ (boundary layer mesh refinement) tại thành ống và vách vòi phun để đảm bảo giá trị $y^+ \approx 1$. Quá trình kiểm tra tính độc lập của lưới được thực hiện qua 4 cấp độ phân giải lưới từ 120.000 đến 450.000 phần tử.
- Thiết lập điều kiện biên và tiêu chuẩn hội tụ: Điều kiện biên đầu vào dòng sơ cấp và thứ cấp sử dụng áp suất tĩnh và nhiệt độ tổng (Pressure-inlet); ngõ ra sử dụng áp suất ngưng tụ (Pressure-outlet). Thuật toán giải COUPLED hoặc SIMPLEC được sử dụng cùng sơ đồ sai phân bậc hai (Second-Order Upwind). "Điều kiện hội tụ của ejector là các sai số giảm đến $10^{-6}$" đối với tất cả các phương trình liên tục, động lượng, năng lượng, $k$ và $\varepsilon$.
- Giao thức đo đạc thực nghiệm: Hệ thống trang bị cảm biến áp suất dải đo $0 - 40\text{ bar}$ (sai số $\pm 0,1%$), cảm biến nhiệt độ Pt100 cấp chính xác Class A (sai số $\pm 0,15^\circ\text{C}$), và đồng hồ đo lưu lượng khối lượng kiểu Coriolis (sai số $\pm 0,2%$). Nguồn nhiệt cấp cho bình phát sinh được giả lập ổn định thông qua hệ thống nước nóng gia nhiệt điều khiển PID chính xác $\pm 0,5^\circ\text{C}$.
Data và phân tích
- Đặc tính mẫu thực nghiệm: Chế tạo 4 mẫu ejector thử nghiệm (EJ1, EJ2, EJ3, EJ4) với tỷ lệ diện tích $TLDT = A_{ht}/A_t = (d_{ht}/d_t)^2$ biến thiên quanh giá trị thiết kế lý thuyết ($TLDT = 8,55$). Cụ thể, mẫu EJ3 sở hữu $TLDT = 8,46$, sát nhất với thông số tối ưu lý thuyết.
- Kiểm soát độ lệch và khoảng tin cậy: Phân tích độ không đảm bảo đo (uncertainty analysis) theo tiêu chuẩn ISO GUM cho thấy sai số tích lũy của phép đo COP không vượt quá $\pm 3,8%$.
- So sánh số học giữa thực nghiệm và mô phỏng: Sai số tương đối về tỷ lệ lôi cuốn ($\omega = m_2/m_1$) giữa mô phỏng ANSYS-FLUENT và thực nghiệm dao động trong dải dung sai cho phép từ $\pm 2,5%$ đến $\pm 11,5%$. Xu hướng dự báo của CFD bám sát thực nghiệm nhưng cho kết quả cao hơn nhẹ do mô phỏng chưa phản ánh toàn diện độ nhám bề mặt thực tế và tổn thất nhiệt cục bộ ra môi trường.
+-------------------------------------------------------------------------+
| QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU TAM GIÁC HÓA |
+-------------------------------------------------------------------------+
| |
| [PHÂN TÍCH NHIỆT ĐỘNG 1D - EES] |
| - Thiết lập cân bằng khối lượng & năng lượng |
| - Phân tích Exergy (Gouy-Stodola) |
| - Xuất bộ thông số hình học tham chiếu (TLDT = 8,55) |
| │ |
| ▼ |
| [MÔ PHỎNG KHÍ ĐỘNG 3D CFD - ANSYS FLUENT] |
| - Phương trình RANS + Mô hình rối k-ε |
| - Hội tụ sai số dư ≤ 10⁻⁶; y⁺ ≈ 1 |
| - Tối ưu hóa: TLDT = 8,88; Lht = (2,38 - 5,08)dht |
| │ |
| ▼ |
| [KIỂM CHỨNG BĂNG THỬ VẬT LÝ 2,8 kW] |
| - Chế tạo 4 cấu hình Ejector: EJ1, EJ2, EJ3 (TLDT = 8,46), EJ4 |
| - Kiểm chứng sai số: ±2,5% đến ±11,5% |
| - Lượng hóa tiết kiệm điện: 21,7% - 30,7% |
| |
+-------------------------------------------------------------------------+
Phát hiện đột phá và implications
Những phát hiện then chốt
Luận án mang lại 5 phát hiện khoa học có giá trị then chốt:
- Cơ chế nghẽn dòng và áp suất ngưng tụ tới hạn ($P_{nt}^*$): Thực nghiệm chứng minh khi áp suất ngưng tụ $P_{nt} < P_{nt}^$, tỷ lệ lôi cuốn $\omega$ duy trì không đổi (chế độ vận hành tới hạn - critical mode). Khi $P_{nt} > P_{nt}^$, sóng xung kích bị đẩy ngược vào buồng hòa trộn làm phá hủy chế độ nghẽn dòng cấp hai, dẫn đến $\omega$ suy giảm phi tuyến tính đột ngột về 0 và gây hiện tượng dòng chảy ngược (back-flow).
- Quy luật tác động của áp suất phát sinh ($P_{ps}$): Tỷ lệ lôi cuốn $\omega$ tăng tuyến tính theo áp suất phát sinh và đạt giá trị bão hòa tối ưu khi tiệm cận điểm tới hạn của R134a. Áp suất phát sinh tối ưu tương ứng cho các ejector được lượng hóa chính xác: EJ1 đạt 23 bar; EJ2 và EJ4 đạt 25 bar; EJ3 đạt 27 bar.
- Độ nhạy của khoảng cách ống phun - buồng hòa trộn ($NXP$): Hiệu quả lôi cuốn phụ thuộc cực kỳ nhạy vào vị trí tương đối giữa vòi phun và buồng trộn. Khoảng cách tối ưu xác định trên băng thử lần lượt là: 2 mm (EJ1), 3 mm (EJ2), 3 mm (EJ3) và 5 mm (EJ4). Nếu $NXP$ quá nhỏ, tiết diện dòng thứ cấp bị bóp hẹp cản trở quá trình hút; nếu $NXP$ quá lớn, độ đồng trục bị suy giảm và dòng sơ cấp mở rộng sớm làm tăng tổn thất va đập.
- Hình học tối ưu hóa bằng mô phỏng số: Kết quả CFD khẳng định tỷ lệ diện tích $TLDT = 8,88$ và chiều dài thân buồng hòa trộn $L_{ht} = (2,38 - 5,08) d_{ht}$ đảm bảo quá trình hòa trộn và hệ thống sóng xung kích hoàn tất trọn vẹn bên trong buồng trộn, triệt tiêu xung động áp suất tại cửa vào ống khuếch tán.
- Hiệu năng năng lượng vượt trội của chu trình ghép tầng: Chu trình tích hợp cải thiện vượt bậc cả hai chỉ số hiệu suất: $COP_c$ tăng 30,8%, $COP_n$ tăng 122%, công nén của máy nén hơi giảm 39,01%, và "Chu trình tích hợp có thể thiết kiệm được 21,7% - 30,7% điện năng tiêu thụ so với chu trình máy lạnh có máy nén hơi thông thường."
| Thông số đánh giá hiệu năng |
Chu trình máy nén hơi đơn (VCRS) |
Chu trình Ejector đơn |
Chu trình tích hợp ghép tầng (Luận án) |
Mức độ cải thiện / Đột phá |
| Môi chất công tác |
R410A |
R134a |
R134a (Ejector) - R410A (Máy nén) |
Tối ưu hóa chuyên biệt từng cấp nhiệt |
| Hệ số hiệu quả nhiệt ($COP_n$) |
Không áp dụng |
0,15 - 0,25 |
Tăng thêm 122% |
Tận dụng nhiệt thải 60 - 100°C |
| Hệ số hiệu quả điện ($COP_c$) |
Chuẩn cơ sở ($100%$) |
Không áp dụng |
Tăng thêm 30,8% |
Giảm áp lực tiêu thụ điện lưới |
| Công nén yêu cầu ($N_{mn}$) |
Chuẩn cơ sở ($100%$) |
Không áp dụng |
Giảm 39,01% |
Giảm dung lượng và kích thước máy nén |
| Tỷ lệ tiết kiệm điện năng |
$0%$ (Gốc so sánh) |
Không áp dụng |
21,7% - 30,7% |
Lượng hóa trên băng thử 2,8 kW |
| Tổn thất Exergy tối thiểu |
Cao tại van tiết lưu |
Rất cao tại Ejector |
1,5 kW (tại $T_{tg} = 26,5^\circ\text{C}$) |
Cực tiểu hóa tính bất thuận nghịch |
Implications đa chiều
- Ý nghĩa học thuật: Cung cấp bộ dữ liệu thực nghiệm chuẩn hóa và công thức hồi quy giải tích kinh nghiệm phục vụ việc phát triển các thuật toán tối ưu hóa hệ thống nhiệt đa nguồn trong tương lai.
- Ý nghĩa kỹ thuật công nghiệp: Mở ra giải pháp chế tạo thiết bị điều hòa không khí lai (hybrid HVAC) cho các nhà máy công nghiệp sẵn có nguồn nhiệt thải (lò hơi, xưởng sấy, nhà máy nhiệt điện) và các tòa nhà xanh tích hợp dàn thu nhiệt mặt trời.
- Ý nghĩa chính sách và môi trường: Đóng góp luận cứ kỹ thuật cho Chương trình quốc gia về sử dụng năng lượng tiết kiệm và hiệu quả (VNEEP), hỗ trợ đạt mục tiêu Net Zero Carbon vào năm 2050 thông qua cắt giảm từ 21,7% đến 30,7% nhu cầu điện năng làm mát.
Limitations và Future Research
Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:
- Môi chất R134a: Dù sở hữu tính chất nhiệt động phù hợp với dải nhiệt độ nguồn nhiệt thấp ($70^\circ\text{C} - 100^\circ\text{C}$) và an toàn, R134a vẫn có chỉ số làm nóng lên toàn cầu tương đối cao ($GWP = 1430$), đối mặt với lộ trình cắt giảm trong tương lai theo Bản sửa đổi Kigali.
- Điều kiện thí nghiệm nguồn nhiệt: Việc sử dụng nguồn nước nóng gia nhiệt bằng điện trở để chủ động ổn định thông số kiểm chứng đã bỏ qua tính chất dao động bất định thời gian thực của bức xạ mặt trời hoặc lưu lượng nhiệt thải công nghiệp.
- Hình học ejector cố định (Fixed-geometry ejector): Các mẫu ejector chế tạo chưa tích hợp cơ cấu kim phun điều chỉnh tiết diện cổ vòi phun (spindle / needle actuator) để tự động thích ứng khi phụ tải lạnh thay đổi lớn.
Chương trình nghiên cứu mở rộng trong giai đoạn tiếp theo tập trung vào 4 định hướng:
- Nghiên cứu ứng dụng các môi chất lạnh hydrofluoroolefin (HFO) và hydrocacbon (HC) thân thiện môi trường thế hệ mới có GWP thấp như R1234ze(E), R1233zd(E), R290 và R600a.
- Khảo sát ảnh hưởng của độ quá nhiệt dòng sơ cấp trước khi vào vòi phun ejector nhằm triệt tiêu hiện tượng ngưng tụ sương cục bộ tại cổ vòi phun siêu âm.
- Phát triển mô hình tích hợp ejector vào hệ thống máy lạnh hấp thụ (absorption refrigeration) đa cấu phần.
- Thiết kế ejector có hình học biến thiên (Variable Geometry Ejector - VGE) điều khiển thông minh theo tải thực tế.
Tác động và ảnh hưởng
Công trình nghiên cứu đã được bảo chứng chất lượng qua hệ thống các bài báo khoa học xuất bản trên tạp chí quốc tế thuộc danh mục ISI/SCIE (Journal of Thermal Analysis and Calorimetry, IF = 4.756, Springer) và các kỷ yếu hội nghị quốc tế uy tín (GTSD 2018, TSME 2016, ASEAN Smart Grid Congress 2017). Về mặt công nghiệp, kết quả của luận án cung cấp hướng dẫn thiết kế cơ khí chính xác có khả năng chuyển giao trực tiếp cho các doanh nghiệp sản xuất thiết bị lạnh và điều hòa không khí tại Việt Nam, cho phép nội địa hóa các giải pháp thu hồi nhiệt thải công nghiệp với chi phí đầu tư hợp lý.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận cơ sở dữ liệu thực nghiệm chi tiết, hệ phương trình hồi quy đa biến và phương pháp luận mô phỏng 3D CFD chuẩn mực để mở rộng các nghiên cứu chuyên sâu về máy nén nhiệt.
- Kỹ sư R&D ngành HVAC/R: Sử dụng trực tiếp quy trình tính toán kích thước hình học ($d_t, d_{1r}, d_{ht}, L_{ht}, NXP$) để chế tạo các module ejector phụ trợ tiết kiệm điện năng cho hệ thống máy lạnh thương mại.
- Chủ đầu tư công nghiệp và Quản lý năng lượng: Nhận diện tiềm năng cắt giảm chi phí vận hành từ 21,7% đến 30,7% thông qua việc tận dụng nguồn nước nóng hoặc khói thải lò hơi sẵn có.
- Cơ quan quản lý nhà nước: Sử dụng các số liệu tiết kiệm năng lượng đã kiểm chứng làm cơ sở xây dựng quy chuẩn kỹ thuật và định mức tiêu hao năng lượng cho ngành công nghiệp nhiệt lạnh.
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ì?
Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Nhiệt động lực học Ghép tầng thông qua việc giải quyết trọn vẹn bài toán ghép nối gián tiếp giữa ejector nhiệt thế thấp và máy nén hơi. Luận án đã giải mã cơ chế phân bố tổn thất exergy bên trong ejector siêu âm và xác lập bộ 3 phương trình hồi quy tuyến tính đa biến đầu tiên tại Việt Nam liên kết trực tiếp kích thước hình học tới hạn ($d_t, d_{1r}, d_{ht}$) với 5 thông số vận hành động học với hệ số xác định $R^2$ đạt từ 84,98% đến 89,71%.
2. Đổi mới phương pháp luận nghiên cứu so với các công bố quốc tế trước đây?
Luận án kết hợp đồng thời mô hình nhiệt động 1D trên nền tảng EES với mô phỏng vi mô 3D CFD trên ANSYS-FLUENT và đối chuẩn trực tiếp với mô hình chuẩn quốc tế của Scott et al. (2008) [25] trên môi chất R141b. Khác biệt cốt lõi là luận án đã chuyển tiếp thành công từ kết quả mô phỏng CFD sang chế tạo thực tế 4 mẫu ejector kim loại chính xác và vận hành kiểm chứng trên băng thử ghép tầng 2,8 kW, thu hẹp sai số giữa lý thuyết và thực nghiệm xuống dải $\pm 2,5% - \pm 11,5%$.
3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực nghiệm?
Phát hiện phản trực giác nhất là việc tăng áp suất phát sinh $P_{ps}$ không phải lúc nào cũng làm tăng tỷ lệ lôi cuốn $\omega$. Tỷ lệ lôi cuốn chỉ tăng đến một giới hạn cực đại gần điểm tới hạn của môi chất (23 - 27 bar đối với R134a) rồi bão hòa; đồng thời, nếu áp suất phát sinh quá cao sẽ kéo theo sự suy giảm của áp suất ngưng tụ tới hạn $P_{nt}^*$, làm cho hệ thống trở nên cực kỳ nhạy cảm và dễ mất ổn định trước sự biến động nhiệt độ của môi trường giải nhiệt.
4. Nghiên cứu có cung cấp giao thức tái lập thực nghiệm chi tiết không?
Có. Luận án công bố chi tiết bảng thông số trạng thái nhiệt động tại 10 điểm nút của chu trình thực nghiệm (Bảng 4.2 bao gồm Áp suất $P$, Nhiệt độ $T$, Enthalpy $i$, Entropy $s$, Thể tích riêng $v$), kích thước hình học chi tiết của 4 mẫu ejector thử nghiệm, sơ đồ nguyên lý thủy lực, danh mục cảm biến đo lường và các điều kiện biên mô phỏng CFD (áp suất, mô hình rối $k-\varepsilon$, chuẩn hội tụ $10^{-6}$).
5. Lộ trình nghiên cứu phát triển 10 năm được vạch ra như thế nào?
Lộ trình 10 năm tập trung vào 3 mũi nhọn: (1) Thay thế môi chất R134a bằng các môi chất lạnh xanh HFO/HC (R1234ze, R290) có $GWP < 1$; (2) Phát triển ejector thông minh có vòi phun biến thiên tiết diện tự động điều khiển bằng thuật toán thích nghi; (3) Tích hợp hoàn chỉnh hệ thống thu nhiệt mặt trời thực tế và mở rộng quy mô công suất lên dải thương mại 50 kW - 100 kW phục vụ điều hòa không khí tòa nhà xanh.
Kết luận
- Xây dựng hoàn chỉnh mô hình toán học giải tích 1D trên EES và phần mềm tính toán nhiệt động cho hệ thống máy lạnh ghép tầng ejector – máy lạnh có máy nén hơi, chứng minh sự vượt trội của phương án ghép nối gián tiếp qua bình trao đổi nhiệt trung gian.
- Xác lập quy luật phân bố tổn thất exergy, chứng minh ejector chiếm tổn thất lớn nhất và xác định được điểm cực tiểu tổn thất exergy tổng thể của hệ thống đạt 1,5 kW tại nhiệt độ trung gian $T_{tg} = 26,5^\circ\text{C}$.
- Thiết lập thành công bộ 3 phương trình hồi quy tuyến tính đa biến tính toán kích thước hình học ejector R134a ($d_t, d_{1r}, d_{ht}$) theo các thông số vận hành với độ tin cậy $R^2$ từ 84,98% đến 89,71%.
- Ứng dụng mô phỏng 3D CFD trên ANSYS-FLUENT xác định được kích thước hình học tối ưu vi mô: tỷ lệ diện tích $TLDT = 8,88$, chiều dài thân buồng hòa trộn $L_{ht} = (2,38 - 5,08) d_{ht}$, và khoảng cách vòi phun tối ưu $NXP = 3\text{ mm}$.
- Chế tạo thành công mô hình thực nghiệm 2,8 kW, chứng minh chu trình tích hợp nâng cao $COP_c$ thêm 30,8%, tăng $COP_n$ lên 122%, giảm công nén 39,01%, và tiết kiệm từ 21,7% đến 30,7% điện năng tiêu thụ so với chu trình máy nén hơi tiêu chuẩn.
- Mở ra hướng tiếp cận bền vững và khả thi cao trong việc tận dụng các nguồn nhiệt thế thấp (60°C - 100°C) từ năng lượng mặt trời và nhiệt thải công nghiệp cho ngành điều hòa không khí tại Việt Nam và khu vực nhiệt đới.
