Bối cảnh và vấn đề nghiên cứu

Trong bối cảnh điều kiện sống ngày càng nâng cao, nhu cầu sử dụng các hệ thống điều hòa không khí và làm lạnh gia tăng nhanh chóng. Phần lớn các hệ thống điều hòa không khí truyền thống—chẳng hạn như hệ thống làm lạnh cục bộ, hệ thống làm lạnh trung tâm sử dụng máy làm lạnh nước (water chiller) hay hệ thống điều hòa biến lưu lượng môi chất (VRV)—hoạt động dựa trên chu trình nén hơi dẫn động bằng điện năng. Việc tiêu thụ lượng lớn điện năng sản xuất từ nhiên liệu hóa thạch không chỉ làm cạn kiệt tài nguyên năng lượng mà còn gây ô nhiễm môi trường và gia tăng phát thải khí nhà kính.

Trước thách thức này, công nghệ làm lạnh và điều hòa không khí sử dụng thiết bị ejector nhận được sự quan tâm nghiên cứu lớn. Ejector hoạt động như một máy nén nhiệt (thermal compressor), vận hành bằng nguồn nhiệt đầu vào như nhiệt bức xạ mặt trời hoặc nhiệt thải công nghiệp thay vì tiêu tốn cơ năng. Chu trình làm lạnh ejector sở hữu nhiều ưu điểm cơ học: kết cấu đơn giản, không có các chi tiết chuyển động cơ học, chi phí vận hành và bảo dưỡng thấp, kích thước nhỏ gọn. Ngoài ứng dụng trong chu trình lạnh, ejector còn được sử dụng rộng rãi trong các lĩnh vực công nghiệp khác như tăng áp chất lỏng, hút tách khí chân không, khử mùi dầu ăn, chế biến hóa dầu và chưng cất dầu mỏ.

Hiệu quả vận hành của chu trình làm lạnh ejector phụ thuộc chủ yếu vào hệ số cuốn (entrainment ratio). Tuy nhiên, dòng môi chất bên trong ejector là dòng chuyển động có tính nén được, đạt vận tốc cao (từ dưới âm đến siêu âm), xuất hiện các hiện tượng thủy khí động lực học phức tạp như sóng xung kích, lớp biên và tương tác hòa trộn rối. Việc tiến hành khảo sát thực nghiệm đối với các trường dòng phức tạp này thường đòi hỏi chi phí lớn và thời gian lắp đặt kéo dài. Do đó, phương pháp mô phỏng số động lực học chất lưu (Computational Fluid Dynamics - CFD) là phương án tối ưu để khảo sát chi tiết trường áp suất, nhiệt độ, vận tốc và hiệu suất của thiết bị.

Đề tài tốt nghiệp "Nghiên cứu số đặc tính dòng chảy và nhiệt bên trong thiết bị ejector ở các điều kiện vận hành khác nhau" (Numerical investigation flow and thermal characteristics inside of the ejector device at different operating conditions) được nhóm sinh viên Trần Đình Trung Đức, Đặng Minh Đức và Võ Tiến Dũng thực hiện dưới sự hướng dẫn của TS. Trần Thanh Tịnh tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh.

