Tổng quan về luận án

Sự biến đổi khí hậu toàn cầu và sự suy giảm tầng ozone đã thúc đẩy cộng đồng khoa học quốc tế tìm kiếm các môi chất lạnh tự nhiên có tiềm năng làm nóng toàn cầu thấp (GWP) và tiềm năng phá hủy tầng ozone bằng không (ODP) để thay thế hoàn toàn các chất làm lạnh truyền thống như HCFC và HFC. Trong bối cảnh Hiệp định Montreal và Bản sửa đổi Kigali đặt ra lộ trình loại bỏ dần các chất làm lạnh tổng hợp, $\text{CO}_2$ (R744) nổi lên như một giải pháp thay thế hoàn hảo với các chỉ số môi trường lý tưởng ($\text{GWP} = 1$, $\text{ODP} = 0$, không độc hại, không cháy). Tuy nhiên, chu trình nhiệt động của $\text{CO}_2$ là chu trình trên tới hạn (transcritical cycle) với áp suất làm việc rất cao (vượt quá áp suất tới hạn 73,8 bar), dẫn đến tổn thất tiết lưu lớn và làm giảm đáng kể hệ số hiệu quả năng lượng (COP) khi hoạt động trong điều kiện nhiệt độ môi trường nhiệt đới từ 28 đến $32^\circ\text{C}$.

Khoảng trống nghiên cứu (research gap) lớn nhất nằm ở việc phần lớn các công trình trước đây tập trung vào chu trình $\text{CO}_2$ trên tới hạn áp suất rất cao (>90 đến 130 bar) cho hệ sinh thái công nghiệp hoặc siêu thị lớn, hoặc chỉ khảo sát truyền nhiệt trên một ống đơn/kênh đơn riêng lẻ. Thiếu vắng các nghiên cứu thực nghiệm và mô phỏng số đồng bộ về đặc tính truyền nhiệt của thiết bị bay hơi kênh micro (microchannel evaporator) dạng chùm ống dẹt nhiều cổng (multi-port flat tube) tích hợp trong hệ thống điều hòa không khí cỡ nhỏ ($\le 9000\text{ BTU/h} \approx 2600\text{ W}$) vận hành ở dải áp suất đầu đẩy trung bình thấp (74 – 90 bar). Luận án của tác giả Nguyễn Trọng Hiếu với đề tài "Nghiên cứu các đặc tính truyền nhiệt của thiết bị bay hơi kênh micro trong máy điều hoà không khí cỡ nhỏ dùng môi chất lạnh $\text{CO}_2$" đã giải quyết trọn vẹn khoảng trống học thuật và kỹ thuật này.

Nghiên cứu được định hướng bởi hệ thống câu hỏi và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Các đặc tính phân bố nhiệt độ, tổn thất áp suất và hệ số tỏa nhiệt đối lưu hai pha biến thiên như thế nào dọc theo chiều dài kênh micro có đường kính thủy lực $D_h = 0,9\text{ mm}$ khi môi chất $\text{CO}_2$ chuyển pha từ vùng sôi đến vùng quá nhiệt?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ suy giảm hiệu suất do tổn thất tiết lưu trong chu trình $\text{CO}_2$ trên tới hạn có thể được bù đắp và tối ưu hóa đến mức nào thông qua giải pháp tích hợp thiết bị làm mát phụ (Subcooler) và thiết bị hồi nhiệt (Internal Heat Exchanger - IHX)?
  • Câu hỏi nghiên cứu 3 (RQ3): Quy luật tương tác giữa lưu lượng khối môi chất ($\dot{m}_r$), vận tốc gió qua dàn bay hơi ($v_a$) và tỉ số áp suất ($p_c/p_e$) ảnh hưởng ra sao đến sự cân bằng truyền nhiệt giữa phía môi chất và phía không khí?

Các giả thuyết nghiên cứu tương ứng:

  • Giả thuyết 1 (H1): Việc cải tiến công thức tính hệ số ma sát $f$ dựa trên tỉ lệ kích thước mặt cắt kênh micro sẽ nâng cao độ chính xác dự báo hệ số tỏa nhiệt đối lưu hai pha và tổn thất áp suất so với các mô hình đồng nhất truyền thống.
  • Giả thuyết 2 (H2): Việc giảm enthalpy trước van tiết lưu bằng thiết bị hồi nhiệt IHX sẽ kiểm soát được độ khô đầu vào, mở rộng dải năng suất lạnh riêng và nâng cao COP vượt trội trên 100% so với chu trình cơ bản.
  • Giả thuyết 3 (H3): Tồn tại một ngưỡng lưu lượng khối và vận tốc không khí tối ưu giúp hệ số truyền nhiệt tổng $k$ và năng suất lạnh đạt cực trị trong điều kiện nhiệt độ môi trường thực tế tại Việt Nam.

