Tổng quan về luận án

Nghiên cứu về máy lạnh hấp thụ $NH_3\text{-}H_2O$ chu trình liên tục phục vụ sản xuất nước đá và làm lạnh sâu đóng vai trò quan trọng trong chiến lược chuyển dịch năng lượng và bảo vệ môi trường toàn cầu. Khác với môi chất $H_2O\text{-}LiBr$ bị giới hạn ở dải nhiệt độ trên $0^\circ\text{C}$, cặp môi chất $NH_3\text{-}H_2O$ cho phép đạt nhiệt độ bay hơi âm sâu (điểm sôi của $NH_3$ ở $1\text{ bar}$ là $-33,4^\circ\text{C}$, nhiệt độ đông đặc $-77,7^\circ\text{C}$), đáp ứng hoàn hảo yêu cầu bảo quản thủy hải sản và chế biến thực phẩm tại các vùng nhiệt đới ven biển. Trong hệ thống lạnh hấp thụ, bình hấp thụ (absorber) là thiết bị phức tạp nhất, chiếm hơn 40% tổng diện tích trao đổi nhiệt toàn chu trình và quyết định trực tiếp đến hiệu suất năng lượng ($COP$) cũng như giá thành chế tạo.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các công trình quốc tế thường tập trung vào cấu trúc vi kênh (micro-channel) đắt đỏ hoặc bình hấp thụ dạng tấm/dạng bọt đòi hỏi công nghệ chế tạo tinh vi vượt quá điều kiện hạ tầng kỹ thuật tại các nước đang phát triển. Đồng thời, các mô hình truyền nhiệt - truyền chất đồng thời trên màng dung dịch chảy ngang (falling film on horizontal tubes) chưa phản ánh chính xác các điều kiện nhiệt động thực tế tại vùng khí hậu nóng ẩm (nhiệt độ nước giải nhiệt dao động $30\text{--}36^\circ\text{C}$, nhiệt độ nước muối làm đá từ $-10^\circ\text{C}$ đến $-19^\circ\text{C}$). Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Hiếu Nghĩa, thực hiện dưới sự hướng dẫn khoa học của GS. Lê Chí Hiệp và PGS. Hoàng An Quốc tại Trường Đại học Bách Khoa – ĐHQG TP. Hồ Chí Minh (2017), đã giải quyết triệt để bài toán khoa học và công nghệ này.

Hệ thống câu hỏi nghiên cứu và giả thuyết được xác lập cụ thể:

  • Câu hỏi nghiên cứu 1 ($RQ_1$): Cơ chế tương tác động học và quy luật biến thiên cục bộ của hệ số truyền nhiệt ($k$) và hệ số truyền chất ($h_m$) trên màng dung dịch lỏng $NH_3\text{-}H_2O$ bao quanh ống tròn nằm ngang diễn ra như thế nào?
  • Câu hỏi nghiên cứu 2 ($RQ_2$): Đâu là nồng độ dung dịch nạp ban đầu và chế độ vận hành tối ưu nhằm cân bằng giữa tổn thất nhiệt bình phát sinh và tốc độ hấp thụ hơi để cực đại hóa $COP$ và hiệu suất làm lạnh nước muối ($COP_u$)?
  • Câu hỏi nghiên cứu 3 ($RQ_3$): Làm thế nào để thiết lập một tương quan thực nghiệm - số trị đáng tin cậy về $k$ và $h_m$ phục vụ thiết kế chế tạo bình hấp thụ kiểu chùm ống ngang với công nghệ gia công nội địa?