Mục tiêu và nhiệm vụ nghiên cứu của đồ án được xác định cụ thể trong văn bản:

  1. Mục tiêu nghiên cứu: Mô phỏng thiết bị ejector với các thông số đầu vào khác nhau, từ đó đánh giá hiệu suất vận hành (thông qua hệ số cuốn) và so sánh, đối chiếu với cơ sở lý thuyết khí động học của thiết bị.
  2. Nhiệm vụ nghiên cứu:
    • Thực hiện mô phỏng số trường dòng chảy và đặc tính nhiệt bên trong ejector ở các chế độ vận hành khác nhau.
    • Trích xuất dữ liệu mô phỏng để xây dựng đồ thị và xác định hệ số cuốn.
    • Phân tích dòng chảy ở các chế độ, đánh giá kết quả từ đồ thị hệ số cuốn và so sánh với kết quả tính toán lý thuyết cũng như các nghiên cứu công bố liên quan.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Thiết bị ejector và cụm vòi phun (nozzle) có biên dạng dòng cải tiến theo đường cong Witozinsky.
  • Phạm vi nghiên cứu:
    • Môi chất làm việc: Khảo sát với môi chất làm lạnh R134a ở trạng thái hơi bão hòa tại nhiệt độ $15^\circ\text{C}$.
    • Điều kiện vận hành: Khảo sát dòng môi chất nén được với lưu lượng khối lượng dòng sơ cấp (motive flow) thay đổi (lên đến $0{,}289219\text{ kg/s}$) trong khi cố định lưu lượng dòng thứ cấp (entrained flow); khảo sát các trường hợp số Mach tại vòi phun ($M = 1$, $M < 1$, $M > 1$); khảo sát sự biến thiên của áp suất nén đầu ra ($P_c$) và phân tích chi tiết 6 trường hợp vận hành (Case 1 đến Case 6) cùng 3 mặt cắt khảo sát đặc trưng (B1, B2, B3 hoặc P1, P2, P3).
    • Thời gian thực hiện: Từ tháng 02/2023 đến tháng 07/2023.

Cơ sở lý thuyết và phương pháp

Đồ án xây dựng nền tảng phân tích dựa trên sự kết hợp giữa lý thuyết khí động học dòng nén được một chiều, các phương trình bảo toàn cơ học chất lưu vi phân và phương pháp thể tích hữu hạn (Finite Volume Method - FVM).

Khung lý thuyết và mô hình toán học bao gồm:

  1. Quan hệ diện tích - vận tốc cho dòng khí nén được: Mối liên hệ vi phân giữa biến thiên diện tích mặt cắt ngang $dA$ và biến thiên vận tốc $du$ dọc theo dòng chảy được mô tả qua phương trình: $$\frac{dA}{A} = (M^2 - 1) \frac{du}{u}$$ Trong đó $A$ là diện tích mặt cắt ngang, $u$ là vận tốc trung bình, $M$ là số Mach ($M = u/a$). Phương trình chỉ ra quy luật hình học:
    • Với dòng dưới âm ($0 \le M < 1$): Để tăng vận tốc ($du > 0$), diện tích ống phải thu hẹp ($dA < 0$, ống hội tụ).
    • Với dòng siêu âm ($M > 1$): Để tăng vận tốc ($du > 0$), diện tích ống phải mở rộng ($dA > 0$, ống phân kỳ).
    • Tại tiết diện tới hạn hay cổ vòi phun ($M = 1$): $dA = 0$, tương ứng với vị trí có diện tích mặt cắt nhỏ nhất (throat).
  2. Hệ phương trình điều khiển dòng chất lưu (Navier-Stokes):
    • Phương trình liên tục (Bảo toàn khối lượng): $$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \vec{v}) = 0$$
    • Phương trình động lượng (Bảo toàn số lượng chuyển động): $$\frac{\partial (\rho \vec{v})}{\partial t} + \nabla \cdot (\rho \vec{v}\vec{v}) = -\nabla P + \nabla \cdot (\tau) + \rho \vec{g}$$ với $\tau$ là ten-sơ ứng suất nhớt, $\mu_e$ là độ nhớt động học hiệu dụng.
    • Phương trình năng lượng: $$\frac{\partial (\rho E)}{\partial t} + \nabla \cdot (\vec{v}(\rho E + P)) = \nabla \cdot \left(k_e \nabla T - \sum h_j J_j + (\tau \cdot \vec{v})\right)$$ với $E$ là năng lượng toàn phần, $k_e$ là độ dẫn nhiệt hiệu dụng, $h$ là entanpi, $J$ là dòng khuếch tán.
  3. Các thông số hiệu suất và thủy động lực học:
    • Số Reynolds: $Re = \frac{\rho v d_h}{\mu}$
    • Vận tốc trung bình theo lưu lượng khối lượng $G$: $v = \frac{4G}{\pi \rho d_h^2}$
    • Hệ số cuốn của thiết bị ejector ($U$): $$U = \frac{m_e}{m_m}$$ với $m_e$ là lưu lượng khối lượng dòng thứ cấp (entrained flow) và $m_m$ là lưu lượng khối lượng dòng sơ cấp (motive flow).
  4. Mô hình hình học vòi phun cải tiến: Đồ án kế thừa nghiên cứu của Xinyue Hao và cộng sự (2021) về vòi phun cải tiến biên dạng theo đường cong Witozinsky cho đoạn hội tụ và buồng hút. Đường cong này tối ưu hóa sự chuyển tiếp của dòng chảy, làm giảm tổn thất ma sát và tổn thất bất thuận nghịch so với thiết kế vòi phun hình nón phẳng truyền thống.

