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

Truyền nhiệt trong quá trình sôi hai pha đóng vai trò thiết yếu trong kỹ thuật Nhiệt - Điện lạnh nhờ hiệu suất trao đổi nhiệt vượt trội so với truyền nhiệt đơn pha. Hiện tượng chuyển pha từ lỏng sang hơi cho phép hệ thống hấp thụ và giải phóng một lượng nhiệt ẩn hóa hơi rất lớn, đồng thời duy trì khả năng kiểm soát nhiệt độ bề mặt ổn định và an toàn. Đặc tính này làm cho quá trình sôi trở thành giải pháp tản nhiệt then chốt trong các hệ thống đòi hỏi mật độ dòng nhiệt cao như thanh nhiên liệu trong lò phản ứng hạt nhân, thiết bị vi điện tử công suất lớn và các chu trình làm lạnh công nghiệp.

Tuy nhiên, việc nghiên cứu thực nghiệm quá trình sôi thường gặp nhiều trở ngại do chi phí chế tạo mô hình tốn kém, khó khăn trong việc đo đạc chi tiết các trường vật lý vi mô (phân bố nhiệt độ, vận tốc, áp suất cục bộ xung quanh bọt hơi) và tiềm ẩn các rủi ro vận hành ở điều kiện nhiệt độ, áp suất cao. Phương pháp mô phỏng số động lực học lưu chất (CFD) trên phần mềm ANSYS Workbench ra đời như một công cụ hỗ trợ cần thiết, cho phép khảo sát sâu cơ chế truyền nhiệt vi mô mà không bị giới hạn về không gian và thời gian. Trên cơ sở đó, nhóm tác giả đã thực hiện đề tài "Nghiên cứu mô phỏng số ảnh hưởng của hình dạng bề mặt đến khả năng trao đổi nhiệt của quá trình sôi bể (pool boiling)".

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

  • Mục tiêu nghiên cứu:

    1. Khảo sát ảnh hưởng của sự chênh lệch độ quá nhiệt trong dải $\Delta T = 5 \div 30^\circ\text{C}$ đến hệ số tỏa nhiệt đối lưu.
    2. Đánh giá ảnh hưởng của sức căng bề mặt $\sigma$ khi tăng từ $0,059\text{ N/m}$ lên $0,072\text{ N/m}$.
    3. Đánh giá ảnh hưởng của góc tiếp xúc bề mặt $\theta$ khi thay đổi từ $75^\circ$ lên $90^\circ$.
    4. Xác định tác động của việc bổ sung hạt nano $\text{SiO}_2$ vào nước với các tỷ lệ nồng độ phần trăm thể tích ($0,01%$; $0,02%$; $0,05%$; $0,06%$) đến hiệu suất trao đổi nhiệt và mức độ quá nhiệt bề mặt vách.
    5. Khảo sát ảnh hưởng của việc thay đổi cấu hình bề mặt gia nhiệt (vị trí, hình dạng và vật liệu) đến hiệu suất truyền nhiệt tổng thể.
    6. Cung cấp dữ liệu mô phỏng hỗ trợ việc cải thiện thiết kế, lựa chọn thông số vận hành tối ưu cho các thiết bị trao đổi nhiệt dạng sôi bể.

  • Nhiệm vụ nghiên cứu:

    1. Thiết lập mô hình hình học và thực hiện mô phỏng số dựa trên các thông số ban đầu từ tài liệu tham khảo.
    2. So sánh, đánh giá kết quả thu được từ mô phỏng CFD với các công bố khoa học liên quan để kiểm chứng độ tin cậy.
    3. Phân tích định lượng ảnh hưởng của các yếu tố hình học, nồng độ thể tích hạt nano và đặc tính lưu chất đến hiệu suất trao đổi nhiệt.
    4. Tổng hợp nhận xét, rút ra kết luận khoa học và đề xuất kiến nghị kỹ thuật.

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

    • Đối tượng nghiên cứu: Hiện tượng sôi bể (pool boiling) của lưu chất ở trạng thái tĩnh tiếp xúc với bề mặt gia nhiệt dạng tấm phẳng và dạng ống trơn nằm ngang; lưu chất nghiên cứu gồm nước nguyên chất và dung dịch hạt nano $\text{SiO}_2$/nước.
    • Phạm vi nghiên cứu: Khảo sát mô phỏng số 2D đối với bề mặt tấm phẳng làm bằng thép đặt ở đáy bể; mô phỏng số 3D đối với bề mặt ống trơn làm bằng đồng đặt ở giữa bể; dải độ quá nhiệt khảo sát $\Delta T = 5 \div 30^\circ\text{C}$; góc tiếp xúc $\theta = 75^\circ \div 90^\circ$; sức căng bề mặt $\sigma = 0,059 \div 0,072\text{ N/m}$; công cụ thực hiện là mô-đun Fluent thuộc phần mềm ANSYS Workbench 19.

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