Hệ thống giả thuyết nghiên cứu tương ứng:

  • Giả thuyết 1 ($H_1$): Sự phân bố của hệ số truyền nhiệt và truyền chất dọc theo nửa chu vi ống không đồng đều; hệ số truyền chất đạt cực đại ngay tại điểm tiếp xúc ban đầu của màng và suy giảm nhanh chóng do sự tích tụ điện trở nhiệt và nồng độ ở bề mặt phân chia pha lỏng - hơi.
  • Giả thuyết 2 ($H_2$): Tồn tại một nồng độ nạp tới hạn ($31%$) giúp tối ưu hóa hiệu quả chu trình làm đá, cho phép sai số giữa tính toán mô phỏng nhiệt động và vận hành thực nghiệm đạt dưới $2%$.
  • Giả thuyết 3 ($H_3$): Xử lý cơ học bề mặt ống (thổi cát tăng độ dẫn ướt) kết hợp với cấu trúc phân phối dung dịch nhỏ giọt ($\Phi 1,2\text{ mm}$, bước $4\text{ mm}$) sẽ đảm bảo chế độ màng ổn định, bù đắp sự thiếu hụt công nghệ gia công chuyên dụng mà không làm suy giảm hiệu năng hấp thụ.

Khung lý thuyết của nghiên cứu tích hợp lý thuyết nhiệt động học chất lưu hai cấu tử, lý thuyết chảy màng Nusselt mở rộng, phương trình Navier-Stokes đơn giản hóa cho lớp biên mỏng và định luật khuếch tán Fick kết hợp dẫn nhiệt Fourier. Nghiên cứu mang lại đóng góp đột phá với mô hình thực nghiệm máy lạnh hấp thụ một cấp hoàn chỉnh công suất lạnh $2\text{ kW}$, xác lập bộ thông số vận hành chuẩn: nồng độ nạp tối ưu $31%$, $COP = 0,436$, $COP_u = 0,262$, nhiệt độ hơi $NH_3$ ra khỏi bình phát sinh $t_5 = 116,5^\circ\text{C}$, nồng độ dung dịch loãng $29,14%$ và dung dịch đặc $34,11%$.

Literature Review và Positioning

Lịch sử phát triển các mô hình truyền nhiệt - truyền chất trong màng chảy hấp thụ ghi nhận nhiều mốc lý thuyết quan trọng. Nakoryakov và Grigor (1977) đặt nền móng bằng mô hình màng chảy tầng trên vách phẳng đứng đẳng nhiệt sử dụng phương pháp khai triển Fourier, bỏ qua biến thiên chiều dày màng. Grossman (1983) nâng cấp mô hình trên vách đứng cho cả hai điều kiện biên đẳng nhiệt và đoạn nhiệt bằng phương pháp sai phân hữu hạn. Andberg và Vliet (1983, 1987) mở rộng bài toán sang ống giải nhiệt nằm ngang và chỉ ra rằng khi dung dịch tiến gần đến trạng thái cân bằng pha, động lực hấp thụ suy giảm đột ngột, trong đó lưu lượng khối lượng dung dịch là biến số chi phối mạnh nhất đến chiều dài yêu cầu của thiết bị. Habib và Wood (1990), Conlisk (1995) cùng Chowdhury và cộng sự (1998) tiếp tục tinh chỉnh mô hình màng chảy ống ngang thông qua lý thuyết lớp biên và các phép chuyển đổi tọa độ để tích hợp sự biến thiên chiều dày màng lỏng dọc theo góc ôm của ống.

Trong y văn tồn tại hai cuộc tranh luận học thuật sâu sắc:

  1. Trở lực truyền pha: Một trường phái cho rằng trở lực pha lỏng (khuếch tán khối lượng của $H_2O$ và $NH_3$) hoàn toàn kiểm soát tốc độ hấp thụ; trong khi trường phái khác (đại diện bởi Meacham & Garimella, 2002; Killion & Garimella, 2001) chứng minh rằng trong điều kiện dòng hơi chuyển động đối lưu cưỡng bức, sự kết hợp giữa hiệu ứng pha hơi và hiện tượng quá lạnh dung dịch cục bộ có thể làm đảo chiều tỷ trọng cản trở nhiệt - chất.
  2. Ảnh hưởng của lưu lượng dung dịch loãng: Meacham & Garimella (2002) và Sangsoo Lee và cộng sự (2002) chứng minh thực nghiệm rằng tăng lưu lượng dung dịch sẽ tăng đáng kể hệ số truyền nhiệt $k$. Ngược lại, Goel và Goswami (2005) khẳng định lưu lượng dung dịch hầu như không làm thay đổi hệ số truyền nhiệt tổng mà chỉ làm tăng công suất nhiệt do bề mặt truyền nhiệt bị giới hạn ở phía nước giải nhiệt.