Phương pháp nghiên cứu thực tế:

  • Phương pháp số: Ứng dụng phương pháp thể tích hữu hạn (FVM) tích hợp trong phần mềm thương mại ANSYS Fluent.
  • Mô hình rối: Sử dụng mô hình nhiễu loạn Realizable $k-\varepsilon$, phù hợp cho các dòng xoáy và dòng tia nén tốc độ cao.
  • Quy trình chia lưới và kiểm tra tính độc lập của lưới: Tiến hành phân chia lưới không gian với 3 cấp độ lưới (M1, M2, M3) bằng công cụ ANSYS Meshing. Đánh giá chất lượng phần tử qua các thông số: tỉ số hướng (Aspect Ratio), độ méo (Skewness), chất lượng trực giao (Orthogonal Quality) và độ trơn tru (Smoothness). Phân tích sự hội tụ của thặng dư sai số (residual) và hồ sơ vận tốc để chọn lưới tính toán tối ưu.

Nguồn dữ liệu:

  • Thông số nhiệt động và tính chất vật lý của môi chất làm lạnh R134a tại trạng thái hơi bão hòa $15^\circ\text{C}$ (khối lượng riêng, độ nhớt, nhiệt dung riêng, độ dẫn nhiệt).
  • Kích thước hình học tiêu chuẩn của thiết bị ejector thử nghiệm dựa theo công bố khoa học thực nghiệm đối chứng của Xinyue Hao và cộng sự (2021).

Nội dung chính theo từng chương

Chương 1: Giới thiệu tổng quan (Introduction)

Chương mở đầu phân tích bối cảnh phát triển của kỹ thuật điều hòa không khí và làm lạnh, nêu bật các hạn chế về mặt tiêu thụ năng lượng và khí thải của các hệ thống nén hơi cơ học dẫn động điện truyền thống. Tác giả trình bày lịch sử của ejector hơi từ chu trình Leblanc (1910) và chỉ ra tiềm năng của ejector khi tận dụng nguồn nhiệt thải hoặc năng lượng mặt trời. Chương này lý giải nguyên nhân chọn phương pháp mô phỏng số CFD nhằm khắc phục các nhược điểm tốn kém chi phí, thời gian của nghiên cứu thực nghiệm. Đồng thời, chương 1 xác lập mục tiêu, 3 nhiệm vụ cụ thể, đối tượng nghiên cứu là thiết bị vòi phun ejector và phạm vi khảo sát điều kiện biên lưu lượng dòng sơ cấp lên đến $0{,}289219\text{ kg/s}$ với môi chất R134a.

Chương 2: Cơ sở lý thuyết (Theoretical Background)