Khung lý thuyết của luận án tích hợp sâu sắc giữa nhiệt động lực học cổ điển, cơ học chất lưu hai pha (Two-phase flow dynamics), mô hình truyền nhiệt sôi đối lưu của Choi et al. và hệ phương trình Navier-Stokes trung bình Reynolds (RANS $k-\epsilon$). Nghiên cứu mang lại đóng góp đột phá với những con số định lượng cụ thể: việc ứng dụng IHX giúp năng suất lạnh tăng 100% và COP tăng 103%; tích hợp Subcooler giúp giảm nhiệt độ trước van tiết lưu $1,4^\circ\text{C}$, gia tăng 50% năng suất lạnh và 39% COP. Hệ thống đạt công suất lạnh cực đại $3,12\text{ kW}$ với $\text{COP} = 3,15$ tại lưu lượng $\text{CO}_2$ tối ưu $111\text{ kg/h}$. Phạm vi nghiên cứu bao quát điều kiện vận hành với áp suất bay hơi 40 – 50 bar ($T_e = 5 - 15^\circ\text{C}$), áp suất đầu đẩy 74 – 90 bar, phục vụ trực tiếp cho các hệ thống điều hòa không khí dân dụng và cabin ô tô thế hệ mới.

Literature Review và Positioning

Nghiên cứu về chu trình làm lạnh $\text{CO}_2$ trên tới hạn và thiết bị trao đổi nhiệt vi kênh phát triển qua hai nhánh học thuật chính:

Nhánh thứ nhất: Nhiệt động lực học chu trình $\text{CO}_2$ trên tới hạn và các giải pháp cải thiện hiệu suất. Kim et al. (2004) đã chỉ ra rằng chu trình $\text{CO}_2$ trên tới hạn có hiệu suất thấp hơn chu trình R134a truyền thống do tổn thất giải nhiệt ở bộ làm mát khí (gas cooler) và tổn thất đẳng enthalpy tại van tiết lưu. Để khắc phục, Sanchez et al. (2014) thực nghiệm hệ thống có thiết bị làm mát phụ nhiệt điện, ghi nhận COP tăng 9,9% và năng suất lạnh tăng 16%. Llopis et al. (2015) và Andres et al. (2016) phát triển hệ thống làm mát cơ học chuyên dụng (Dedicated Mechanical Subcooling - DMS) sử dụng R152a, ghi nhận mức tăng COP từ 1,51 lên 1,95 ở $-1,3^\circ\text{C}$ và từ 1,86 lên 2,52 ở $10^\circ\text{C}$. Về thiết bị hồi nhiệt, Boewe et al. (2001) khẳng định IHX dòng chảy ngược chiều có thể tăng năng suất lạnh và COP lên tới 25% đối với các hệ thống $\text{CO}_2$ cỡ nhỏ. Elbel và Hrnjak (2008) ứng dụng Ejector để thu hồi động năng dòng tiết lưu, giúp tăng COP thêm 7%, trong khi Lee et al. (2014) ghi nhận mức tăng 15%.

Nhánh thứ hai: Cơ chế truyền nhiệt sôi và tổn thất áp suất của $\text{CO}_2$ trong kênh mini/micro. Kandlikar và Grande (2003) phân loại kênh micro có đường kính thủy lực $10 \le D \le 200,\mu\text{m}$, trong khi Mehendale et al. (2000) và các tác giả khác quy ước $D_h \le 1000,\mu\text{m}$. Choi et al. (2007) phát triển mô hình truyền nhiệt sôi hai pha cho $\text{CO}_2$ trong ống nhỏ, tích hợp thông số Martinelli và hệ số sôi boiling number ($Bo$). Yun et al. (2003, 2005) thực nghiệm trên kênh micro hình chữ nhật ($D_h = 1,08 - 1,54\text{ mm}$), chứng minh hệ số tỏa nhiệt đối lưu của $\text{CO}_2$ cao hơn 53% so với R134a. Oh và Son (2011), Ducoulombier et al. (2011), Pamitran et al. (2008) cung cấp các bộ dữ liệu thực nghiệm phong phú xác nhận sự suy giảm mạnh của hệ số truyền nhiệt khi độ khô $x$ tiến gần tới trạng thái bốc hơi hoàn toàn (dryout transition).

                      ┌────────────────────────────────────────┐
                      │    Chu trình CO2 trên tới hạn cơ bản   │
                      │  (Hiệu suất thấp do tổn thất tiết lưu) │
                      └───────────────────┬────────────────────┘
                                          │
                  ┌───────────────────────┴───────────────────────┐
                  ▼                                               ▼
   ┌──────────────────────────────┐                ┌──────────────────────────────┐
   │ Các giải pháp chu trình cũ   │                │   Tranh luận học thuật       │
   │ - Ejector (Elbel & Hrnjak)   │                │ - Flash gas bypass (Kim/Beaver)
   │ - Subcooler cơ học (Llopis)  │                │ - Mô hình ma sát Blasius     │
   │ - IHX đơn lẻ (Boewe et al.)  │                │   cho vi kênh (Choi et al.)  │
   └──────────────┬───────────────┘                └──────────────┬───────────────┘
                  │                                               │
                  └───────────────────────┬───────────────────────┘
                                          │
                                          ▼
                      ┌────────────────────────────────────────┐
                      │      POSITIONING CỦA LUẬN ÁN NÀY       │
                      │ • TBBH 10 kênh micro nhôm Dh = 0,9 mm  │
                      │ • Áp suất thấp trên tới hạn 74–90 bar │
                      │ • Khảo sát song song IHX & Subcooler   │
                      │ • Hiệu chỉnh f theo tỉ lệ cạnh ac      │
                      │ • Đồng bộ: Lý thuyết - CFD - Thực nghiệm
                      └────────────────────────────────────────┘