Luận án của tác giả Nguyễn Hiếu Nghĩa định vị chính xác vào giao điểm của hai cuộc tranh luận này. Nghiên cứu đã chứng minh rằng khi đường kính ống nhỏ ($\Phi 9,6\text{ mm}$) và dòng chảy ở chế độ giọt ($Re = 4\Gamma/\mu < 100$), trở lực phía nước giải nhiệt bên trong ống và độ ướt bề mặt ngoài đóng vai trò quyết định tương đương với trở lực khuếch tán màng lỏng.

Khi đối chiếu với các công trình quốc tế tiêu biểu:

  • So với nghiên cứu của Sangsoo Lee và cộng sự (2002) trên chùm ống ngang $\Phi 9,5\text{ mm}$ sử dụng dung dịch nồng độ $25%$ ở áp suất $2,5\text{ bar}$, công bố giá trị $k = 0,88\text{ kW}/(\text{m}^2\cdot\text{K})$ và $h_m = 1,65\times 10^{-5}\text{ m/s}$, kết quả của luận án tại nồng độ $30,8%$, mật độ phân phối $\Gamma = 0,008\text{ kg}/(\text{m}\cdot\text{s})$, nhiệt độ nước $306,3\text{ K}$ đạt $k = 0,863\text{ kW}/(\text{m}^2\cdot\text{K})$ và $h_m = 1,45\times 10^{-5}\text{ m/s}$. Sự tương đồng chặt chẽ này khẳng định tính đúng đắn của mô hình.
  • So với nghiên cứu của Brice Le Lostec và cộng sự (2008) tại Canada trên chiller $10\text{ kW}$ nồng độ nạp $66,66%$ đạt $COP = 0,6$ cho mục đích điều hòa không khí, nghiên cứu của luận án làm việc ở vùng nhiệt độ bay hơi âm sâu phục vụ sản xuất đá ($-10^\circ\text{C}$ đến $-19^\circ\text{C}$) đạt $COP = 0,436$, minh chứng sự phù hợp nhiệt động trong điều kiện vận hành khắc nghiệt hơn nhiều.

Đó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 lý thuyết màng chảy tầng rơi tự do của Nusselt và lý thuyết hai màng (Two-film theory) của Whitman sang bài toán truyền nhiệt - truyền chất phi dừng, phi đẳng nhiệt của hỗn hợp hai cấu tử biến tính liên tục trên bề mặt cong:

  • Mở rộng lý thuyết lớp biên nhiệt - nồng độ: Thiết lập quy luật biến thiên không gian hai chiều của trường vận tốc $u(x, y)$, nhiệt độ $T(x, y)$ và nồng độ $C(x, y)$ dọc theo tọa độ cung $\varepsilon(x)$ của nửa chu vi ống tròn.
  • Mệnh đề lý thuyết 1: Tại $1/4$ cung đầu tiên tính từ đỉnh ống ($\varepsilon = 0$ đến $\pi/4$), hệ số tỏa nhiệt của màng $\alpha_{iw}$ và hệ số truyền nhiệt tổng $k$ tăng dần đạt đỉnh do sự xáo trộn màng lỏng khi va chạm rơi tự do, sau đó suy giảm liên tục ở $3/4$ cung còn lại do sự phát triển dày lên của lớp biên nhiệt và sự suy giảm chênh lệch nồng độ động lực.
  • Mệnh đề lý thuyết 2: Hệ số truyền chất $h_m$ có sự tăng vọt cục bộ (spike) ngay tại điểm tiếp xúc lỏng - hơi ở đỉnh ống, sau đó rơi dốc đứng và duy trì dạng đường nằm ngang tiệm cận trạng thái bão hòa trước khi rời khỏi đáy ống.

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 hệ phương trình đạo hàm riêng phân loại parabolic mô tả ba định luật bảo toàn:

  1. Phương trình bảo toàn động lượng màng lỏng: $$\mu \frac{\partial^2 u}{\partial y^2} + \rho g \sin\theta = 0$$
  2. Phương trình bảo toàn năng lượng có đối lưu và khuếch tán: $$u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} = a_q \frac{\partial^2 T}{\partial y^2}$$
  3. Phương trình bảo toàn nồng độ cấu tử $NH_3$: $$u \frac{\partial C}{\partial x} + v \frac{\partial C}{\partial y} = D \frac{\partial^2 C}{\partial y^2}$$