Chương 2 cung cấp cơ sở lý thuyết toàn diện về ejector hơi:

  • Cấu tạo tiêu chuẩn gồm 4 bộ phận chính: vòi phun sơ cấp (nozzle), buồng hút (suction chamber), buồng hòa trộn (mixing chamber) và ống khuếch tán (diffuser).
  • Phân tích nguyên lý làm việc dựa trên nguyên lý Bernoulli và hiệu ứng Venturi: dòng sơ cấp áp suất cao giãn nở qua vòi phun để chuyển hóa áp năng thành động năng siêu âm, tạo ra vùng áp suất chân không tại buồng hút để cuốn dòng thứ cấp; hai dòng chất lưu hòa trộn và tiếp tục hồi phục áp suất khi giảm tốc qua ống khuếch tán.
  • Thiết lập hệ phương trình vi phân chi phối dòng nén được, mối quan hệ diện tích - vận tốc theo số Mach ($dA/A = (M^2 - 1)du/u$), phân tích sự khác biệt giữa dòng dưới âm ($M < 1$), tới hạn ($M = 1$) và siêu âm ($M > 1$).
  • Trình bày hệ phương trình Navier-Stokes (bảo toàn khối lượng, động lượng, năng lượng), công thức xác định số Reynolds ($Re$) và công thức hệ số cuốn ($U = m_e / m_m$).
  • Phân loại ejector theo pha (đơn pha, hai pha), theo đặc tính dòng (nén được, không nén được, dưới âm, siêu âm) và theo kết cấu hình học (tiết diện tròn, chữ nhật; vòi phun đơn, đa vòi phun; hòa trộn đẳng áp hoặc hòa trộn đẳng tích).
  • Phân tích ưu - nhược điểm của ejector và tổng thuật các nghiên cứu tiền bối: Chunnanond & Aphornratana (tổng quan ứng dụng lạnh), Pradeep Gupt et al. (mô hình mạng nơ-ron nhân tạo), Yulei Huang et al. (mô hình vòi phun chuẩn 2D) và trọng tâm là nghiên cứu vòi phun biên dạng Witozinsky của Xinyue Hao et al.

Chương 3: Tổng quan phần mềm ANSYS và mô phỏng CFD (ANSYS Software and CFD Simulation Overview)

Chương 3 trình bày tổng quan về kỹ thuật mô phỏng CFD và các công cụ ứng dụng:

  • Định nghĩa CFD và liệt kê các lĩnh vực công nghiệp ứng dụng (hàng không, ô tô, HVAC/R, hóa chất, thiết bị công nghiệp).
  • Chu trình 8 bước thực hiện mô phỏng CFD chuẩn: (1) Mô tả bài toán; (2) Xây dựng hình học; (3) Tạo lưới; (4) Thiết lập thông số mô hình; (5) Chọn phương pháp giải; (6) Chạy mô phỏng; (7) Kiểm tra tính độc lập của lưới; (8) Khảo sát các trường hợp, xử lý kết quả và xuất báo cáo.
  • Giới thiệu ba phương pháp số rời rạc hóa hệ phương trình: Phương pháp phần tử hữu hạn (FEM), Phương pháp thể tích hữu hạn (FVM) và Phương pháp sai phân hữu hạn (FDM), nhấn mạnh ưu thế bảo toàn thông lượng cục bộ của FVM.
  • Cấu trúc nền tảng ANSYS Workbench kết nối 4 mô-đun: DesignModeler (dựng hình 2D/3D), Meshing (chia lưới), Fluent (cài đặt điều kiện biên và giải số) và CFD-Post (hậu xử lý kết quả).
  • Các tiêu chí kỹ thuật kiểm soát chất lượng lưới: Aspect Ratio (tỉ số giữa cạnh dài nhất và ngắn nhất của phần tử, giá trị lý tưởng bằng 1), Skewness, Orthogonal Quality và Smoothness; phân biệt lưới có cấu trúc (structured grid) và không cấu trúc (unstructured grid).

Chương 4: Xây dựng mô hình và mô phỏng (Model Building and Simulations)