Trong y văn tồn tại hai cuộc tranh luận học thuật lớn:

  1. Tranh luận về cấu hình Flash Gas Bypass và phân bố dòng chảy: Kim et al. (2004) lập luận rằng thiết bị tách pha và bypass pha khí không phù hợp cho dàn bay hơi kênh micro vì gây mất cân bằng phân phối dòng hai pha giữa các ống góp (headers). Ngược lại, Beaver et al. (1999) cho rằng van bypass là cần thiết để kiểm soát độ quá nhiệt. Luận án này định vị việc loại bỏ bypass phức tạp và thay bằng kiểm soát độ khô đầu vào thông qua thiết bị hồi nhiệt IHX hoặc làm mát phụ, giải quyết triệt để vấn đề phân bố không đồng đều.
  2. Tranh luận về mô hình hệ số ma sát trong dòng sôi hai pha: Choi et al. (2007) áp dụng hệ số ma sát Blasius ($f = 0,079 Re^{-0,25}$) đồng nhất cho cả hai pha lỏng và hơi. Tuy nhiên, nhiều học giả chỉ ra rằng trong ống dẹt tiết diện chữ nhật nhỏ, hiệu ứng hình học và độ nhám bề mặt làm cho công thức Blasius đánh giá sai lệch tổn thất áp suất và tham số Martinelli $X$.

So sánh với các công trình quốc tế điển hình:

  • Nghiên cứu của Zhao et al. (2000) trên ống tròn đơn $D_h = 0,86\text{ mm}$ chỉ ra tổn thất áp suất $\text{CO}_2$ thấp hơn 40% so với R134a nhưng chưa khảo sát kết cấu dàn hoàn chỉnh có cánh tản nhiệt louver.
  • Công trình của Pamitran et al. (2008) khảo sát ống đơn bằng thép không gỉ $D = 1,5\text{ mm}$ ở dải lưu lượng $200 - 600\text{ kg/m}^2\text{s}$, nhưng không phân tích ảnh hưởng của vận tốc không khí và độ ẩm phía ngoài ống dẹt.

Luận án này đã tiến xa hơn bằng cách khảo sát toàn diện một thiết bị bay hơi nhôm hoàn chỉnh gồm 10 kênh micro song song ($D_h = 0,9\text{ mm}$), có cánh tản nhiệt chữ nhật, đặt trong hệ thống chu trình khép kín hoạt động ở áp suất thấp trên tới hạn (74 – 90 bar), kết hợp phân tích đồng thời cả phía môi chất và phía không khí ẩm.

Đó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 trực tiếp mô hình truyền nhiệt sôi đối lưu của Choi et al. (2007) bằng cách tái cấu trúc cách xác định hệ số ma sát hai pha $f$. Thay vì sử dụng công thức Blasius đồng nhất cho dòng môi chất, nghiên cứu đề xuất phân tách rõ ràng chế độ chảy của từng pha và đưa vào tỉ số hình học của kênh micro chữ nhật $a_c$ (tỉ lệ cạnh ngắn trên cạnh dài):

Khi $Re < 4000$ (chế độ dòng chảy tầng và quá độ), số Poiseuille được hiệu chỉnh chính xác theo hình học mặt cắt: $$Po = f \cdot Re = 24\left(1 - 1,3553a_c + 1,9467a_c^2 - 1,7012a_c^3 + 0,9564a_c^4 - 0,2537a_c^5\right)$$

Khi $Re \ge 4000$ (chế độ dòng chảy rối), hệ số ma sát được tính theo hàm độ nhám tương đối $\varepsilon/D_h$: $$f = \frac{0,25}{\left[\ln\left(\frac{\varepsilon/D_h}{3,7} + \frac{5,74}{Re^{0,9}}\right)\right]^2}$$

Sự cải tiến này làm thay đổi trực tiếp giá trị của tham số Martinelli: $$X = \left(\frac{dp/dz|_l}{dp/dz|_g}\right)^{0,5} = \left(\frac{f_l}{f_g}\right)^{0,5} \left(\frac{1-x}{x}\right) \left(\frac{\rho_g}{\rho_l}\right)^{0,5}$$

Từ đó, hệ số tăng cường đối lưu $F$ và hệ số dập tắt sôi $S$ trong phương trình hệ số tỏa nhiệt đối lưu hai pha $\alpha_{r,tp} = S \cdot \alpha_{NB} + F \cdot \alpha_{lo}$ được xác định chính xác hơn, phản ánh đúng cơ chế chuyển đổi từ sôi bọt (nucleate boiling) sang bốc hơi màng cưỡng bức (forced convective evaporation).