Điều kiện biên tiếp xúc lỏng - hơi ($y = \delta$) được xác lập dựa trên cân bằng nhiệt - chất liên tục:

  • Cân bằng thông lượng nhiệt hấp thụ: $$q_f = -\lambda \left.\frac{\partial T}{\partial y}\right|{y=\delta} = j{NH_3} \cdot i_{ab} - \alpha_v (T_{int} - T_v)$$
  • Cân bằng dòng khối lượng hấp thụ: $$j_{NH_3} = -\rho D \left.\frac{\partial C}{\partial y}\right|{y=\delta} = \rho \cdot h_m (C{int} - C_{al})$$

Điều kiện biên vách ống ($y = 0$): Dẫn nhiệt liên tục qua vách đồng sang dòng nước giải nhiệt chảy cưỡng bức có hệ số đối lưu $\alpha_w$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng (Positivism) kết hợp mô hình hóa giải tích - số trị (Empirical-Analytical Paradigm). Quy trình tích hợp đa phương pháp bao gồm:

  1. Mô phỏng động học nhiệt động toàn hệ thống trên môi trường máy tính nhằm thu hẹp miền khảo sát thông số.
  2. Giải hệ phương trình vi phân đạo hàm riêng (PDE) phi tuyến mô tả quá trình truyền nhiệt - truyền chất 2D cho phần tử thể tích màng ống bằng phương pháp sai phân hữu hạn (Finite Difference Method).
  3. Thí nghiệm thực nghiệm kiểm chứng trên mô hình máy lạnh hấp thụ một cấp tải nhiệt thực $2\text{ kW}$.

Quy trình nghiên cứu rigorous

Hệ thống thực nghiệm được thiết kế và chế tạo với độ chính xác cao:

  • Cấu trúc bình hấp thụ: Vỏ hộp chữ nhật chứa chùm ống song song gồm 6 cột, 28 hàng ống đồng trơn $\Phi 9,6\text{ mm}$, chiều dài hiệu dụng $292\text{ mm}$, bước ống đứng $S_v = 20\text{ mm}$, bước ống ngang $S_h = 13\text{ mm}$.
  • Xử lý bề mặt và phân phối: Toàn bộ bề mặt ngoài ống được thổi cát kỹ thuật để tăng năng lượng bề mặt và độ dẫn ướt ($WR$). Ống phân phối dung dịch đỉnh khoan hàng lỗ $\Phi 1,2\text{ mm}$ với bước $4\text{ mm}$ đảm bảo chế độ chảy giọt ($Re < 100$) đều khắp các hàng ống.
  • Hệ thống đo lường:
    • Cảm biến nhiệt độ Pt100 độ chính xác $\pm 0,05^\circ\text{C}$ tại tất cả các điểm nút trạng thái (nhiệt độ phát sinh $t_g$, ngưng tụ $t_c$, bay hơi $t_e$, hấp thụ $t_a$, hơi ra $t_5$, nước muối vào/ra).
    • Đồng hồ đo lưu lượng thể tích dạng phao và turbine độ chính xác cao cho dòng nước giải nhiệt ($10\text{--}20\text{ l/phút}$) và dung dịch tuần hoàn.
    • Lấy mẫu hóa nghiệm nồng độ dung dịch loãng và đặc bằng phương pháp chuẩn độ axit - bazơ kết hợp đo tỷ trọng kế với độ không đảm bảo đo $< 0,3%$.