Chương 4 mô tả chi tiết quy trình thiết lập mô hình tính toán:

  • Sử dụng phần mềm SolidWorks và DesignModeler để xây dựng mô hình hình học 3D của thiết bị ejector thử nghiệm dựa trên các kích thước hình học chuẩn (Bảng 4.1), áp dụng hàm giải tích đường cong Witozinsky cho biên dạng đoạn hội tụ vòi phun và buồng hút.
  • Tiến hành chia lưới trong ANSYS Meshing, tạo lập 3 mô hình lưới có mật độ khác nhau (M1, M2, M3) để khảo sát tính độc lập của lưới. Đánh giá chất lượng lưới qua Bảng 4.2 và kiểm tra mức độ hội tụ của các giá trị thặng dư (residual) cùng sự ổn định của hồ sơ vận tốc.
  • Thiết lập mô hình vật lý trong Fluent: kích hoạt phương trình năng lượng (Energy Equation), lựa chọn mô hình dòng rối Realizable $k-\varepsilon$.
  • Khai báo thuộc tính môi chất R134a ở trạng thái hơi bão hòa tại $15^\circ\text{C}$ trong Cell Zone Conditions (Bảng 4.3).
  • Cài đặt điều kiện biên: đầu vào dòng sơ cấp và thứ cấp theo lưu lượng khối lượng (Mass flow inlet: $m_m$ và $m_e$), đầu ra theo áp suất (Pressure outlet: $P_c$), điều kiện thành rắn không trượt (No-slip wall).

Chương 5: Kết quả và thảo luận (Results and Discussion)

Chương 5 là nội dung trọng tâm trình bày toàn bộ số liệu và phân tích trường dòng:

  • Phân tích trường vận tốc, áp suất tĩnh và nhiệt độ của môi chất R134a dọc theo trục tọa độ x và tại khu vực cổ vòi phun.
  • Khảo sát đường dòng (streamlines) và trường vectơ vận tốc, thể hiện rõ quá trình cuốn hút dòng thứ cấp và hiện tượng cuộn xoáy trong buồng hòa trộn.
  • Phân tích chi tiết đặc tính dòng tại 3 mặt cắt ngang khảo sát (B1, B2, B3 ứng với các mặt cắt P1, P2, P3), vẽ đồ thị phân bố vận tốc theo phương hướng kính Y.
  • Trình bày kết quả mô phỏng cho 6 trường hợp vận hành cụ thể (Case 1 đến Case 6), so sánh đối chiếu sâu giữa Case 4 và Case 6 về các đường đẳng trị vận tốc (velocity contour), áp suất tĩnh và nhiệt độ.
  • Khảo sát ảnh hưởng của áp suất nén hỗn hợp tại cửa ra ($P_c$) và lưu lượng dòng sơ cấp lên hệ số cuốn ($U$), biểu diễn qua đồ thị biến thiên $U - P_c$.
  • Phân tích dòng đẳng entropy qua vòi phun dưới âm và siêu âm, đối chiếu quan hệ diện tích - số Mach giữa mô phỏng và lý thuyết khí động học giải tích.

Chương 6: Kết luận và khuyến nghị (Conclusion and Recommendation)

Chương cuối cùng tổng kết các kết quả đạt được từ đồ án: khẳng định tính khả thi của mô hình mô phỏng số trong việc mô tả hiện tượng động lực học và nhiệt bên trong ejector, xác định mối liên hệ giữa các thông số vận hành với hệ số cuốn $U$. Đồng thời, chương này đưa ra các đề xuất về việc tiếp tục tối ưu hóa hình học vòi phun (như vị trí thoát vòi phun - Nozzle Exit Position NXP = 9 mm) và mở rộng phạm vi nghiên cứu thực nghiệm.

Kết quả và đóng góp

Đồ án đã thu thập và xử lý hệ thống số liệu mô phỏng CFD hoàn chỉnh về phân bố nhiệt - thủy động lực học của môi chất R134a trong thiết bị ejector cải tiến.

Dưới đây là bảng tổng hợp các trường hợp mô phỏng và thông số vận hành chính được khảo sát trong đồ án:

Trường hợp Lưu lượng dòng sơ cấp $m_m$ (kg/s) Lưu lượng dòng thứ cấp $m_e$ (kg/s) Áp suất đầu ra $P_c$ (Pa) Khảo sát số Mach ($M$) Vị trí vòi phun NXP (mm)
Case 1 Khảo sát biến thiên Cố định theo điều kiện biên Biến thiên theo tải $M < 1$ (Dưới âm) 9
Case 2 Khảo sát chuẩn Cố định $P_c$ định mức $M < 1$ (Dưới âm) 9
Case 3 Khảo sát tăng áp Cố định Tăng áp suất lưng $M = 1$ (Tới hạn) 9
Case 4 $0{,}289219$ Cố định ($0$ hoặc $0{,}044896$) Áp suất nén tới hạn $M \ge 1$ (Siêu âm) 9
Case 5 Khảo sát đối chứng Cố định Thay đổi áp suất $P_c$ Khảo sát mở rộng 9
Case 6 $0{,}089219$ (Đối chứng bài báo [1]) $0{,}044896$ Biến thiên kiểm tra $U$ Chế độ dòng nén 9