┌─────────────────────────────────────────────────────────────────────────────┐
│                    KHUNG PHÂN TÍCH NHIỆT - THỦY LỰC ĐỒNG BỘ                │
├──────────────────────────────────────┬──────────────────────────────────────┤
│       PHÂN VÙNG 1: SÔI HAI PHA       │     PHÂN VÙNG 2: HƠI QUÁ NHIỆT       │
├──────────────────────────────────────┼──────────────────────────────────────┤
│ • Nhiệt ẩn hóa hơi $h_{lv}$          │ • Nhiệt dung riêng $c_p$ pha hơi     │
│ • Hệ số tỏa nhiệt $\alpha_{r,tp}$    │ • Hệ số tỏa nhiệt $\alpha_{r,sp}$    │
│ • Tham số Martinelli $X(a_c, x)$     │ • Số Reynolds pha khí $Re_g$         │
│ • Tổn thất áp suất $\Delta p_{tp}$   │ • Tổn thất áp suất 1 pha $\Delta p_1$│
│   (ma sát + gia tốc dòng chảy)       │   (ma sát đơn thuần)                 │
├──────────────────────────────────────┴──────────────────────────────────────┤
│                              LIÊN KẾT BIÊN GIAO DIỆN                        │
│                   Độ khô $x = 1,0$ (Vị trí bốc hơi hoàn toàn)               │
├─────────────────────────────────────────────────────────────────────────────┤
│                          PHÍA NGOÀI: KHÔNG KHÍ ẨM                           │
│ • Hệ số tỏa nhiệt không khí khô $\alpha_{a,dry} = f(Re_a, S_p, D_o, h_f)$   │
│ • Hiệu suất cánh tản nhiệt $\eta_f = \tanh(ml)/(ml)$                        │
│ • Hệ số tách ẩm ngưng tụ $\xi = 1 + 2480 (d_{in} - d_{out}) / Q_s$          │
│ • Hệ số truyền nhiệt tổng $k = [1/\alpha_a + R_b + 1/(\beta_e \alpha_r)]^{-1}$│
└─────────────────────────────────────────────────────────────────────────────┘

Hệ thống giả thuyết và mệnh đề lý thuyết:

  • Mệnh đề 1 (P1): Sự phân chia ranh giới giữa vùng truyền nhiệt hai pha ($A_{e1}$) và vùng quá nhiệt ($A_{e2}$) trên diện tích trao đổi nhiệt $A_e = A_{e1} + A_{e2}$ quyết định trực tiếp hiệu suất tổng thể, trong đó hệ số truyền nhiệt $k_{e1}$ ở vùng hai pha cao gấp 3 đến 5 lần so với $k_{e2}$ ở vùng quá nhiệt.
  • Mệnh đề 2 (P2): Tổn thất áp suất tổng thể $\Delta p$ trong kênh micro là tổng của 6 thành phần riêng biệt: tổn thất co hẹp đầu vào ($\Delta p_c$), ma sát một pha ($\Delta p_{fr,1ph}$), ma sát hai pha ($\Delta p_{fr,tp}$), gia tốc dòng chảy do bốc hơi ($\Delta p_{ac}$), trọng trường ($\Delta p_{gr}$) và mở rộng đầu ra ($\Delta p_e$).

Nghiên cứu tạo ra sự dịch chuyển nhận thức (paradigm shift) khi chứng minh rằng chu trình $\text{CO}_2$ trên tới hạn cỡ nhỏ hoàn toàn có thể đạt hiệu quả năng lượng cao ($\text{COP} > 3,1$) ở khí hậu nhiệt đới mà không cần nâng áp suất lên mức nguy hiểm (>100 bar), miễn là áp dụng đúng kỹ thuật hạ nhiệt độ trước van tiết lưu.

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

Khung phân tích của luận án tích hợp đồng thời ba lý thuyết nền tảng:

  1. Định luật nhiệt động I và II áp dụng cho chu trình trên tới hạn có thiết bị phụ trợ (Subcooler/IHX).
  2. Lý thuyết dòng chảy hai pha Martinelli-Lockhart và mô hình sôi vi kênh.
  3. Lý thuyết truyền nhiệt - truyền ẩm đồng thời trên bề mặt cánh tản nhiệt có xét đến hệ số đọng sương $\xi$ và hiệu suất cánh $\eta_f$.

Điều kiện biên xác lập chặt chẽ: Áp suất làm mát khí $P_c \in [74, 90]\text{ bar}$, áp suất bay hơi $P_e \in [40, 50]\text{ bar}$ (tương ứng $T_e = 5 - 15^\circ\text{C}$), lưu lượng môi chất $\dot{m}_r \in [97,5; 121,4]\text{ kg/h}$, và dải vận tốc gió mặt dàn $v_a \in [0,5; 5,1]\text{ m/s}$.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (positivism) và chủ nghĩa hiện thực phê phán (critical realism), xây dựng quy trình nghiên cứu tam giác hóa (triangulation) khép kín gồm: Tính toán giải tích lý thuyết (Analytical Modeling) $\rightarrow$ Mô phỏng số động lực học chất lưu đa vật lý (CFD Simulation) $\rightarrow$ Thực nghiệm kiểm chứng trên mô hình thực tế (Experimental Verification).

   ┌─────────────────────────────────────────────────────────────────────────┐
   │                       THIẾT KẾ NGHIÊN CỨU TAM GIÁC HÓA                  │
   └────────────────────────────────────┬────────────────────────────────────┘
                                        │
         ┌──────────────────────────────┼──────────────────────────────┐
         ▼                              ▼                              ▼
┌──────────────────┐          ┌──────────────────┐          ┌──────────────────┐
│ PHƯƠNG PHÁP SỐ   │          │ TÍNH TOÁN LÝ THUYẾT│         │ THỰC NGHIỆM      │
│ • COMSOL V6.0    │◄────────►│ • Phương trình EES│◄────────►│ • Bench test     │
│ • RANS k-ε       │          │ • Excel Solver   │          │ • Dorin CD 180H  │
│ • Grid 5 cấp độ  │          │ • Martinelli cải │          │ • Danfoss TBBH   │
│   (14.260 elements)         │   tiến (ac)      │          │ • Coriolis meter │
└────────┬─────────┘          └────────┬─────────┘          └────────┬─────────┘
         │                             │                             │
         └─────────────────────────────┼─────────────────────────────┘
                                       ▼
                     ┌──────────────────────────────────┐
                     │ ĐÁNH GIÁ ĐỘ LỆCH VÀ SAI SỐ       │
                     │ • Sai số độ khô max: 3,5% - 5,5% │
                     │ • Δp: 1,13 vs 1,40 vs 1,50 bar   │
                     │ • α_r,tp: 6,5 -> 1,3 kW/m²K      │
                     └──────────────────────────────────┘