Data và phân tích

Thuật toán số trị được lập trình trên ngôn ngữ MATLAB để giải lặp hệ phương trình màng lỏng 2D:

  • Miền tính toán được rời rạc hóa với lưới phi đồng nhất tập trung mật độ cao tại bề mặt tự do $y = \delta$.
  • Kiểm tra tính ổn định và hội tụ của thuật toán số qua tiêu chuẩn sai số tương đối nhỏ hơn $10^{-6}$.
  • Đánh giá sai số hệ thống giữa mô phỏng và thực nghiệm:
    • Sai số công suất nhiệt bình hấp thụ giữa thuật toán tính hệ thống ($Q_{a_compute}$) và đo đạc thực nghiệm ($Q_{a_meas}$) là $4,3%$.
    • Sai số giữa mô phỏng số chi tiết từ hàm quan hệ ($Q_{a_sim}$) so với thực nghiệm ($Q_{a_meas}$) là $12,3%$.
    • Sai số trung bình giữa hệ số $COP$ lý thuyết và $COP$ thực nghiệm trong dải làm nước đá ($-10^\circ\text{C}$ đến $-19^\circ\text{C}$) là $1,2%$.

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

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

  1. Xác lập nồng độ dung dịch nạp tối ưu cho mục đích sản xuất đá: Khảo sát chuỗi nồng độ nạp từ $29,5%$ đến $32,5%$ đã chứng minh $31%$ là điểm tối ưu tuyệt đối. Tại điểm này, hệ thống đạt $COP = 0,436$, hiệu suất làm lạnh nước muối $COP_u = 0,262$, nhiệt độ hơi $NH_3$ rời bình phát sinh đạt $t_5 = 116,5^\circ\text{C}$, nồng độ dung dịch loãng $29,14%$ và nồng độ dung dịch đặc $34,11%$.
  2. Quy luật phân bố cục bộ của hệ số truyền nhiệt và truyền chất: Trên nửa chu vi ống giải nhiệt, hệ số tỏa nhiệt màng $\alpha_{iw}$ và hệ số truyền nhiệt tổng $k$ tăng ở $1/4$ góc ôm đầu tiên rồi giảm dần về phía đáy ống. Hệ số truyền chất $h_m$ đạt giá trị đỉnh tại vị trí tiếp xúc tức thời rồi suy giảm nhanh chóng do hiệu ứng bão hòa bề mặt.
  3. Xây dựng tương quan thực nghiệm - số trị đa biến: Thiết lập thành công hai phương trình hồi quy xác định hệ số truyền nhiệt $k$ $[\text{W}/(\text{m}^2\cdot\text{K})]$ và hệ số truyền chất $h_m$ $[\text{m/s}]$ dưới dạng hàm của nồng độ dung dịch $C_{al}$ ($28%\text{--}31%$), mật độ phân phối dòng $\Gamma$ ($0,001\text{--}0,03\text{ kg}/(\text{m}\cdot\text{s})$) và nhiệt độ nước giải nhiệt $T$ ($301\text{--}311\text{ K}$): $$k = f(C_{al}, \Gamma, T) = f(0,308;, 0,008;, 306,3) = 0,863\text{ kW}/(\text{m}^2\cdot\text{K})$$ $$h_m = f(C_{al}, \Gamma, T) = f(0,308;, 0,008;, 306,3) = 1,45\times 10^{-5}\text{ m/s}$$
  4. Tương quan nhiệt độ phát sinh tối ưu: Xác lập phương trình đa thức biểu diễn quan hệ giữa nhiệt độ phát sinh tối ưu $t_{g_opt}$ theo nhiệt độ bay hơi $t_e$, nhiệt độ ngưng tụ $t_c$ và nhiệt độ dung dịch ra khỏi bình hấp thụ $t_a$, cho phép tự động hóa điều khiển chu trình khi thông số thời tiết biến động.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp dữ liệu thực nghiệm chuẩn xác khẳng định mức độ chi phối của trở lực nhiệt phía nước giải nhiệt bên trong ống đối với tổng trở lực truyền chất màng bên ngoài, hoàn thiện lý thuyết truyền nhiệt - truyền chất đồng thời trong chất lưu hai cấu tử.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình thiết kế tích hợp từ phần tử thể tích 2D PDE đến mô phỏng chu trình hệ thống và kiểm chứng thực nghiệm quy mô công nghiệp nhỏ ($2\text{ kW}$).
  • Về mặt kỹ thuật ứng dụng: Chứng minh rằng bình hấp thụ ống đồng trơn kết hợp công nghệ thổi cát và cơ cấu phân phối nhỏ giọt hoàn toàn có thể chế tạo tại các xưởng cơ khí trong nước mà không cần đầu tư dây chuyền nhập ngoại đắt tiền.
  • Về chính sách và môi trường: Cung cấp luận cứ kỹ thuật để triển khai các hệ thống máy đá hấp thụ sử dụng nguồn nhiệt thải công nghiệp, nhiệt lò hơi sinh khối (trấu, mùn cưa) hoặc năng lượng mặt trời tại vùng sâu, vùng xa và hải đảo Việt Nam.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, nghiên cứu vẫn tồn tại một số giới hạn học thuật cần ghi nhận khách quan:

  1. Giả thiết mô hình hóa hình học 2D: Mô hình số đơn giản hóa phần tử thể tích ống thành bài toán 2D đối xứng dọc theo nửa chu vi, chưa tính toán hết hiện tượng xáo trộn 3D và sự hình thành giọt rơi không ổn định tại các khoảng hở giữa các hàng ống trong không gian ba chiều.
  2. Bỏ qua hiệu ứng biến dạng động của sóng màng: Mô hình coi màng chảy là bề mặt phẳng trơn ổn định, trong khi trên thực tế sự gợn sóng của màng lỏng (wave roll) có thể làm gia tăng đối lưu tự nhiên và hệ số truyền chất cục bộ.
  3. Phạm vi quy mô công suất: Các kết quả thực nghiệm mới được kiểm chứng trên mô hình công suất $2\text{ kW}$ (quy mô phòng thí nghiệm/bán công nghiệp), đòi hỏi thêm các hiệu chỉnh hệ số khi nhân rộng lên các tổ máy công nghiệp quy mô $20\text{--}100\text{ kW}$.