Các kết quả phân tích trường dòng định lượng chính:

  • Phân bố vận tốc và sóng giãn nở: Tại cửa ra vòi phun sơ cấp, vận tốc dòng môi chất R134a tăng vọt khi tiết diện thu hẹp về phía cổ vòi phun (throat), đạt vận tốc âm thanh ($M = 1$) tại cổ và tiếp tục giãn nở đạt vận tốc siêu âm ($M > 1$) trong đoạn phân kỳ trước khi va đập hòa trộn với dòng thứ cấp.
  • Biến thiên áp suất và nhiệt độ: Áp suất tĩnh và nhiệt độ tĩnh giảm mạnh đến giá trị cực tiểu tại khu vực cửa ra vòi phun, tạo độ chân không cần thiết trong buồng hút để hút dòng môi chất thứ cấp. Sau đó, áp suất được hồi phục dần dọc theo buồng hòa trộn và ống khuếch tán.
  • Quy luật biến thiên của hệ số cuốn ($U$): Đồ án xác lập đồ thị đặc tính làm việc của ejector. Khi áp suất đầu ra $P_c$ thấp hơn áp suất tới hạn (critical pressure), ejector làm việc ở chế độ dòng nghẽn kép (double choking) hoặc tới hạn, hệ số cuốn $U$ giữ nguyên giá trị không đổi. Khi $P_c$ vượt quá áp suất tới hạn, dòng thứ cấp bị suy giảm đột ngột và hệ số cuốn $U$ giảm mạnh về 0 (chế độ hỏng/nghẽn dòng).
  • So sánh Case 4 và Case 6: Thể hiện rõ ảnh hưởng của lưu lượng dòng sơ cấp lên khả năng cuốn và khả năng chịu áp suất nén đầu ra $P_c$. Khi lưu lượng $m_m$ tăng lên $0{,}289219\text{ kg/s}$, động năng dòng tia lớn hơn giúp cải thiện áp suất nén cực đại của thiết bị.

Đóng góp khoa học và kỹ thuật tự nêu của đồ án:

  1. Xây dựng thành công quy trình mô phỏng số 3D/đối xứng trục hoàn chỉnh cho thiết bị ejector sử dụng môi chất R134a trên bộ công cụ ANSYS Fluent kết hợp mô hình rối Realizable $k-\varepsilon$.
  2. Đánh giá tính độc lập của lưới một cách hệ thống qua 3 cấp lưới (M1, M2, M3), chứng minh tính hội tụ của giải thuật số thông qua chỉ số thặng dư và sự trùng khớp của hồ sơ vận tốc.
  3. Kiểm chứng sự tương thích giữa lý thuyết khí động học dòng nén một chiều (quan hệ diện tích - số Mach $dA/A = (M^2-1)du/u$) với kết quả số mô phỏng thực tế trên vòi phun có biên dạng cong Witozinsky.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của đồ án:

  • Nghiên cứu giới hạn ở bài toán mô phỏng số đơn pha với giả định môi chất R134a ở trạng thái hơi bão hòa ($15^\circ\text{C}$), chưa tính toán đến hiện tượng chuyển pha ngưng tự cục bộ (non-equilibrium condensation) thường xảy ra khi dòng hơi giãn nở siêu âm giảm áp và giảm nhiệt độ cực nhanh trong vòi phun.
  • Đồ án chưa tiến hành chế tạo mô hình thực nghiệm vật lý trực tiếp để đo đạc thông số thực tế tại phòng thí nghiệm, mà độ tin cậy của mô hình phụ thuộc vào việc so sánh gián tiếp với số liệu công bố trong bài báo thực nghiệm của Xinyue Hao và cộng sự (2021).
  • Phạm vi khảo sát điều kiện biên chủ yếu tập trung vào việc thay đổi lưu lượng dòng sơ cấp và áp suất đầu ra, trong khi lưu lượng dòng thứ cấp được giữ cố định.