Khóa luận xây dựng hệ thống cơ sở lý thuyết dựa trên các nguyên lý nhiệt động lực học và thủy động lực học của quá trình sôi bọt và sôi màng:

  • Cơ chế hình thành và phát triển bọt hơi: Điều kiện cân bằng cơ học của bọt hơi hình cầu được mô tả qua phương trình Laplace: $$\Delta p = p_h - p_f = \frac{2\sigma}{R_o}$$ Trong đó $p_h$ là áp suất hơi bão hòa bên trong bọt, $p_f$ là áp suất chất lỏng bao quanh, $\sigma$ là sức căng bề mặt và $R_o$ là bán kính tới hạn của bọt hơi. Kết hợp với phương trình Clausius-Clapeyron, bán kính tới hạn được xác định theo độ quá nhiệt $\Delta T$: $$R_o = \frac{2\sigma T_s}{r \rho_h \Delta T}$$ Theo tính toán lý thuyết ở áp suất khí quyển ($T_s = 100^\circ\text{C}$), khi $\Delta T = 5^\circ\text{C}$ thì $R_o = 6,7,\mu\text{m}$; khi $\Delta T = 25^\circ\text{C}$ thì $R_o$ giảm xuống còn $1,3,\mu\text{m}$. Tốc độ phát triển bán kính bọt hơi theo thời gian $\tau$ tuân theo quan hệ $R = 3,5 \sqrt{\frac{\lambda \Delta T \tau}{r \rho_h}}$. Đường kính tách ly bọt hơi $D_t$ được tính theo công thức cân bằng lực nâng và sức căng bề mặt: $D_t = 0,0208,\theta \sqrt{\frac{\sigma}{g(\rho_f - \rho_h)}}$. Chu kỳ sản sinh bọt hơi của nước ở áp suất tiêu chuẩn là $\tau_o = 0,05\text{ s}$, với tốc độ nổi của bọt hơi lên mặt thoáng dao động từ $0,17\text{ m/s}$ đến $0,35\text{ m/s}$.

  • Đặc tính đường cong sôi (Nukiyama, 1934): Đường cong biểu diễn mối liên hệ giữa mật độ dòng nhiệt $q$, hệ số tỏa nhiệt đối lưu $\alpha$ và độ quá nhiệt $\Delta T$ phân chia thành 4 vùng rõ rệt: Vùng đối lưu tự nhiên ($\Delta T < 5^\circ\text{C}$, $T_w < T_{IB}$); Vùng sôi bọt ổn định ($\Delta T > 5^\circ\text{C}$, $T_{IB} < T_w < T_{DNB}$); Vùng chuyển tiếp không ổn định ($22,2^\circ\text{C} < \Delta T < 100^\circ\text{C}$, $T_{DNB} < T_w < T_{MFB}$) nơi bọt hơi liên kết tạo màng cản trở truyền nhiệt; và Vùng sôi màng ổn định ($\Delta T > 110^\circ\text{C}$, $T_w > T_{MFB}$) nơi bức xạ nhiệt đóng vai trò chủ đạo.

  • Mô hình dòng chảy đa pha và mô hình rối: Đề tài sử dụng mô hình đa pha Eulerian (Eulerian Multiphase Model) thay vì mô hình Volume of Fluid (VOF) dựa trên sự phân tích so sánh kỹ thuật:

Tiêu chí Mô hình Eulerian Mô hình VOF
Ưu điểm Dễ thực hiện phép tính toán học; hiệu quả cao cho dòng đa pha không nén được; mô phỏng ổn định tương tác giữa các pha liên tục và phân tán. Theo dõi mặt phân chia pha (interface) rõ ràng; phù hợp cho các bài toán mặt thoáng và dòng nén được.
Nhược điểm Khó theo dõi sắc nét bề mặt phân chia pha; chi phí bộ nhớ tăng theo số pha. Khối lượng tính toán lớn; dễ xảy ra hiện tượng "xé lưới" tại bề mặt phân chia pha phức tạp.
Giới hạn kỹ thuật Mô hình ứng suất Reynolds không hỗ trợ riêng cho từng pha; hạt phân tán chỉ tương tác với pha chính. Yêu cầu bộ giải Pressure-Based; không cho phép vùng trống không có chất lỏng; chỉ một pha được định nghĩa khí lý tưởng.

Trong mô hình Eulerian, hệ thống phương trình điều khiển bao gồm phương trình bảo toàn khối lượng (liên tục), phương trình bảo toàn động lượng và phương trình bảo toàn năng lượng được giải riêng biệt cho từng pha $q$, đồng thời chia sẻ chung một trường áp suất $p$. Trạng thái chảy rối trong hỗn hợp sôi được mô hình hóa bằng mô hình rối $k-\varepsilon$ hỗn hợp (Mixture Turbulence Model) để tối ưu hóa thời gian tính toán và tài nguyên phần cứng.

  • Phương pháp giải số và kiểm tra hội tụ: Miền tính toán được rời rạc hóa bằng phương pháp thể tích hữu hạn (Finite Volume Method - FVM) kết hợp bộ giải dựa trên áp suất (Pressure-Based Solver). Nhóm nghiên cứu tiến hành kiểm nghiệm tính độc lập của lưới (Mesh-independent test) thông qua 3 cấu trúc lưới M1, M2, M3; đồng thời theo dõi độ dư hội tụ (Residuals) và so sánh profile phân bố vận tốc qua các mặt cắt chuẩn để chọn mật độ lưới tối ưu.

  • Nguồn dữ liệu: Dữ liệu hình học, điều kiện biên và tính chất nhiệt vật lý được trích xuất và đối chiếu từ hai công bố khoa học quốc tế: bài báo của tác giả nghiên cứu về mô phỏng sôi bọt với chất lỏng nano $\text{SiO}_2$/nước ("Numerical Simulation of Nucleate Boiling Heat Transfer Performance using SiO2/Water Nanofluid") và nghiên cứu thực nghiệm sôi bể trên ống trơn của Abhishek Kumar và Shou Yin Yang ("Pool boiling heat transfer performance of low-GWP refrigerant R-513A on smooth tube"), kết hợp bảng tra cứu nhiệt vật lý chuẩn của nước bão hòa ở $100^\circ\text{C}$, hơi nước, thép và đồng.

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

Chương 1. Tổng quan

Chương 1 trình bày lý do lựa chọn đề tài, tính cấp thiết của việc nghiên cứu truyền nhiệt hai pha trong bối cảnh phát triển của ngành Kỹ thuật Nhiệt - Điện lạnh. Tác giả phân tích bản chất vật lý của quá trình sôi, so sánh chi tiết giữa sôi bể (pool boiling - chất lỏng tĩnh, bọt hơi lớn hình thành trên bề mặt gia nhiệt) và sôi dòng (flow boiling - chất lỏng chuyển động cưỡng bức trong ống, bọt hơi nhỏ bị cuốn theo dòng). Chương này tổng hợp các công trình nghiên cứu trong nước (nghiên cứu kênh vi mô, ảnh hưởng của gờ truyền nhiệt, nano chất lỏng) và quốc tế (các công trình của M. Mahian, Alessandro Sciacca, Abhishek Kumar). Đồng thời, tác giả luận giải việc lựa chọn mô hình đa pha Eulerian và mô hình rối $k-\varepsilon$ hỗn hợp dựa trên ưu thế về độ ổn định tính toán đối với dòng sôi hai pha.

Chương 2. Cơ sở lý thuyết