Chương trình nghiên cứu mở rộng trong tương lai bao gồm 5 định hướng:

  • Ứng dụng hạt nano dẫn nhiệt (nanofluids) vào dung dịch $NH_3\text{-}H_2O$ nhằm cải thiện độ dẫn nhiệt riêng của màng lỏng.
  • Nghiên cứu cơ chế tăng cường truyền chất Marangoni thông qua việc bổ sung các phụ gia hoạt tính bề mặt bậc cao.
  • Mở rộng mô hình số 3D Dynamic CFD mô phỏng chính xác sự va đập của giọt lỏng từ hàng ống trên xuống hàng ống dưới.
  • Tích hợp hệ thống điều khiển thông minh tự động điều chỉnh nhiệt độ phát sinh theo thuật toán bám điểm cực đại ($MPPT$) kết hợp nguồn nhiệt mặt trời.
  • Thử nghiệm trên các cấu trúc ống có cánh vi mô nội tại (micro-fin tubes) được gia công bằng công nghệ in 3D kim loại.

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

Công trình nghiên cứu tạo ra những tác động khoa học và kinh tế - xã hội rõ nét:

  • Tác động học thuật: Đóng góp cơ sở dữ liệu thực nghiệm chuẩn xác về hệ số truyền nhiệt và truyền chất của cặp $NH_3\text{-}H_2O$ trên chùm ống ngang, là tài liệu tham khảo cốt lõi cho các nghiên cứu tiếp theo về kỹ thuật nhiệt lạnh và năng lượng tái tạo tại Việt Nam và khu vực Đông Nam Á.
  • Chuyển dịch công nghiệp: Mở ra giải pháp khả thi để các doanh nghiệp cơ điện lạnh trong nước tự chủ thiết kế, gia công các tổ máy sản xuất nước đá cây, đá vảy phục vụ tàu đánh bắt cá xa bờ sử dụng nguồn nhiệt thải từ khí xả động cơ diesel chính, giúp tiết kiệm $100%$ điện năng chạy máy nén lạnh truyền thống.
  • Lợi ích xã hội và môi trường: Sử dụng $NH_3$ (R717) – môi chất tự nhiên có chỉ số phá hủy tầng ozone $ODP = 0$ và tiềm năng gây hiệu ứng nhà kính $GWP = 0$, góp phần thực hiện Nghị định thư Montreal và Thỏa thuận Paris về giảm phát thải khí nhà kính.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận giải số PDE màng chảy 2D và quy trình thiết kế mô hình thực nghiệm chuẩn mực.
  • Các nhà khoa học nhiệt - lạnh: Sử dụng bộ tương quan thực nghiệm $k = f(C, \Gamma, T)$ và $h_m = f(C, \Gamma, T)$ làm chuẩn đối sánh cho các mô hình lý thuyết mới.
  • Kỹ sư R&D công nghiệp: Ứng dụng trực tiếp bản vẽ cấu tạo, kích thước ống $\Phi 9,6\text{ mm}$, khoảng cách bước ống $20\times 13\text{ mm}$ và thông số bộ phân phối $\Phi 1,2\text{ mm}$ để chế tạo thiết bị hấp thụ thương mại.
  • Nhà hoạch định chính sách năng lượng: Có căn cứ khoa học vững chắc để xây dựng các chương trình tài trợ dự án cấp nhiệt lạnh phi lưới điện cho ngư dân và đồng bào hải đảo.

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 mở rộng Lý thuyết màng chảy tầng Nusselt và Lý thuyết truyền chất hai màng sang bài toán truyền nhiệt - truyền chất đồng thời phi đẳng nhiệt của dung dịch hai cấu tử $NH_3\text{-}H_2O$ trên chùm ống tròn nằm ngang trong điều kiện đối lưu ngược chiều của dòng hơi. Luận án đã mô hình hóa và kiểm chứng thành công quy luật biến thiên không đồng nhất của hệ số truyền nhiệt $k$ và truyền chất $h_m$ dọc theo góc ôm nửa chu vi ống, giải quyết điểm nghẽn lý thuyết về sự suy giảm động lực hấp thụ khi lớp biên màng lỏng dày lên và tiến tới trạng thái cân bằng pha.