Hướng nghiên cứu tiếp được đề xuất:

  • Mở rộng bài toán mô phỏng số sang mô hình dòng hai pha (two-phase flow model) hoặc mô hình hơi thực (real gas model) để xét đến hiện tượng ngưng tụ cục bộ và giọt lỏng trong ejector.
  • Thực hiện thiết kế, gia công chế tạo mô hình ejector thực nghiệm với biên dạng vòi phun Witozinsky để tiến hành đo đạc đối chứng trực tiếp các thông số áp suất, nhiệt độ và hệ số cuốn thực tế.
  • Khảo sát tối ưu hóa đa mục tiêu đối với các thông số hình học: tỉ số diện tích buồng hòa trộn, góc côn ống khuếch tán và vị trí đặt vòi phun (Nozzle Exit Position - NXP) trong các điều kiện tải lạnh biến đổi.
  • Nghiên cứu thử nghiệm bổ sung với các dòng môi chất làm lạnh thế hệ mới có chỉ số GWP thấp nhằm đáp ứng các tiêu chuẩn bảo vệ môi trường khắt khe hơn.

Giá trị tham khảo

Đồ án tốt nghiệp là tài liệu tham khảo có giá trị học thuật và kỹ thuật thực tiễn đối với:

  • Sinh viên và học viên cao học ngành Kỹ thuật Nhiệt, Nhiệt - Lạnh, Năng lượng và Kỹ thuật Cơ khí: Cung cấp tài liệu tham khảo chi tiết về chu trình làm lạnh ejector, các công thức giải tích khí động học dòng nén được và phương pháp luận thiết lập bài toán mô phỏng dòng tia.
  • Kỹ sư và nhà nghiên cứu ứng dụng CFD: Cung cấp quy trình từng bước xây dựng mô hình trên SolidWorks/DesignModeler, kỹ thuật chia lưới và tiêu chí đánh giá chất lượng lưới (Aspect Ratio, Skewness, Orthogonal Quality), thiết lập điều kiện biên trong ANSYS Fluent và phương pháp xử lý dữ liệu trường dòng trong CFD-Post.
  • Người làm đề tài về tối ưu hóa thiết bị trao đổi nhiệt - khối và thu hồi nhiệt thải: Tham khảo phương pháp ứng dụng đường cong Witozinsky vào thiết kế vòi phun nhằm giảm thiểu tổn thất áp suất bất thuận nghịch.

Các phần có giá trị tham khảo cao nhất trong đồ án bao gồm:

  • Chương 2 (Mục 2.3 và 2.8): Phần diễn giải trực quan về mối quan hệ giữa biến thiên diện tích và vận tốc theo số Mach ($dA/A = (M^2-1)du/u$), cùng hệ phương trình Navier-Stokes chi tiết cho dòng nén được.
  • Chương 3 (Mục 3.2): Tiêu chuẩn định lượng chất lượng lưới và quy trình kiểm tra tính độc lập của lưới.
  • Chương 5 (Mục 5.2 và 5.3): Phương pháp trích xuất dữ liệu, vẽ đồ thị phân bố vận tốc dọc trục, phân tích trường dòng tại các mặt cắt ngang và đánh giá đặc tính làm việc của ejector theo hệ số cuốn $U$.

Câu hỏi thường gặp

1. Đồ án sử dụng môi chất làm việc nào trong các mô phỏng CFD và ở trạng thái nhiệt động ra sao?
Môi chất được sử dụng trong toàn bộ các tính toán mô phỏng của đồ án là R134a ở trạng thái hơi bão hòa tại nhiệt độ $15^\circ\text{C}$ (các thông số vật lý được tra cứu và thiết lập chi tiết tại Bảng 4.3 của đồ án). Mặc dù nước được nhắc đến trong chương tổng quan như một môi chất phổ biến của ejector hơi truyền thống, mô phỏng chuyên sâu của đồ án áp dụng cho môi chất lạnh R134a.

