Tổng quan về luận án
Sự phát triển của ngành công nghiệp nhiệt - điện lạnh toàn cầu đang đứng trước bước ngoặt chuyển dịch công nghệ sâu sắc dưới tác động của Nghị định thư Montreal và Bản sửa đổi, bổ sung Kigali. Việc loại bỏ các môi chất lạnh nhân tạo nhóm HCFC và HFC vốn có tiềm năng gây nóng lên toàn cầu (GWP) cao đang thúc đẩy mạnh mẽ các nghiên cứu về môi chất tự nhiên. Trong bối cảnh đó, Carbon Dioxide ($CO_2$ - R744) nổi lên như một giải pháp thay thế tối ưu với các đặc tính thân thiện môi trường vượt trội ($ODP = 0$, $GWP = 1$), không độc hại, không cháy và có sẵn trong tự nhiên. Tuy nhiên, tính chất nhiệt động của chu trình $CO_2$ trên tới hạn (transcritical cycle) đặt ra thách thức lớn: áp suất làm việc rất cao (vượt qua áp suất tới hạn $73,8\text{ bar}$) và hiệu suất năng lượng (COP) suy giảm đáng kể khi nhiệt độ môi trường giải nhiệt tăng cao. Song song đó, các thiết bị trao đổi nhiệt dạng ống - cánh truyền thống thường cồng kềnh và có hiệu quả truyền nhiệt hữu hạn. Việc ứng dụng công nghệ kênh micro với đường kính thủy lực tương đương nhỏ ($D_h \le 1\text{ mm}$) để chế tạo thiết bị bay hơi compact là hướng đi đột phá nhằm tối ưu hóa mật độ dòng nhiệt và giảm kích thước thiết bị.
Luận án tiến sĩ kỹ thuật cơ khí mang tên "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 $CO_2$" do nghiên cứu sinh Nguyễn Trọng Hiếu thực hiện dưới sự hướng dẫn khoa học của PGS.TS Đặng Thành Trung và GS.TS Jyh Tong Teng tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh đã giải quyết trọn vẹn khoảng trống học thuật này. Nghiên cứu tập trung giải đáp hai câu hỏi nghiên cứu cốt lõi:
- Cơ chế truyền nhiệt hai pha khi sôi và tổn thất áp suất của $CO_2$ diễn biến như thế nào trong ống dẹt đa kênh micro dưới các dải thông số vận hành thực tế?
- Những giải pháp nhiệt động nào giúp tối ưu hóa hệ số COP và công suất lạnh của chu trình $CO_2$ trên tới hạn cỡ nhỏ khi vận hành ở dải áp suất đầu đẩy trung bình ($74 - 90\text{ bar}$)?
Trên cơ sở nền tảng của Nhiệt động lực học kỹ thuật và Truyền nhiệt - Truyền chất đa pha, luận án thiết lập phạm vi nghiên cứu tập trung vào hệ thống điều hòa không khí cỡ nhỏ có công suất danh định $9000\text{ BTU/h}$ ($\approx 2600\text{ W}$), nhiệt độ môi chất bay hơi $T_e = 5 - 15^\circ\text{C}$ (tương ứng áp suất bay hơi $40 - 50\text{ bar}$), áp suất thiết bị làm mát $74 - 90\text{ bar}$, và thiết bị bay hơi bằng nhôm dạng ống dẹt 10 kênh micro song song với đường kính thủy lực $D_h = 0,9\text{ mm}$. Nghiên cứu đã tạo ra đóng góp mang tính định lượng rõ rệt: chứng minh việc bổ sung thiết bị làm mát phụ (Subcooler) giúp hạ nhiệt độ trước van tiết lưu $1,4^\circ\text{C}$, làm tăng năng suất lạnh $50%$ và nâng hệ số COP thêm $39%$; đặc biệt, tích hợp thiết bị hồi nhiệt (Internal Heat Exchanger - IHX) giúp năng suất lạnh tăng tới $100%$ và COP tăng vượt bậc $103%$.
Literature Review và Positioning
Tổng quan y văn quốc tế ghi nhận ba dòng nghiên cứu chính định hình lĩnh vực trao đổi nhiệt $CO_2$ trên tới hạn: nghiên cứu chu trình nhiệt động nâng cao hiệu suất, đặc tính truyền nhiệt sôi trong kênh thu nhỏ, và động học dòng chảy đa pha. Về phương diện chu trình, Kim và cộng sự [2] đã chỉ ra nguyên nhân chu trình $CO_2$ trên tới hạn cơ bản có hiệu quả năng lượng thấp hơn chu trình R134a truyền thống là do tổn thất giải nhiệt ở bộ làm mát khí (gas cooler) và tổn thất tiết lưu sinh entropy lớn. Để bù đắp khiếm khuyết này, Boewe và cộng sự [13] chứng minh thiết bị hồi nhiệt IHX có thể nâng COP lên đến $25%$ trên hệ thống $CO_2$ nhỏ. Song song đó, các nghiên cứu của Sanchez và cộng sự [7] hay Llopis và cộng sự [9] ứng dụng công nghệ làm mát cơ học chuyên dụng (Dedicated Mechanical Subcooling - DMS) giúp cải thiện năng suất lạnh từ $12,7%$ đến $20,7%$.
Về đặc tính truyền nhiệt kênh micro, định nghĩa và phân loại kênh đã trải qua nhiều tranh luận học thuật giữa Kandlikar & Grande [18] ($10 \le D \le 200\text{ }\mu\text{m}$), Mehendale và cộng sự [19] ($1 \le D \le 100\text{ }\mu\text{m}$), và các học giả chấp nhận dải kênh micro mở rộng $D \le 1000\text{ }\mu\text{m}$ [20, 21]. Trong dòng chảy sôi hai pha, Choi và cộng sự [68] đã đề xuất mô hình tương quan tổng quát cho hệ số tỏa nhiệt đối lưu, trong khi Yun và cộng sự [41] khẳng định hệ số tỏa nhiệt đối lưu của $CO_2$ trong kênh micro cao hơn $53%$ so với R134a. Về tổn thất áp suất, Pamitran và cộng sự [34] chỉ ra tổn thất áp suất của $CO_2$ trong ống $1,5\text{ mm}$ thấp hơn $16%$ so với R-22, và Zhao và cộng sự [35] chứng minh trong ống $0,86\text{ mm}$, tổn thất áp suất $CO_2$ thấp hơn $40%$ so với R134a nhờ khối lượng riêng pha hơi cao và độ nhớt động học thấp của R744.
Trong y văn tồn tại hai cuộc tranh luận học thuật lớn. Thứ nhất là tính khả thi giữa công nghệ Flash Gas Bypass và IHX trong hệ thống kênh micro: Kim và cộng sự [2] cho rằng Flash Gas Bypass gây khó khăn nghiêm trọng trong việc phân chia dòng môi chất hai pha đồng đều vào các kênh micro, trong khi Beaver và cộng sự [17] nhấn mạnh sự bất ổn định của độ quá nhiệt đầu ra. Thứ hai là giả định dòng đồng nhất (homogeneous model) khi tính hệ số ma sát hai pha của Choi và cộng sự [68].
Luận án của Nguyễn Trọng Hiếu định vị chính xác vào khoảng trống này: đa số các công bố quốc tế như Cho & Kim [48], Oh & Son [51] ($D = 4,57 - 7,75\text{ mm}$), hay Ducoulombier và cộng sự [53] ($D = 0,529\text{ mm}$) chỉ khảo sát ống đơn hoặc ống thép không gỉ trần ở dải áp suất trên $90\text{ bar}$. Nghiên cứu này tiên phong tiếp cận cấu trúc dàn bay hơi ống dẹt 10 kênh micro bằng nhôm có cánh tản nhiệt ($D_h = 0,9\text{ mm}$) trong điều hòa không khí dân dụng, đồng thời giải quyết bài toán vận hành ở áp suất đầu đẩy an toàn $74 - 90\text{ bar}$ với sự kiểm chứng đa chiều giữa lý thuyết, mô phỏng số 3D và thực nghiệm.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đã đóng góp vào lý thuyết truyền nhiệt sôi hai pha bằng cách hiệu chỉnh và mở rộng mô hình truyền nhiệt đối lưu của Choi và cộng sự [68]. Mô hình gốc của Choi giả định hệ số ma sát dòng chảy tuân theo công thức Blasius đồng nhất cho cả hai pha lỏng và hơi:
$$f = 0,079 \cdot Re^{-0,25}$$
Nhận thấy hạn chế của giả định này đối với kênh micro tiết diện phi tròn có tỷ lệ cạnh hẹp, tác giả đã đề xuất cấu trúc phân tách hệ số ma sát riêng biệt theo chế độ chảy thực của từng pha dựa trên lý thuyết dòng chảy trong vi kênh:
- Với dòng chảy tầng hoặc quá độ ($Re < 4000$), áp dụng nghiệm đa thức bậc 5 phụ thuộc vào tỷ số cạnh $a_c$ của kênh micro:
$$Po = f \cdot Re = 24 \cdot \left(1 - 1,3553a_c + 1,9467a_c^2 - 1,7012a_c^3 + 0,9564a_c^4 - 0,2537a_c^5\right)$$
- Với dòng chảy rối ($Re \ge 4000$), phương trình ma sát được tích hợp độ nhám tương đối $\varepsilon/D_h$:
$$f = \frac{0,25 \cdot \left[\ln\left(\frac{\varepsilon/D_h}{3,7} + \frac{5,74}{Re^{0,9}}\right)\right]^{-2}}{4}$$
Sự hiệu chỉnh này cho phép tham số Martinelli $X$ và hệ số tăng cường sôi hai pha $S, F$ phản ánh chính xác tương tác thủy lực vi mô trong kênh $D_h = 0,9\text{ mm}$, nâng cao độ chính xác dự báo hệ số tỏa nhiệt đối lưu $\alpha_{r,tp}$ từ vùng sôi bọt đến vùng quá nhiệt hơi.
MÔ HÌNH TRUYỀN NHIỆT VÀ THỦY LỰC HIỆU CHỈNH
Khung phân tích độc đáo
Khung phân tích của đề tài được xây dựng trên sự tích hợp của ba trụ cột lý thuyết:
- Lý thuyết Chu trình Nhiệt động trên tới hạn: Phân tích các quá trình nén đẳng entropy, làm mát đẳng áp trên điểm tới hạn, tiết lưu đẳng enthalpy ($h_3 = h_4$) và bay hơi đẳng áp thu nhiệt $Q_e = \dot{m}r(h{out} - h_{in})$.
- Lý thuyết Truyền nhiệt Hai pha và Chuyển pha vi mô: Tách biệt diện tích trao đổi nhiệt thành vùng hai pha ($A_{e1}$) và vùng quá nhiệt ($A_{e2}$), kết hợp hệ số truyền nhiệt tổng $k_{e1}, k_{e2}$ có tính đến nhiệt trở dẫn nhiệt vách nhôm $R_b$ và hệ số làm cánh $\beta_e$.
- Lý thuyết Nhiệt - Ẩm Không khí: Xác định hệ số tỏa nhiệt đối lưu phía không khí $\alpha_{a,dry}$ qua công thức thực nghiệm phụ thuộc bước cánh $S_p$ ($1 - 1,2\text{ mm}$), kết hợp hiệu suất cánh $\eta_f = \tanh(ml)/(ml)$ và hệ số tách ẩm $\xi$ khi không khí bị làm lạnh dưới nhiệt độ đọng sương:
$$\xi = 1 + 2480 \cdot \frac{d_{a,in} - d_{a,out}}{t_{a,in} - t_{a,out}} \cdot \frac{Q_t}{Q_s}$$
Điều kiện biên xác lập chặt chẽ: áp suất đầu đẩy $P_c \in [74, 90]\text{ bar}$, áp suất bay hơi $P_e \in [40, 50]\text{ bar}$, lưu lượng $CO_2$ $\dot{m} \in [75,6; 121,4]\text{ kg/h}$, vận tốc không khí qua 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
Nghiên cứu tuân thủ thế giới quan thực chứng (positivism) với phương pháp tiếp cận đa tầng kết hợp chặt chẽ giữa tính toán giải tích lý thuyết (Analytical Calculation), mô phỏng số động lực học dòng chảy 3D (Computational Fluid Dynamics - CFD) và thực nghiệm kiểm chứng (Experimental Validation) trên hệ thống thực tế.
THIẾT KẾ NGHIÊN CỨU ĐA TẦNG (TRIANGULATION)
Quy trình nghiên cứu rigorous
Quy trình nghiên cứu được triển khai theo các bước chuẩn hóa nghiêm ngặt:
- Mô hình hóa giải tích: Thiết lập thuật toán giải trên phần mềm EES (Engineering Equation Solver) và bảng tính chuyên dụng để xác định các điểm nút nhiệt động chu trình và kích thước hình học dàn bay hơi.
- Mô phỏng số 3D bằng COMSOL Multiphysics v6.0: Thiết lập bài toán đa vật lý (Multiphysics) giải đồng thời phương trình dòng chảy rối RANS $k-\varepsilon$ (Reynolds-Averaged Navier-Stokes), phương trình bảo toàn năng lượng và truyền chất chuyển pha. Mô hình hình học tái lập chính xác ống dẹt đa kênh micro Danfoss với 10 đường dẫn song song. Lưới phần tử hữu hạn được khảo sát tính độc lập với 5 cấp độ lưới tự động kết hợp chia lưới thủ công tại lớp biên thủy lực và lớp biên nhiệt, kiểm tra độ hội tụ nghiêm ngặt của số dư phân giải.
- Hệ thống thực nghiệm thực tế: Hệ thống được chế tạo hoàn chỉnh 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. Cấu hình bao gồm:
- Máy nén bán kín chuyên dụng cho $CO_2$ transcritical: Dorin CD 180H.
- Thiết bị làm mát (Gas cooler): Panasonic.
- Thiết bị bay hơi kênh micro (MCHX): Danfoss (vật liệu nhôm, ống dẹt $D_h = 0,9\text{ mm}$).
- Thiết bị hồi nhiệt (IHX) và thiết bị làm mát phụ (Subcooler).
- Hệ thống van tiết lưu tự động, van an toàn bảo vệ dải áp suất cao $74 - 90\text{ bar}$.
Data và phân tích
Toàn bộ thiết bị đo lường được kiểm định nghiêm ngặt: cảm biến áp suất dải đo $0 - 160\text{ bar}$, đầu đo nhiệt độ RTD PT100 độ chính xác cao, lưu lượng kế khối lượng kiểu Coriolis và máy đo vận tốc gió màng nhiệt.
Dữ liệu thực nghiệm được đối soát chéo (cross-validation) với mô phỏng và lý thuyết. Mức độ tương thích và độ tin cậy được thể hiện qua các chỉ số định lượng cụ thể trích xuất trực tiếp từ luận án:
- Độ khô của môi chất ($x$): Sai số độ khô cực đại giữa mô phỏng số COMSOL 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ỉ là $3,5%$.
- Tổn thất áp suất ($\Delta p$): Tổn thất áp suất tổng trên dàn bay hơi thu được từ ba phương pháp lần lượt là: tính toán lý thuyết đạt $1,13\text{ bar}$, mô phỏng số đạt $1,40\text{ bar}$, và đo đạc thực nghiệm đạt $1,50\text{ bar}$.
- Hệ số tỏa nhiệt đối lưu hai pha ($\alpha_{r,tp}$): Giá trị tính toán và mô phỏng số nằm hoàn toàn trong dải biến thiên thực nghiệm từ $6,5\text{ kW/m}^2\text{K}$ (vùng hai pha đầu vào) giảm dần xuống $1,3\text{ kW/m}^2\text{K}$ (vùng hơi quá nhiệt đầu ra), với biên độ sai số tuyệt đối trong khoảng $\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
- Hiệu ứng nhiệt động đột phá từ Thiết bị làm mát phụ (Subcooler): Việc lắp đặt bộ làm mát phụ trước van tiết lưu giúp hạ nhiệt độ dòng $CO_2$ lỏng thêm $1,4^\circ\text{C}$. Độ giảm nhiệt độ này trực tiếp làm giảm độ khô đầu vào dàn bay hơi sau tiết lưu, giải phóng enthalpy hữu ích và giúp năng suất lạnh tăng $50%$, đồng thời cải thiện hệ số COP thêm $39%$.
- Bước nhảy hiệu suất từ Thiết bị hồi nhiệt (IHX): Tích hợp IHX cho phép dòng gas lạnh đầu ra dàn bay hơi làm mát sâu dòng gas nóng sau thiết bị làm mát. Kết quả thực nghiệm ghi nhận sự gia tăng vượt bậc: năng suất lạnh tăng $100%$ và hệ số COP tăng tới $103%$, đồng thời triệt tiêu hoàn toàn nguy cơ lọt lỏng về máy nén Dorin CD 180H.
- Quy luật đáp ứng của Lưu lượng khối lượng $CO_2$: Khi điều chỉnh lưu lượng $CO_2$ từ $97,5\text{ kg/h}$ lên $121,4\text{ kg/h}$, nhiệt độ bay hơi tăng tuyến tính từ $8,2^\circ\text{C}$ lên $14,5^\circ\text{C}$. Hệ thống đạt điểm cực trị tối ưu tại lưu lượng $\dot{m} = 111\text{ kg/h}$, ghi nhận công suất lạnh đạt đỉnh $3,12\text{ kW}$ và hệ số COP đạt giá trị cao nhất $3,15$.
- Điểm cân bằng nhiệt động phía không khí: Khảo sát dải vận tốc không khí thổi qua dàn bay hơi từ $v_a = 0,5\text{ m/s}$ đến $5,1\text{ m/s}$ trong điều kiện lưu lượng $CO_2$ cố định $75,6\text{ kg/h}$. Kết quả xác định tại vận tốc $v_a = 5,1\text{ m/s}$, công suất thu nhiệt phía không khí cân bằng chính xác với năng suất lạnh hấp thụ phía môi chất $CO_2$ ở mức $2,09\text{ kW}$.
- Cơ chế suy giảm hệ số tỏa nhiệt đối lưu theo độ khô: Quá trình mô phỏng và thực nghiệm chứng minh rằng hệ số tỏa nhiệt $\alpha_{r,tp}$ đạt cực đại ở vùng sôi bọt ẩm ($x \approx 0,2 - 0,4$) đạt xấp xỉ $6,5\text{ kW/m}^2\text{K}$, nhưng suy giảm mạnh khi độ khô $x \to 1$ và chuyển sang vùng quá nhiệt chỉ còn $1,3\text{ kW/m}^2\text{K}$ do hiện tượng cạn màng chất lỏng (dryout) sát vách vi kênh.
Implications đa chiều
- Về mặt học thuật: Luận án cung cấp bộ dữ liệu thực nghiệm và phân bố trường nhiệt độ - áp suất có giá trị tham chiếu cao cho môi chất $CO_2$ trong cấu trúc ống dẹt đa kênh micro, mở rộng hiểu biết về hiện tượng sôi cưỡng bức trong không gian hẹp dưới áp suất cận tới hạn.
- Về phương pháp luận: Xác lập quy trình mô hình hóa liên hoàn từ tính toán nhiệt động cơ sở đến mô phỏng số 3D đa vật lý trên COMSOL Multiphysics, tạo chuẩn mực cho việc thiết kế bộ trao đổi nhiệt compact chịu áp lực cao.
- Về ứng dụng thực tiễn: Chứng minh tính khả thi kỹ thuật của việc ứng dụng $CO_2$ cho máy điều hòa không khí dân dụng và thương mại cỡ nhỏ ($9000\text{ BTU/h}$) hoạt động an toàn ở áp suất $74 - 90\text{ bar}$ mà vẫn đạt COP cạnh tranh ($> 3,0$) nhờ tích hợp IHX và Subcooler.
- Về chính sách và môi trường: Cung cấp cơ sở khoa học vững chắc hỗ trợ các nhà hoạch định chính sách và doanh nghiệp HVAC xây dựng lộ trình loại bỏ hoàn toàn môi chất HFCs theo cam kết giảm phát thải khí nhà kính quốc gia.
Limitations và Future Research
Mặc dù đạt được những kết quả đột phá, công trình nghiên cứu vẫn tồn tại một số giới hạn khoa học cần được tiếp tục hoàn thiện:
- Giới hạn về thủy động lực học màng chất lỏng: Luận án tập trung chuyên sâu vào các đặc tính truyền nhiệt (nhiệt độ, áp suất, hệ số đối lưu, mật độ dòng nhiệt) mà chưa đi sâu khảo sát các đặc tính dòng vi mô chi tiết như lý thuyết lớp biên động học, độ nhám bề mặt thực tế bên trong lòng kênh và hiện tượng phân bố không đều dòng hai pha giữa các ống nhánh từ ống góp (header maldistribution).
- Phạm vi điều kiện biên môi trường: Thử nghiệm thực nghiệm chủ yếu thực hiện trong dải nhiệt độ môi trường xung quanh $28 - 32^\circ\text{C}$ tại khu vực TP. Hồ Chí Minh, chưa mở rộng sang các điều kiện khí hậu khắc nghiệt ($> 40^\circ\text{C}$ hoặc $< 0^\circ\text{C}$).
- Đặc tính quá độ và điều khiển tự động: Nghiên cứu khảo sát hệ thống ở trạng thái làm việc xác lập (steady-state), chưa phân tích động học đáp ứng quá độ khi tải lạnh thay đổi đột ngột.
Từ những hạn chế trên, các hướng phát triển nghiên cứu tiếp theo bao gồm:
- Hướng 1: Nghiên cứu tối ưu hóa hình học ống góp (manifold/header) nhằm triệt tiêu hiện tượng phân bố lệch pha môi chất vào các kênh micro song song.
- Hướng 2: Tích hợp công nghệ van tiết lưu điện tử (EEV) kết hợp thuật toán điều khiển mờ (Fuzzy Logic) hoặc AI để tự động tối ưu hóa áp suất đầu đẩy $P_c$ theo nhiệt độ môi trường ngoài trời.
- Hướng 3: Khảo sát ứng dụng bộ làm mát phụ kết hợp năng lượng địa nhiệt hoặc hệ thống bay hơi gián tiếp phục vụ các trạm sạc xe điện và điều hòa không khí ô tô thế hệ mới.
- Hướng 4: Nghiên cứu chuyên sâu ảnh hưởng của dầu bôi trơn chuyên dụng hòa tan trong $CO_2$ đến hệ số truyền nhiệt sôi trong vi kênh.
Tác động và ảnh hưởng
Công trình nghiên cứu đã tạo ra những tác động học thuật và thực tiễn sâu rộng:
- Giá trị công bố quốc tế: Kết quả của luận án đã được công bố qua 9 công trình khoa học uy tín, bao gồm 01 bài báo trên tạp chí SCIE (Impact Factor 3.5), 01 bài báo trên tạp chí WoS-ESCI (Scimago Q3), và 02 chương sách (Book Chapters) thuộc hệ thống Scopus xuất bản từ kỷ yếu hội nghị quốc tế chuyên ngành.
- Thúc đẩy chuyển giao công nghệ công nghiệp: Đồ án mở đường cho các doanh nghiệp sản xuất thiết bị lạnh tại Việt Nam và khu vực tiếp cận công nghệ chế tạo dàn trao đổi nhiệt microchannel nhôm gọn nhẹ, giảm lượng nạp môi chất đến $40 - 60%$ so với dàn ống đồng truyền thống.
- Đóng góp cho phát triển bền vững: Góp phần hiện thực hóa mục tiêu giảm phát thải ròng bằng 0 (Net Zero) vào năm 2050 thông qua việc ứng dụng môi chất tự nhiên $CO_2$ ($GWP = 1$), loại bỏ nguy cơ phát thải khí nhà kính từ các môi chất flo nhân tạo.
Đối tượng hưởng lợi
MA TRẬN ĐỐI TƯỢNG HƯỞNG LỢI
- Nghiên cứu sinh và Giới nghiên cứu học thuật: Tiếp cận phương pháp luận chuẩn mực kết hợp mô phỏng số COMSOL Multiphysics với thực nghiệm nhiệt động chính xác; kế thừa công thức hiệu chỉnh hệ số ma sát và truyền nhiệt sôi hai pha trong vi kênh.
- Kỹ sư R&D ngành Nhiệt - Lạnh và Ô tô: Nhận được các thông số thiết kế thực tế để chế tạo hệ thống điều hòa $CO_2$ cỡ nhỏ cho cabin xe điện, kho lạnh mini và điều hòa dân dụng đạt hiệu suất tối ưu ở dải áp suất $74 - 90\text{ bar}$.
- Các nhà quản lý năng lượng và môi trường: Có cơ sở dữ liệu định lượng để xây dựng quy chuẩn kỹ thuật quốc gia về thử nghiệm an toàn và đánh giá hiệu suất hệ thống lạnh sử dụng môi chất tự nhiên áp suất cao.
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 lý thuyết cốt lõi là việc hiệu chỉnh công thức xác định hệ số ma sát hai pha $f$ trong mô hình truyền nhiệt sôi của Choi và cộng sự [68]. Bằng cách thay thế giả định dòng đồng nhất Blasius bằng mô hình phân tách trạng thái chảy (đa thức bậc 5 theo tỷ số cạnh kênh $a_c$ cho vùng chảy tầng/quá độ $Re < 4000$ và mô hình có xét độ nhám tương đối cho vùng chảy rối $Re \ge 4000$), nghiên cứu đã nâng cao độ chính xác dự báo hệ số tỏa nhiệt đối lưu $\alpha_{r,tp}$ với sai số thực nghiệm được kiểm soát chặt chẽ trong biên độ $\pm 1,5\text{ kW/m}^2\text{K}$.
2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm?
So với các nghiên cứu của Cho & Kim [48] hay Oh & Son [51] vốn chỉ tập trung vào ống đơn bằng thép không gỉ trần, luận án đã thiết lập quy trình kiểm chứng chéo 3 chiều đồng thời: Giải tích nhiệt động (EES/Excel) $\leftrightarrow$ Mô phỏng số 3D đa vật lý RANS $k-\varepsilon$ (COMSOL v6.0) $\leftrightarrow$ Thực nghiệm Test-bench hoàn chỉnh trên dàn bay hơi nhôm ống dẹt 10 kênh micro Danfoss ($D_h = 0,9\text{ mm}$).
3. Phát hiện thực nghiệm nào mang tính bất ngờ nhất?
Kết quả bất ngờ nhất là mức độ khuếch đại hiệu suất nhiệt động vượt trội của thiết bị hồi nhiệt IHX: chỉ bằng việc trao đổi nhiệt nội bộ giữa dòng gas hút về máy nén và dòng gas lỏng sau thiết bị làm mát, công suất lạnh đã tăng $100%$ và hệ số COP tăng tới $103%$, cao hơn rất nhiều so với mức tăng $25%$ từng được công bố trong các nghiên cứu lý thuyết trước đây của Boewe và cộng sự [13].
4. Khả năng tái lập (Replication Protocol) của nghiên cứu được đảm bảo như thế nào?
Nghiên cứu công bố chi tiết toàn bộ thông số hình học kênh micro ($D_h = 0,9\text{ mm}$, bước cánh $S_p = 1 - 1,2\text{ mm}$), mã hiệu thiết bị chuẩn hóa (máy nén Dorin CD 180H, dàn Danfoss, gas cooler Panasonic), thuật toán phân chia lưới phần tử hữu hạn, các phương trình điều khiển trong phụ lục EES và bảng dữ liệu thực nghiệm chi tiết.
5. Tầm nhìn chương trình nghiên cứu 10 năm tiếp theo?
Chương trình 10 năm hướng tới việc hoàn thiện công nghệ chu trình $CO_2$ siêu tới hạn tích hợp Ejector điều khiển điện tử đa cấp, tối ưu hóa vật liệu phủ nano chống bám tuyết trên cánh tản nhiệt micro, và thương mại hóa module điều hòa nhiệt độ tích hợp cho toàn bộ các dòng xe điện thế hệ mới tại Đông Nam Á.
Kết luận
Công trình nghiên cứu của NCS. Nguyễn Trọng Hiếu đã hoàn thành xuất sắc các mục tiêu khoa học đề ra, đúc kết thành 5 đóng góp học thuật và thực tiễn cốt lõi:
- Xác lập thành công bộ dữ liệu thực nghiệm toàn diện về đặc tính truyền nhiệt và tổn thất áp suất của môi chất $CO_2$ trong thiết bị bay hơi ống dẹt đa kênh micro ($D_h = 0,9\text{ mm}$) ở dải nhiệt độ bay hơi $5 - 15^\circ\text{C}$ và áp suất đầu đẩy $74 - 90\text{ bar}$.
- Hiệu chỉnh hoàn thiện mô hình toán học tính toán hệ số ma sát và truyền nhiệt sôi hai pha, đạt độ tin cậy vượt trội với sai số độ khô cực đại giữa mô phỏng và thực nghiệm chỉ $5,5%$, giữa lý thuyết và thực nghiệm chỉ $3,5%$.
- Chứng minh định lượng giải pháp nâng cao hiệu năng đột phá: thiết bị làm mát phụ giúp tăng $50%$ năng suất lạnh và $39%$ COP; thiết bị hồi nhiệt IHX giúp tăng $100%$ năng suất lạnh và $103%$ COP.
- Xác định chính xác chế độ vận hành tối ưu của hệ thống tại lưu lượng $CO_2$ $\dot{m} = 111\text{ kg/h}$ (đạt công suất lạnh $3,12\text{ kW}$, COP đạt $3,15$) và vận tốc không khí qua dàn $v_a = 5,1\text{ m/s}$ (công suất trao đổi nhiệt đạt $2,09\text{ kW}$).
- Mở ra 3 hướng nghiên cứu ứng dụng mũi nhọn: tối ưu hóa phân bố dòng hai pha trong ống góp micro, điều khiển thông minh chu trình $CO_2$ transcritical, và phát triển hệ thống điều hòa không khí xanh cho phương tiện giao thông điện hóa toàn cầu.