Nội dung chương 2 tập trung vào hệ thống phương trình toán học mô tả cơ chế sôi. Tác giả phân tích các yếu tố ảnh hưởng trực tiếp đến quá trình truyền nhiệt:

  • Độ quá nhiệt $\Delta T$: Tác động đến sự chuyển đổi giữa các chế độ đối lưu tự nhiên, sôi bọt, sôi màng không ổn định và sôi màng ổn định.
  • Góc tiếp xúc $\theta$: Phân biệt bề mặt dính ướt ($\theta < 90^\circ$, bọt hơi dễ tách ly) và không dính ướt ($\theta > 90^\circ$, bọt hơi trải rộng, dễ tạo màng hơi cản nhiệt).
  • Sức căng bề mặt $\sigma$ và độ nhớt $\mu$: Mối quan hệ tỷ lệ nghịch với hệ số tỏa nhiệt đối lưu qua các công thức thực nghiệm $\alpha \sim \sigma^{-0,33}$ và $\alpha \sim \mu^{-0,45}$.
  • Trạng thái bề mặt và cấu trúc dòng: Ảnh hưởng của độ nhám bề mặt đến số lượng tâm hóa hơi; cấu trúc dòng và sự phân bố nhiệt độ trong ống nằm ngang và ống đứng. Chương này cũng xác lập hệ thống phương trình vi phân nhiều pha Eulerian (phương trình liên tục, động lượng, năng lượng) và phương trình truyền vận động năng rối $k$ cùng tỷ lệ tiêu tán $\varepsilon$.

Chương 3. Giới thiệu về mô phỏng CFD

Chương 3 khái quát nguyên lý tính toán động lực học lưu chất, phân loại các bài toán CFD cơ bản (trao đổi nhiệt, dòng nén được, quá trình đốt cháy). Tác giả chi tiết hóa quy trình mô phỏng CFD gồm quy trình cơ bản cho người dùng (8 bước: từ xây dựng hình học, thiết lập thông số, chạy mô phỏng đến phân tích kết quả) và quy trình nâng cao cho kỹ sư gồm hai khâu trọng yếu:

  • Quy trình kiểm tra (Verification): Khảo sát độc lập lưới nhằm triệt tiêu sai số do phương pháp số.
  • Quy trình kiểm nghiệm (Validation): So sánh số liệu mô phỏng với dữ liệu thực nghiệm hoặc kết quả giải tích. Chương này cũng giới thiệu cấu trúc các mô-đun của ANSYS Fluent 19 (Geometry Design Modeler, Meshing, Setup, Solution, Results), phương pháp thể tích hữu hạn (FVM) và các tiêu chuẩn đánh giá chất lượng lưới (độ xiên - skewness, tỷ lệ co giãn - aspect ratio).

Chương 4. Xây dựng mô hình

Chương 4 mô tả toàn bộ quá trình thiết lập bài toán mô phỏng cho hai mô hình hình học:

  • Mô hình 1 (2D): Hình học dạng bể chữ nhật, bề mặt gia nhiệt dạng tấm phẳng bằng thép đặt tại đáy; lưu chất là nước hoặc nanofluid $\text{SiO}_2$/nước.
  • Mô hình 2 (3D): Bể chứa dạng khối hộp, bề mặt gia nhiệt là ống trơn bằng đồng đặt nằm ngang ở vị trí giữa bể; lưu chất là nước. Tác giả trình bày phương pháp chia lưới cho cả hai mô hình, thiết lập các bộ lưới thử nghiệm M1, M2, M3 để kiểm tra tính độc lập của lưới. Các thông số vật lý của nước lỏng ở $25^\circ\text{C}$ và khi sôi ở $100^\circ\text{C}$, hơi nước bão hòa, thép, đồng và các nồng độ hạt nano $\text{SiO}_2$ ($0,01% \div 0,06%$) được cài đặt chi tiết trong môi trường Fluent Setup cùng các điều kiện biên nhiệt độ và áp suất tương ứng.

Chương 5. Kết quả mô phỏng