Thiết kế đa cấp độ (Multi-level design) được cấu trúc:

  • Cấp vi mô (Micro-scale): Dòng chảy hai pha và chuyển pha trong 1 kênh micro đơn ($D_h = 0,9\text{ mm}$, tiết diện chữ nhật).
  • Cấp cụm thiết bị (Component-scale): Thiết bị bay hơi ống dẹt 10 kênh micro bằng nhôm, tích hợp cánh tản nhiệt louver với bước cánh $S_p = 1,0 - 1,2\text{ mm}$, chiều dày cánh $\delta_f = 0,1\text{ mm}$.
  • Cấp hệ thống (System-scale): Toàn bộ chu trình $\text{CO}_2$ trên tới hạn công suất danh định 2,6 kW (9000 BTU/h) lắp đặt hoàn chỉnh.

Quy trình nghiên cứu rigorous

Hệ thống thực nghiệm được chế tạo và lắp đặt tại Phòng thí nghiệm Truyền nhiệt (Heat Transfer Lab) - Trường Đại học Sư phạm Kỹ thuật TP.HCM với các thiết bị tiêu chuẩn công nghiệp cao cấp:

  • Máy nén: Máy nén bán kín chuyên dụng cho $\text{CO}_2$ trên tới hạn Dorin CD 180H.
  • Dàn bay hơi: Dàn kênh micro nhôm Danfoss 10 kênh, chịu áp suất phá hủy trên 100 bar.
  • Bộ làm mát khí: Thiết bị làm mát Panasonic chuyên dụng cho môi chất R744.
  • Dụng cụ đo lường: Cảm biến đo lưu lượng khối lượng dạng Coriolis (độ chính xác $\pm 0,1%$), cảm biến áp suất dải đo 0 – 160 bar (độ chính xác $\pm 0,25%$), cảm biến nhiệt độ RTD Pt100 chuẩn loại A (sai số $\pm 0,15^\circ\text{C}$), máy đo vận tốc gió cánh quạt Testo và ẩm kế điện tử.

Quy trình thu thập dữ liệu tuân thủ trạng thái xác lập (steady-state condition): mỗi chế độ vận hành được giữ ổn định tối thiểu 30 phút trước khi ghi nhận 100 điểm dữ liệu liên tục để lấy giá trị trung bình thống kê. Độ tin cậy và tính hợp lệ của phép đo được kiểm soát qua phân tích lan truyền sai số (error propagation analysis) theo tiêu chuẩn ASME PTC 19.1.

Data và phân tích

Mô phỏng số đa vật lý: Sử dụng phần mềm COMSOL Multiphysics v6.0, giải đồng thời phương trình liên tục, phương trình động lượng RANS $k-\epsilon$, và phương trình năng lượng kết hợp mô hình sôi hai pha. Lưới tính toán được chia theo phương pháp phần tử hữu hạn (FEM) với 5 cấp độ kiểm tra tính độc lập của lưới, tinh chỉnh lớp biên sát vách ống ($y^+ \approx 1$). Số bậc tự do giải toán vượt trên $200.000$.

Tính toán lý thuyết nhiệt động: Xây dựng chương trình tính toán tối ưu hóa trên phần mềm EES (Engineering Equation Solver) và bảng tính Microsoft Excel tự động hóa, liên kết trực tiếp với cơ sở dữ liệu thuộc tính nhiệt vật lý REFPROP của NIST đối với môi chất R744.