2. Vòi phun cải tiến trong đồ án được thiết kế theo biên dạng hình học nào và nhằm mục đích gì?
Vòi phun cải tiến được thiết kế với đoạn hội tụ và buồng hút tuân theo biên dạng đường cong Witozinsky (Witozinsky curve), kế thừa từ nghiên cứu của Xinyue Hao và cộng sự (2021). Thiết kế này giúp biên dạng vòi phun tương thích tối ưu với các đặc tính động lực học của dòng chất lưu nén được, từ đó giảm thiểu các tổn thất thủy lực bất thuận nghịch bên trong vòi phun so với kết cấu vòi phun hình nón phẳng truyền thống.

3. Mô hình dòng rối và phương pháp số nào được lựa chọn để giải hệ phương trình vi phân trong ANSYS Fluent?
Đồ án ứng dụng Phương pháp Thể tích Hữu hạn (Finite Volume Method - FVM) để rời rạc hóa hệ phương trình vi phân đạo hàm riêng Navier-Stokes. Về mô hình dòng rối, nhóm tác giả kích hoạt phương trình năng lượng và lựa chọn mô hình Realizable $k-\varepsilon$ (Realizable k-epsilon turbulence model), là mô hình phù hợp để mô phỏng chính xác các dòng chảy có gradien áp suất lớn, dòng cuộn xoáy và dòng tia nén tốc độ cao.

4. Hiệu suất của thiết bị ejector được định lượng bằng đại lượng nào và được tính toán ra sao?
Hiệu suất làm việc của ejector được đo lường thông qua hệ số cuốn (Entrainment ratio, ký hiệu là $U$). Đại lượng này được xác định bằng tỉ số giữa lưu lượng khối lượng của dòng thứ cấp (dòng bị cuốn, $m_e$) và lưu lượng khối lượng của dòng sơ cấp (dòng động lực, $m_m$) theo công thức: $$U = \frac{m_e}{m_m}$$

5. Quy trình kiểm tra tính độc lập của lưới tính toán (mesh independence study) được thực hiện như thế nào?
Nhóm nghiên cứu thiết lập và khảo sát 3 mô hình lưới với mật độ phần tử khác nhau, ký hiệu là M1, M2 và M3. Chất lượng từng bộ lưới được định lượng thông qua các chỉ số Aspect Ratio, Skewness, Orthogonal Quality và Smoothness. Tính độc lập của lưới được xác nhận khi thặng dư sai số tính toán (residual index) đạt mức hội tụ và biên dạng phân bố vận tốc dọc trục (velocity profile) giữa các cấp lưới không còn sự sai lệch đáng kể.

Kết luận

Đồ án tốt nghiệp của nhóm tác giả Trần Đình Trung Đức, Đặng Minh Đức và Võ Tiến Dũng đã thực hiện thành công nghiên cứu mô phỏng số đặc tính dòng chảy và truyền nhiệt bên trong thiết bị ejector cải tiến sử dụng môi chất R134a. Bằng việc ứng dụng phương pháp thể tích hữu hạn trên phần mềm ANSYS Fluent với mô hình rối Realizable $k-\varepsilon$ và biên dạng vòi phun Witozinsky, nghiên cứu đã làm sáng tỏ cấu trúc trường dòng dưới âm và siêu âm, quy luật phân bố áp suất, nhiệt độ cũng như sự biến thiên của hệ số cuốn $U$ theo các điều kiện vận hành. Kết quả đạt được không chỉ kiểm chứng sự phù hợp với cơ sở lý thuyết khí động học dòng nén được mà còn cung cấp cơ sở dữ liệu kỹ thuật hữu ích phục vụ việc tính toán, thiết kế và tối ưu hóa hệ thống làm lạnh ejector tận dụng nguồn nhiệt thải.