Chương 5 phân tích toàn bộ dữ liệu định lượng thu được từ các ca mô phỏng:

  • Đánh giá sự hội tụ lưới thông qua đồ thị Residuals và đường phân bố vector vận tốc của 3 bộ lưới M1, M2, M3, qua đó xác định bộ lưới tối ưu đảm bảo độ chính xác tính toán.
  • Khảo sát trường phân bố nhiệt độ, áp suất, vận tốc và đường biểu diễn thông lượng nhiệt $q$, hệ số tỏa nhiệt $\alpha$ cho mô hình 1 ở các độ quá nhiệt $\Delta T = 25^\circ\text{C}$ và $\Delta T = 30^\circ\text{C}$.
  • Đánh giá định lượng ảnh hưởng của góc tiếp xúc $\theta$ ($75^\circ \to 90^\circ$) và sức căng bề mặt $\sigma$ ($0,059 \to 0,072\text{ N/m}$) trên mô hình 1.
  • Phân tích hiệu quả tăng cường truyền nhiệt khi thay đổi nồng độ phần trăm thể tích hạt nano $\text{SiO}_2$ ($0,01%$, $0,02%$, $0,05%$, $0,06%$).
  • Trình bày kết quả mô phỏng 3D trên mô hình 2 tại $\Delta T = 15^\circ\text{C}$ qua hai mặt cắt khảo sát, phân tích sự dao động của vector vận tốc, phân bố nhiệt độ, áp suất và so sánh độ lệch hiệu suất giữa mô hình 1 và mô hình 2.

Chương 6. Kết luận và kiến nghị

Chương cuối cùng tổng kết các quy luật nhiệt - thủy động đã được làm sáng tỏ qua mô phỏng số, khẳng định tính khả thi của phương pháp CFD trong nghiên cứu sôi bể và đề xuất các khuyến nghị về việc hoàn thiện mô hình hình học cũng như mở rộng hướng nghiên cứu thực nghiệm.

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

  • Kết quả mô phỏng chính: Các kết quả định lượng từ đồ án chỉ ra rõ rệt các quy luật trao đổi nhiệt trong quá trình sôi bể:
Yếu tố khảo sát Dải biến thiên Ảnh hưởng đến hiệu suất trao đổi nhiệt ($q$ và $\alpha$)
Độ quá nhiệt ($\Delta T$) $5^\circ\text{C} \div 30^\circ\text{C}$ Hiệu suất trao đổi nhiệt, thông lượng nhiệt và hệ số tỏa nhiệt đối lưu tăng mạnh theo sự gia tăng của $\Delta T$.
Góc tiếp xúc ($\theta$) $75^\circ \to 90^\circ$ Hiệu suất trao đổi nhiệt giảm khi góc tiếp xúc tăng, do bọt hơi khó tách ly khỏi bề mặt gia nhiệt.
Sức căng bề mặt ($\sigma$) $0,059\text{ N/m} \to 0,072\text{ N/m}$ Hiệu suất trao đổi nhiệt có sự thay đổi nhưng mức độ biến thiên không đáng kể.
Nồng độ nanofluid $\text{SiO}_2$ $0,01% \div 0,06%$ thể tích Hiệu suất trao đổi nhiệt tăng đáng kể so với việc sử dụng nước nguyên chất ở cùng điều kiện làm việc.
Cấu hình bề mặt gia nhiệt Tấm phẳng thép đáy bể (2D) $\to$ Ống đồng giữa bể (3D) Thay đổi hình dạng sang ống trơn, đổi vật liệu sang đồng và nâng vị trí lên giữa bể làm tăng đáng kể hiệu suất trao đổi nhiệt.

  • Giải pháp và kiến nghị kỹ thuật:

    • Trong thiết kế thiết bị trao đổi nhiệt dạng sôi bể, nên ưu tiên sử dụng bề mặt gia nhiệt có tính dính ướt cao (góc tiếp xúc $\theta$ nhỏ) để bọt hơi dễ dàng tách ly, giải phóng tâm sôi cho chu kỳ tiếp theo.
    • Ứng dụng dung dịch nanofluid $\text{SiO}_2$ với nồng độ thể tích phù hợp để cải thiện đáng kể hệ số dẫn nhiệt hiệu dụng của môi chất công tác.
    • Bố trí bề mặt gia nhiệt dạng ống đặt ở giữa bể và sử dụng vật liệu có độ dẫn nhiệt cao như đồng thay cho thép tấm phẳng ở đáy nhằm thúc đẩy đối lưu tự nhiên và giải tỏa bọt hơi nhanh hơn.

  • Đóng góp của đồ án:

    • Thiết lập hoàn chỉnh quy trình mô phỏng số dòng sôi hai pha sử dụng mô hình Eulerian kết hợp mô hình rối $k-\varepsilon$ hỗn hợp trên phần mềm ANSYS Fluent 19.
    • Mở rộng phạm vi nghiên cứu từ mô hình 2D tấm phẳng sang mô hình 3D ống trơn nằm ngang, đồng thời khảo sát toàn diện tác động đồng thời của các thông số nhiệt động vi mô (sức căng bề mặt, góc tiếp xúc, độ quá nhiệt) và vật liệu mới (nanofluid $\text{SiO}_2$).

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

  • Hạn chế của đề tài:

    • Nghiên cứu thuần túy dựa trên phương pháp mô phỏng số CFD, các thông số đối chứng được lấy gián tiếp từ bài báo khoa học công bố trước đó mà chưa có điều kiện thực hiện đo đạc thực nghiệm song song để kiểm chứng trực tiếp.
    • Mô hình nanofluid mới chỉ khảo sát đến nồng độ thể tích tối đa là $0,06%$, chưa mở rộng đến dải nồng độ cao hơn để xác định điểm suy giảm truyền nhiệt do hiện tượng lắng đọng hạt nano.
    • Phạm vi khảo sát giới hạn ở trạng thái sôi bể tĩnh, chưa phân tích tương tác trong điều kiện có dòng chảy cưỡng bức (flow boiling).

  • Hướng nghiên cứu tiếp:

    • Thiết kế hệ thống thực nghiệm thực tế để kiểm chứng và hiệu chỉnh độ chính xác của các hệ số trong mô hình mô phỏng số.
    • Mở rộng dải khảo sát nồng độ hạt nano $\text{SiO}_2$ và thử nghiệm thêm các loại hạt nano oxit kim loại khác (như $\text{Al}_2\text{O}_3$, $\text{CuO}$).
    • Nghiên cứu mô phỏng trên các bề mặt gia nhiệt có cấu trúc vi mô phức tạp (bề mặt có vi rãnh, vi gờ, cấu trúc xốp) nhằm tìm kiếm giải pháp nâng cao hơn nữa thông lượng nhiệt tới hạn (CHF).

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:

  • Đối tượng tham khảo: Sinh viên, học viên cao học và kỹ sư thuộc các chuyên ngành Công nghệ Kỹ thuật Nhiệt, Kỹ thuật Năng lượng, Cơ khí Động lực và các kỹ sư R&D thiết bị trao đổi nhiệt, hệ thống làm mát trung tâm dữ liệu hoặc lò phản ứng.
  • Nội dung đáng giá:
    • Chương 3 và Chương 4 cung cấp quy trình từng bước thiết lập mô hình đa pha Eulerian và kiểm tra tính độc lập của lưới trong ANSYS Fluent 19.
    • Hệ thống công thức nhiệt động lực học bọt hơi tại Chương 2 giúp người đọc nắm vững phương pháp tính toán bán kính tới hạn, đường kính tách ly và chu kỳ phát triển bọt hơi.
    • Các số liệu định lượng tại Chương 5 đóng vai trò làm tài liệu đối sánh cho các nghiên cứu mô phỏng số liên quan đến sôi bể và chất lỏng nano.

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

1. Vì sao nhóm tác giả lựa chọn mô hình đa pha Eulerian thay vì mô hình VOF cho bài toán sôi bể?
Mô hình Eulerian có khả năng giải quyết hiệu quả các bài toán dòng chảy nhiều pha không nén được với sự tương tác mạnh mẽ giữa pha lỏng liên tục và pha hơi phân tán. Mặc dù VOF theo dõi mặt phân chia pha rõ nét hơn nhưng lại đòi hỏi khối lượng tính toán rất lớn và dễ gặp hiện tượng "xé lưới" tại các giao diện bọt hơi phức tạp. Mô hình Eulerian giải hệ phương trình động lượng, liên tục và năng lượng riêng cho từng pha trên cùng một trường áp suất, mang lại độ ổn định tính toán và khả năng hội tụ tốt hơn cho bài toán truyền nhiệt sôi.