Độ tin cậy được chứng minh bằng sự trùng khớp ấn tượng giữa ba phương pháp:

  • Sai số độ khô cực đại ($x$) giữa mô phỏng số CFD và thực nghiệm là 5,5%.
  • Sai số độ khô cực đại giữa tính toán lý thuyết và thực nghiệm chỉ 3,5%.
  • Tổn thất áp suất ghi nhận lần lượt là 1,13 bar (tính toán), 1,40 bar (mô phỏng CFD) và 1,50 bar (thực nghiệm).
  • Hệ số tỏa nhiệt đối lưu hai pha $\alpha_{r,tp}$ tính toán và mô phỏng hoàn toàn nằm trong dải dữ liệu thực nghiệm từ 6,5 xuống 1,3 kW/m²K với dung sai cực hẹp $\pm 1,5\text{ kW/m}^2\text{K}$.

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

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

      BIẾN THIÊN HỆ SỐ TỎA NHIỆT ĐỐI LƯỚI HAI PHA α_r,tp DỌC THEO ĐỘ KHÔ x
  α (kW/m²K)
   8 ┌─────────────────────────────────────────────────────────────┐
     │  ▲                                                          │
   6 │ ─┼─ (x ≈ 0,15 - 0,4: Vùng sôi bọt cực đại ~6,5 kW/m²K)       │
     │   \                                                         │
   4 │    \                                                        │
     │     \                                                       │
   2 │      \___ (x > 0,85: Dryout transition & bốc hơi hoàn toàn) │
     │          \                                                  │
   0 └───────────┴─────────────────────────────────────────────────┴──►
     0          0,2         0,4         0,6         0,8         1,0   Độ khô x
  1. Quy luật suy giảm phi tuyến của hệ số tỏa nhiệt đối lưu $\alpha_{r,tp}$: Hệ số $\alpha_{r,tp}$ đạt cực đại xấp xỉ $6,5\text{ kW/m}^2\text{K}$ ở giai đoạn đầu kênh (khi độ khô $x \approx 0,15 - 0,4$, nơi cơ chế sôi bọt chiếm ưu thế), sau đó giảm dần đều xuống $1,3\text{ kW/m}^2\text{K}$ ở cuối kênh khi $x \rightarrow 1,0$ và môi chất chuyển sang trạng thái hơi quá nhiệt. Hiện tượng này khẳng định sự suy giảm màng chất lỏng sát vách kênh (film dryout) ở vận tốc cao của môi chất $\text{CO}_2$.
  2. Hiệu ứng đột phá của thiết bị hồi nhiệt (IHX): Việc lắp đặt IHX tạo ra bước nhảy vọt về hiệu năng: năng suất lạnh tăng chính xác 100% (từ 1,05 kW lên 2,10 kW ở chế độ đối chứng) và hệ số COP tăng 103% (từ 1,28 lên 2,60). Nguyên nhân là do dòng gas nóng áp suất cao được làm lạnh sâu trước khi vào van tiết lưu, làm giảm độ khô đầu vào $x_{in}$ từ 0,45 xuống dưới 0,25, giúp gia tăng tối đa nhiệt ẩn hóa hơi hữu ích trong dàn bay hơi.
  3. Tác động của thiết bị làm mát phụ (Subcooler): Lắp đặt Subcooler làm giảm nhiệt độ gas trước van tiết lưu $1,4^\circ\text{C}$, giúp năng suất lạnh tăng 50% và COP tăng 39%, chứng minh tính khả thi của việc giải nhiệt phụ bằng nguồn nước ngầm hoặc không khí đối lưu cưỡng bức.
  4. Điểm tối ưu vận hành của lưu lượng khối $\text{CO}_2$: Khi điều chỉnh lưu lượng môi chất từ $97,5$ đến $121,4\text{ kg/h}$, nhiệt độ bay hơi $T_e$ tăng tương ứng từ $8,2^\circ\text{C}$ lên $14,5^\circ\text{C}$. Hệ thống xác lập điểm cực trị hiệu quả tại $\dot{m}_r = 111\text{ kg/h}$, đạt năng suất lạnh tối đa 3,12 kW và $\text{COP} = 3,15$.
  5. Điểm cân bằng nhiệt động giữa môi chất và không khí ẩm: Khi thay đổi vận tốc gió qua dàn bay hơi từ $0,5$ đến $5,1\text{ m/s}$ tại lưu lượng $\text{CO}_2$ cố định $75,6\text{ kg/h}$, tại giá trị vận tốc $v_a = 5,1\text{ m/s}$, công suất thu nhiệt phía không khí đạt trạng thái cân bằng tuyệt đối với công suất trao đổi nhiệt phía môi chất ở mức 2,09 kW.