2. Điểm cải tiến phương pháp luận của nghiên cứu so với các công trình trước đây là gì?

So với nghiên cứu của Grossman (1983) chỉ giới hạn ở vách phẳng hoặc Chowdhury và cộng sự (1998) chỉ mô phỏng ống đơn cô lập với nhiệt độ vách không đổi, luận án đã: (1) Thiết lập mô hình phần tử thể tích liên kết trên chùm ống 28 hàng có xét đến sự tăng dần nhiệt độ của dòng nước giải nhiệt ngược chiều; (2) Tích hợp mô hình số 2D PDE vào chương trình mô phỏng chu trình toàn hệ thống; (3) Kiểm chứng đồng thời từ cấp độ phần tử đến cấp độ hệ thống thực $2\text{ kW}$ với sai số công suất nhiệt chỉ $4,3%$.

3. Phát hiện nào trong dữ liệu thực nghiệm mang tính bất ngờ và có ý nghĩa nhất?

Phát hiện về nồng độ nạp tối ưu $31%$ mang tính bước ngoặt. Khác với các giả định truyền thống cho rằng tăng nồng độ nạp sẽ luôn làm tăng năng suất lạnh, thực nghiệm chỉ ra nồng độ vượt quá $31%$ làm tăng đột biến áp suất bay hơi và nhiệt độ sôi dung dịch, làm giảm chênh lệch nhiệt độ truyền nhiệt trong bình phát sinh và giảm $COP_u$. Ngược lại, nồng độ dưới $31%$ làm nghèo dòng môi chất sinh ra, giảm mạnh sản lượng đá. Mức $31%$ là điểm cân bằng hoàn hảo giữa hiệu suất nhiệt động và công suất hữu ích ($COP = 0,436; COP_u = 0,262$).