2. Bán kính tới hạn của bọt hơi thay đổi như thế nào khi độ quá nhiệt tăng từ $5^\circ\text{C}$ lên $25^\circ\text{C}$?
Theo công thức tính bán kính tới hạn $R_o = \frac{2\sigma T_s}{r \rho_h \Delta T}$ tại áp suất khí quyển ($T_s = 100^\circ\text{C}$), bán kính tới hạn tỷ lệ nghịch với độ quá nhiệt $\Delta T$. Khi $\Delta T = 5^\circ\text{C}$, bán kính tới hạn là $R_o = 6,7,\mu\text{m}$. Khi tăng độ quá nhiệt lên $\Delta T = 25^\circ\text{C}$, bán kính tới hạn giảm mạnh xuống chỉ còn $1,3,\mu\text{m}$, giúp các tâm hóa hơi kích thước nhỏ trên bề mặt dễ dàng kích hoạt và tạo mầm sôi nhanh hơn.

3. Việc tăng góc tiếp xúc từ $75^\circ$ lên $90^\circ$ ảnh hưởng thế nào đến hiệu suất truyền nhiệt?
Khi góc tiếp xúc $\theta$ tăng từ $75^\circ$ lên $90^\circ$, tính dính ướt của chất lỏng trên bề mặt gia nhiệt giảm đi. Đường kính tách ly của bọt hơi tỷ lệ thuận với góc tiếp xúc ($D_t \sim \theta$), khiến chân bọt hơi trải rộng trên bề mặt đốt nóng, bọt hơi khó tách ly và có xu hướng liên kết lại tạo thành màng hơi cục bộ. Lớp hơi này có độ dẫn nhiệt kém hơn nước rất nhiều, dẫn đến hiệu suất trao đổi nhiệt tổng thể bị suy giảm.

4. Hai mô hình hình học trong đồ án khác nhau ở những điểm cốt lõi nào?
Mô hình 1 được thiết lập dưới dạng 2D, sử dụng bề mặt gia nhiệt dạng tấm phẳng làm bằng thép đặt tại đáy bể chứa, khảo sát cho cả nước và nanofluid $\text{SiO}_2$ ở dải $\Delta T = 5 \div 30^\circ\text{C}$. Mô hình 2 được xây dựng dưới dạng không gian 3D, thay đổi bề mặt gia nhiệt thành dạng ống trơn nằm ngang, đổi vật liệu thành đồng (có hệ số dẫn nhiệt cao hơn thép), chuyển vị trí gia nhiệt lên giữa bể và khảo sát truyền nhiệt với môi chất nước.

5. Việc bổ sung hạt nano $\text{SiO}_2$ vào nước tác động ra sao đến quá trình sôi?
Hỗn hợp nanofluid nước/$\text{SiO}_2$ với các tỷ lệ thể tích từ $0,01%$ đến $0,06%$ giúp tăng cường các đặc tính nhiệt vật lý của môi chất, đặc biệt là hệ số dẫn nhiệt hiệu dụng. Kết quả mô phỏng số khẳng định việc sử dụng nanofluid làm tăng đáng kể hệ số tỏa nhiệt đối lưu và thông lượng nhiệt so với nước nguyên chất ở cùng độ chênh lệch nhiệt độ.

Kết luận

Đồ án tốt nghiệp của nhóm tác giả tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh đã ứng dụng thành công công cụ mô phỏng số CFD (ANSYS Fluent 19) để phân tích chi tiết cơ chế truyền nhiệt sôi bể hai pha. Công trình đã làm sáng tỏ mối liên hệ định lượng giữa hiệu suất trao đổi nhiệt với độ quá nhiệt, góc tiếp xúc, sức căng bề mặt, nồng độ nanofluid $\text{SiO}_2$ và cấu hình hình học bề mặt gia nhiệt. Các phát hiện thu được không chỉ củng cố cơ sở lý thuyết nhiệt - thủy động mà còn cung cấp cơ sở dữ liệu quan trọng phục vụ công tác thiết kế và tối ưu hóa các thiết bị trao đổi nhiệt hiệu suất cao trong công nghiệp.