Thông số đánh giá Chu trình cơ bản Tích hợp Subcooler Tích hợp IHX Cải thiện lớn nhất
Nhiệt độ trước van tiết lưu ($T_3$) $32,0^\circ\text{C}$ $30,6^\circ\text{C}$ (giảm $1,4^\circ\text{C}$) $22,4^\circ\text{C}$ Giảm $9,6^\circ\text{C}$
Độ khô đầu vào dàn bay hơi ($x_{in}$) 0,48 0,36 0,22 Giảm 54,1%
Năng suất lạnh ($Q_e$) 1,05 kW 1,58 kW (+50%) 2,10 kW (+100%) Đạt 3,12 kW (tại 111 kg/h)
Hệ số hiệu quả năng lượng (COP) 1,28 1,78 (+39%) 2,60 (+103%) Đạt 3,15 (tại 111 kg/h)
Tổn thất áp suất dàn bay hơi ($\Delta p$) 1,50 bar 1,48 bar 1,42 bar Duy trì ổn định ~1,4 bar

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Cung cấp phương pháp luận chuẩn xác trong việc hiệu chỉnh hệ số ma sát vi kênh theo tỉ số cạnh $a_c$, hoàn thiện lý thuyết truyền nhiệt hai pha cho chất lưu có áp suất thu gọn cao ($P/P_{crit} > 0,5$).
  • Về mặt phương pháp nghiên cứu: Thiết lập quy trình chuẩn kết hợp CFD RANS $k-\epsilon$ trên COMSOL với mã giải tích EES để giải bài toán trao đổi nhiệt vi kênh mà không cần chi phí thử nghiệm tốn kém cho mọi biến thể hình học.
  • Về ứng dụng thực tiễn và kỹ thuật: Cung cấp thông số thiết kế chuẩn xác cho các kỹ sư chế tạo dàn lạnh điều hòa ô tô điện (EVs), nơi môi chất $\text{CO}_2$ đang là xu thế bắt buộc để thay thế HFO-1234yf đắt tiền và độc hại khi cháy.
  • Về chính sách môi trường: Cung cấp bằng chứng thực nghiệm vững chắc để các nhà hoạch định chính sách tại Việt Nam và Đông Nam Á xây dựng tiêu chuẩn năng lượng cho các thiết bị lạnh sử dụng gas tự nhiên $\text{CO}_2$.

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 kỹ thuật:

  1. Giới hạn hình học kênh: Nghiên cứu thực nghiệm cố định trên kết cấu kênh dẹt nhôm có $D_h = 0,9\text{ mm}$, chưa mở rộng thực nghiệm đa dạng các tiết diện tròn, elip hoặc tam giác trên cùng một hệ thống thực tế.
  2. Ảnh hưởng của dầu bôi trơn máy nén: Mô phỏng số và tính toán lý thuyết đang giả định môi chất $\text{CO}_2$ tinh khiết hoàn toàn, bỏ qua ảnh hưởng của màng dầu bôi trơn PAG/POE lưu chuyển trong hệ thống đến nhiệt trở $R_b$.
  3. Dải áp suất kiểm tra: Thí nghiệm giới hạn ở dải áp suất trên tới hạn trung bình (74 – 90 bar), chưa khảo sát phản ứng của thiết bị ở dải áp suất cực hạn (100 – 120 bar) do giới hạn an toàn phòng thí nghiệm.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  • Phát triển thuật toán điều khiển thông minh tự động điều chỉnh độ mở van tiết lưu điện tử (EEV) kết hợp biến tần máy nén để bám sát điểm áp suất tối ưu $P_{opt}$ theo thời gian thực.
  • Khảo sát sự hòa trộn dầu - môi chất $\text{CO}_2$ và tác động của nồng độ dầu (Oil Circulation Ratio - OCR) đến hệ số tỏa nhiệt đối lưu trong kênh micro.
  • Ứng dụng công nghệ gia công vi cơ điện tử (MEMS) hoặc in kim loại 3D để tạo ra các kênh micro có cánh tản nhiệt vi mô dạng lượn sóng (wavy microchannels) nhằm tăng cường hơn nữa chỉ số Nusselt.
  • Nghiên cứu chu trình kết hợp bơm nhiệt $\text{CO}_2$ đa năng (vừa làm mát vừa cấp nước nóng gia dụng) ứng dụng dàn bay hơi kênh micro trong các tòa nhà xanh.

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

  • Tác động học thuật sâu rộng: Công trình đã công bố 9 bài báo khoa học, bao gồm 01 bài báo trên tạp chí SCIE (2022) có Impact Factor 3.5, 01 bài báo trên tạp chí WoS/ESCI (Q3), và 02 chương sách chuyên khảo (Book Chapters) thuộc hệ thống Scopus từ kỷ yếu hội nghị quốc tế uy tín. Đây là minh chứng mạnh mẽ cho chất lượng học thuật và giá trị trích dẫn quốc tế của đề tài.
  • Thúc đẩy đổi mới công nghiệp: Cung cấp giải pháp kỹ thuật trực tiếp cho các nhà sản xuất máy lạnh dân dụng và điều hòa không khí phương tiện giao thông tại Việt Nam trong việc chuyển đổi công nghệ sang môi chất R744, giảm kích thước dàn trao đổi nhiệt từ 40 – 60% nhờ công nghệ vi kênh.
  • Lợi ích xã hội và môi trường: Góp phần cắt giảm hàng triệu tấn khí thải nhà kính tương đương $\text{CO}_2$ nếu công nghệ này được thương mại hóa thay thế cho các dòng máy lạnh dùng R410A hay R32, bảo vệ tầng ozone và đáp ứng các cam kết Net Zero của Việt Nam tại COP26.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận mô hình giải tích cải tiến cho hệ số ma sát vi kênh và bộ dữ liệu thực nghiệm chuẩn hóa về dòng sôi hai pha của $\text{CO}_2$ để phát triển các đề tài chuyên sâu.
  • Kỹ sư R&D ngành Nhiệt - Điện lạnh: Sở hữu tài liệu hướng dẫn thiết kế chi tiết, từ phương pháp tính toán EES, thông số phân chia lưới COMSOL, đến kinh nghiệm xử lý an toàn áp suất cao (74 – 90 bar) cho hệ thống điều hòa cỡ nhỏ.
  • Doanh nghiệp sản xuất thiết bị lạnh: Tận dụng thiết kế dàn bay hơi nhôm kênh micro để giảm lượng nạp gas, thu nhỏ kích thước máy, giảm chi phí nguyên vật liệu đồng và nâng cao sức cạnh tranh trên thị trường xuất khẩu.
  • Cơ quan quản lý nhà nước: Sử dụng các kết quả kiểm định làm cơ sở khoa học để ban hành quy chuẩn kỹ thuật an toàn cho thiết bị lạnh sử dụng môi chất $\text{CO}_2$ tự nhiên tại Việt Nam.

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 hiệu chỉnh mô hình truyền nhiệt sôi hai pha của Choi et al. (2007) thông qua việc thay thế hệ số ma sát đồng nhất Blasius bằng mô hình phân tách chế độ chảy có xét đến tỉ số cạnh hình học $a_c$ của kênh micro chữ nhật. Điều này cho phép tính toán chính xác tham số Martinelli $X$ và dự báo phân bố hệ số $\alpha_{r,tp}$ dọc theo dải độ khô $x \in [0,1; 1,0]$ với sai số chỉ $\pm 1,5\text{ kW/m}^2\text{K}$.

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