4. Nghiên cứu có cung cấp quy trình nhân bản (replication protocol) đầy đủ không?

Có. Luận án cung cấp chi tiết toàn bộ bản vẽ chế tạo cơ khí bình hấp thụ, kích thước hình học chùm ống ($\Phi 9,6\text{ mm}$, bước dọc $20\text{ mm}$, bước ngang $13\text{ mm}$), thông số bộ phân phối ($\Phi 1,2\text{ mm}$, bước $4\text{ mm}$), mã nguồn chương trình mô phỏng trên MATLAB (Phụ lục 2.1, 2.2, 2.3) và toàn bộ nhật ký số liệu đo đạc thực nghiệm các điểm nút trạng thái.

5. Lộ trình nghiên cứu phát triển 10 năm tiếp theo được định hình ra sao?

Lộ trình 10 năm tập trung vào 3 giai đoạn: (1) Tối ưu hóa bề mặt trao đổi nhiệt bằng vật liệu mới và phụ gia tăng cường truyền chất Marangoni; (2) Module hóa hệ thống lạnh hấp thụ công suất $10\text{--}50\text{ kW}$ tích hợp bộ thu nhiệt mặt trời ống chân không và tận dụng nhiệt khí thải tàu cá; (3) Thương mại hóa máy làm đá hấp thụ tự động hóa hoàn toàn phục vụ chiến lược kinh tế biển bền vững.

Kết luận

  1. Nghiên cứu đã giải quyết thành công bài toán cốt lõi về truyền nhiệt - truyền chất đồng thời trong bình hấp thụ kiểu màng chảy trên chùm ống ngang của máy lạnh hấp thụ $NH_3\text{-}H_2O$ liên tục, hoàn toàn làm chủ thiết kế và công nghệ chế tạo trong điều kiện Việt Nam.
  2. Xây dựng thành công mô hình toán - số 2D PDE trên MATLAB giải các phương trình bảo toàn động lượng, năng lượng và nồng độ, giải thích rõ quy luật phân bố cục bộ của trường nhiệt độ, nồng độ và các hệ số $k$, $h_m$ trên bề mặt ống.
  3. Thiết lập chính xác tương quan thực nghiệm - số trị xác định hệ số truyền nhiệt $k = 0,863\text{ kW}/(\text{m}^2\cdot\text{K})$ và hệ số truyền chất $h_m = 1,45\times 10^{-5}\text{ m/s}$, đạt độ tương đồng cao với các công bố quốc tế uy tín (Sangsoo Lee và cộng sự, 2002).
  4. Chế tạo và vận hành thành công hệ thống thực nghiệm máy lạnh hấp thụ $2\text{ kW}$ sản xuất nước đá (nhiệt độ nước muối $-10^\circ\text{C}$ đến $-19^\circ\text{C}$), xác định nồng độ nạp tối ưu $31%$, đạt $COP = 0,436$ và $COP_u = 0,262$ với sai số mô phỏng - thực nghiệm chỉ $1,2%$.
  5. Mở ra hướng đi mang tính ứng dụng thực tiễn cao trong việc sử dụng năng lượng nhiệt thải và năng lượng tái tạo phục vụ chuỗi cung ứng lạnh và chế biến thủy hải sản, đóng góp thiết thực vào mục tiêu phát triển kinh tế biển và bảo vệ môi trường bền vững.