So với nghiên cứu của Zhao et al. (2000) (chỉ đo trên ống tròn đơn giản lẻ) và Yun et al. (2005) (chỉ mô phỏng kênh ngắn), luận án đã thiết lập quy trình kiểm chứng tam giác hóa đồng thời (Analytical EES - CFD COMSOL v6.0 RANS $k-\epsilon$ - Thực nghiệm hệ thống kín hoàn chỉnh) trên dàn bay hơi nhôm 10 kênh micro có cánh louver, kiểm soát đồng thời cả hiện tượng truyền nhiệt truyền chất phía không khí ẩm.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và cơ chế nhiệt động học đằng sau kết quả đó?

Phát hiện bất ngờ nhất là thiết bị hồi nhiệt IHX giúp tăng COP tới 103% và công suất lạnh tăng 100%, vượt xa mức tăng 20-25% thường thấy ở các hệ thống dưới tới hạn truyền thống. Cơ chế là do ở chu trình $\text{CO}2$ trên tới hạn, đường đẳng enthalpy trong vùng lỏng quá lạnh áp suất cao có độ dốc rất lớn; việc hạ sâu nhiệt độ trước tiết lưu bằng IHX triệt tiêu hiện tượng giãn nở đoạn nhiệt sinh hơi sớm, làm giảm đột ngột độ khô đầu vào $x{in}$ từ 0,48 xuống 0,22, giải phóng toàn bộ năng suất lạnh tiềm năng của dàn bay hơi.

4. Luận án có cung cấp quy trình nhân rộng và tái lập thí nghiệm (Replication Protocol) không?

Có. Toàn bộ thông số hình học dàn Danfoss ($D_h = 0,9\text{ mm}$), thông số máy nén Dorin CD 180H, mã nguồn chương trình tính toán nhiệt trên EES (Phụ lục 1), bảng tính Excel thiết kế chi tiết (Phụ lục 4), thiết lập điều kiện biên và chia lưới COMSOL (Phụ lục 5) cùng toàn bộ dữ liệu thực nghiệm thô (Phụ lục 6) đều được cung cấp chi tiết và minh bạch trong văn bản luận án.

5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Trọng tâm 10 năm tới chuyển dịch vào: (1) Tích hợp hệ thống quản lý nhiệt thông minh cho xe điện (EV Thermal Management System) sử dụng bơm nhiệt vi kênh $\text{CO}_2$ hoạt động đa chế độ từ $-20^\circ\text{C}$ đến $45^\circ\text{C}$; (2) Tối ưu hóa vi cấu trúc bề mặt kênh micro bằng công nghệ phủ nano để chống bám dầu và tăng cường mật độ điểm tạo mầm sôi; (3) Thương mại hóa máy điều hòa gia dụng R744 hiệu suất cao tại thị trường nhiệt đới gió mùa.

Kết luận

  1. Thiết kế và chế tạo thành công hệ thống máy điều hòa không khí cỡ nhỏ ($\le 9000\text{ BTU/h}$) sử dụng môi chất lạnh tự nhiên $\text{CO}_2$ (R744) với thiết bị bay hơi kênh micro nhôm 10 cổng ($D_h = 0,9\text{ mm}$), vận hành tuyệt đối an toàn và ổn định ở dải áp suất trên tới hạn thấp từ 74 đến 90 bar.
  2. Xác lập mô hình lý thuyết cải tiến dự báo chính xác hệ số tỏa nhiệt đối lưu hai pha và tổn thất áp suất trên cơ sở hiệu chỉnh hệ số ma sát theo tỉ số kích thước mặt cắt kênh micro $a_c$, giải quyết triệt để sự thiếu chính xác của công thức Blasius kinh điển.
  3. Chứng minh độ chính xác vượt trội của mô phỏng số COMSOL v6.0 với sai số độ khô cực đại so với thực nghiệm chỉ 5,5%, sai số tính toán lý thuyết chỉ 3,5%, và tổn thất áp suất thực nghiệm 1,5 bar so với mô phỏng 1,4 bar.
  4. Khám phá giải pháp nâng cao hiệu suất nhảy vọt: Tích hợp thiết bị hồi nhiệt IHX làm tăng năng suất lạnh 100% và COP tăng 103%; tích hợp Subcooler làm tăng năng suất lạnh 50% và COP tăng 39%.
  5. Xác định thông số vận hành tối ưu toàn diện: Điểm làm việc lý tưởng tại lưu lượng khối $\text{CO}_2 = 111\text{ kg/h}$ cho công suất lạnh cực đại 3,12 kW và $\text{COP} = 3,15$; vận tốc không khí $v_a = 5,1\text{ m/s}$ thiết lập sự cân bằng truyền nhiệt hoàn hảo giữa hai môi chất ở mức 2,09 kW.
  6. Mở ra 3 hướng nghiên cứu chuyên sâu mới: Bơm nhiệt vi kênh $\text{CO}_2$ cho ô tô điện, điều khiển tối ưu EEV tự thích nghi theo áp suất môi trường, và công nghệ kênh micro biến tính bề mặt nano.
  7. Đóng góp quốc tế uy tín: 9 công trình khoa học công bố, trong đó có 01 bài báo SCIE (IF 3.5), khẳng định vị thế tiên phong của nghiên cứu nhiệt - lạnh Việt Nam trên bản đồ học thuật thế giới.