Giới thiệu dự án

Trong bối cảnh các ngành công nghiệp hiện đại như trung tâm dữ liệu (Data Center), hệ thống vi mạch điện tử công suất cao (GaN/SiC), lò phản ứng hạt nhân và hệ thống lạnh công nghiệp đang đối mặt với bài toán mật độ nhiệt phát sinh vượt ngưỡng $100\text{ W/cm}^2$, các kỹ thuật tản nhiệt đơn pha truyền thống dần chạm tới giới hạn vật lý. Quá trình trao đổi nhiệt sôi hai pha (two-phase boiling heat transfer), đặc biệt là cơ chế sôi bể (pool boiling), đóng vai trò then chốt nhờ khả năng hấp thụ và vận chuyển lượng nhiệt khổng lồ thông qua nhiệt ẩn hóa hơi ($h_{fg}$) với độ chênh lệch nhiệt độ tương đối nhỏ. Tuy nhiên, việc tối ưu hóa hiệu suất truyền nhiệt khi sôi đòi hỏi kiểm soát chặt chẽ các thông số bề mặt và tính chất vi mô của lưu chất.

Phương pháp nghiên cứu thực nghiệm nhiệt - thủy động lực học truyền thống thường tốn kém chi phí chế tạo mẫu thử, khó đo đạc chi tiết trường vận tốc, trường nhiệt độ cục bộ tại lớp biên sát bề mặt gia nhiệt và tiềm ẩn nguy cơ mất an toàn khi tiến gần đến giới hạn thông lượng nhiệt tới hạn (Critical Heat Flux - CHF). Nhằm giải quyết các thách thức trê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)" được thực hiện nhằm ứng dụng công nghệ mô phỏng động lực học dòng chảy (CFD) để khảo sát toàn diện các quy luật biến thiên nhiệt - thủy động.

+-----------------------------------------------------------------------------------+
|                              MỤC TIÊU DỰ ÁN TỔNG QUAN                             |
+-----------------------------------------------------------------------------------+
|  1. Xây dựng mô hình số 2D và 3D mô phỏng hiện tượng sôi bể đa pha trên ANSYS.   |
|  2. Khảo sát ảnh hưởng của độ quá nhiệt vách: Delta T = 5°C đến 30°C.             |
|  3. Đánh giá tác động của góc tiếp xúc (75° -> 90°) và sức căng bề mặt lưu chất.  |
|  4. Phân tích hiệu quả tăng cường truyền nhiệt của Nanofluid SiO2 (0.01% - 0.06%). |
|  5. Tối ưu hóa cấu hình bề mặt gia nhiệt: Tấm phẳng thép 2D vs Ống trơn đồng 3D.  |
+-----------------------------------------------------------------------------------+

Phương pháp tiếp cận của đề tài dựa trên việc ứng dụng phần mềm ANSYS Fluent giải hệ phương trình bảo toàn đa pha Eulerian kết hợp mô hình rối hỗn hợp ($k-\varepsilon$ Mixture Turbulence Model). Kết quả kỳ vọng cung cấp cơ sở dữ liệu định lượng về mật độ dòng nhiệt ($q$), hệ số tỏa nhiệt đối lưu ($\alpha$), động lực học phát triển bọt hơi phục vụ trực tiếp công tác thiết kế thiết bị trao đổi nhiệt (TBTĐN).

Phạm vi và giới hạn nghiên cứu:

  • Môi chất khảo sát: Nước nguyên chất bão hòa tại áp suất khí quyển ($p = 1\text{ atm}$, $T_s = 100^\circ\text{C}$) và hỗn hợp hạt nano $\text{SiO}_2$/nước với nồng độ thể tích từ $0.01%$ đến $0.06%$.
  • Chế độ sôi: Tập trung vào vùng đối lưu tự nhiên và vùng sôi bọt (nucleate boiling) với dải độ quá nhiệt $\Delta T = T_w - T_s = 5 \div 30^\circ\text{C}$.
  • Hình học gia nhiệt: Tấm phẳng đáy bể (2D, vật liệu thép) và ống trơn nằm ngang đặt giữa bể (3D, vật liệu đồng).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Nghiên cứu quá trình sôi bể hiện nay chủ yếu sử dụng ba hướng tiếp cận: Đo đạc thực nghiệm, mô phỏng thể tích chất lỏng (Volume of Fluid - VOF), và mô hình đa pha liên tục (Eulerian Multiphase Model).

Phương pháp Ưu điểm Nhược điểm Giới hạn ứng dụng
Thực nghiệm nhiệt truyền thống Cung cấp dữ liệu vật lý thực tế, độ tin cậy cao. Chi phí chế tạo cao, khó quan sát vi cấu trúc bọt hơi sát vách, nguy cơ cháy nổ khi đạt CHF. Đo đạc thông số vĩ mô tích phân.
Mô phỏng VOF (Volume of Fluid) Theo dõi sắc nét ranh giới tiếp xúc pha lỏng - hơi. Chi phí tính toán cực lớn, dễ mất ổn định lưới khi mật độ bọt hơi cao, khó scale-up hệ 3D. Mô phỏng đơn bọt hoặc nhóm bọt nhỏ.
Mô hình Eulerian (Được chọn) Giải phương trình bảo toàn riêng biệt cho từng pha, ổn định số cao, mô phỏng tốt dòng hạt nano. Không nắm bắt trực tiếp biến dạng ranh giới giao pha vi mô. Hệ thống công nghiệp đa pha, sôi bể diện tích lớn.

Phân loại yêu cầu hệ thống mô phỏng theo mô hình MoSCoW:

  • Must have (Bắt buộc): Mô hình hóa chính xác sự chuyển pha lỏng - hơi; giải thuật hội tụ cho dải độ quá nhiệt $\Delta T = 5 \div 30^\circ\text{C}$; kiểm nghiệm độc lập lưới (Grid Independence Test).
  • Should have (Nên có): Tích hợp tính chất nhiệt vật lý của hạt nano $\text{SiO}_2$ ở các tỷ lệ thể tích khác nhau; khảo sát thay đổi góc dính ướt $\theta = 75^\circ \div 90^\circ$.
  • Could have (Có thể có): Mở rộng khảo sát hình học 3D từ tấm phẳng đáy sang ống tròn nằm ngang vật liệu đồng.
  • Won't have (Chưa thực hiện): Mô phỏng chi tiết độ nhám bề mặt fractal ngẫu nhiên và phản ứng hóa học bám cặn bề mặt lâu dài.

Thiết kế hệ thống

Kiến trúc quy trình mô phỏng số CFD được chuẩn hóa theo quy trình chuẩn của ANSYS Workbench:

Technology Stack:

  • Phần mềm điều phối: ANSYS Workbench 19.0
  • Bộ tiền xử lý hình học: ANSYS DesignModeler 19.0
  • Bộ tiền xử lý chia lưới: ANSYS Meshing 19.0 (Cấu trúc Lưới đa tầng Quad/Hexahedral kết hợp Inflation Layer)
  • Bộ giải số: ANSYS Fluent 19.0 (Pressure-Based Solver, Phase-Coupled SIMPLE)
  • Trình biên dịch User-Defined Functions (UDF): Microsoft Visual Studio C++ Compiler

Phương trình điều khiển và Cơ sở toán học:

  1. Phương trình liên tục cho pha $q$ ($q = l$ cho lỏng, $q = v$ cho hơi): $$\frac{\partial}{\partial t}(\alpha_q \rho_q) + \nabla \cdot (\alpha_q \rho_q \vec{v}q) = \sum{p=1}^{n} (\dot{m}{pq} - \dot{m}{qp})$$

  2. Phương trình động lượng: $$\frac{\partial}{\partial t}(\alpha_q \rho_q \vec{v}q) + \nabla \cdot (\alpha_q \rho_q \vec{v}q \vec{v}q) = -\alpha_q \nabla p + \nabla \cdot \bar{\bar{\tau}}q + \alpha_q \rho_q \vec{g} + \sum{p=1}^{n} (\vec{R}{pq} + \dot{m}{pq}\vec{v}{pq} - \dot{m}{qp}\vec{v}{qp}) + \vec{F}{lift,q} + \vec{F}{vm,q}$$

  3. Động lực học mầm bọt hơi (Phương trình Laplace): $$\Delta p = p_h - p_f = \frac{2\sigma}{R_o} \implies R_o = \frac{2\sigma T_s}{r \rho_h \Delta T}$$ Trong đó $R_o$ là bán kính tới hạn của mầm bọt, $\sigma$ là sức căng bề mặt ($\text{N/m}$), $r$ ($h_{fg}$) là nhiệt ẩn hóa hơi ($\text{J/kg}$), $\rho_h$ là khối lượng riêng của hơi ($\text{kg/m}^3$).

  4. Đường kính bọt hơi khi tách ly bề mặt: $$D_t = 0.0208 \cdot \theta \sqrt{\frac{\sigma}{g(\rho_f - \rho_h)}}$$ Với $\theta$ là góc tiếp xúc (độ dính ướt), $g$ là gia tốc trọng trường ($9.81\text{ m/s}^2$).

/* Đoạn mã UDF C-language định nghĩa nguồn truyền khối chuyển pha sôi lỏng-hơi */
#include "udf.h"

#define LATENT_HEAT 2.257e6 /* Nhiet an hoa hoi cua nuoc (J/kg) */
#define T_SAT 373.15        /* Nhiet do bao hoa tai 1 atm (K) */

DEFINE_SOURCE(mass_transfer_boil, cell, thread, dS, eqn)
{
    real temp = C_T(cell, thread);
    real source = 0.0;
    
    if (temp > T_SAT) {
        /* Toc do sinh hoi ti le voi do qua nhiet */
        source = 0.1 * C_R(cell, thread) * (temp - T_SAT) / T_SAT;
        dS[eqn] = 0.1 * C_R(cell, thread) / T_SAT;
    } else {
        source = 0.0;
        dS[eqn] = 0.0;
    }
    return source;
}

Methodology

Quy trình nghiên cứu áp dụng khung Verification & Validation (V&V) tiêu chuẩn công nghiệp:

  1. Kiểm nghiệm mô hình số (Verification): Thiết lập 3 bộ lưới phân giải tăng dần ($M_1$: thô, $M_2$: trung bình, $M_3$: mịn). Khảo sát profin vận tốc và nhiệt độ dọc theo đường chuẩn để loại trừ hiện tượng phụ thuộc lưới.
  2. Kế hoạch thực hiện (Milestones):
    • Tháng 1-2: Xây dựng mô hình hình học 2D tấm phẳng, phân tích lý thuyết, thiết lập bài toán cơ sở với nước nguyên chất.
    • Tháng 3-4: Chạy mô phỏng khảo sát ảnh hưởng của góc tiếp xúc $\theta$, sức căng bề mặt $\sigma$, và nồng độ thể tích hạt nano $\text{SiO}_2$ ($0.01% \div 0.06%$).
    • Tháng 5-6: Nâng cấp mô hình lên 3D ống trơn nằm ngang vật liệu đồng, phân tích đối chiếu và tối ưu hóa thông số.
  3. Quản trị rủi ro số: Sử dụng hệ số giãn sai số (Under-Relaxation Factors) phù hợp (Pressure: 0.3, Momentum: 0.7, Energy: 1.0) nhằm ngăn ngừa phân kỳ nghiệm do số Courant tăng cao trong pha hơi.

Implementation và kết quả

Development process

Quá trình mô phỏng được triển khai qua các bước cấu hình vật lý chi tiết:

  • Thiết lập không gian tính toán:
    • Mô hình 1 (2D): Bể chứa kích thước $100\text{ mm} \times 150\text{ mm}$, đáy bể là tấm phẳng gia nhiệt bằng thép ($k = 16.27\text{ W/m.K}$), biên bên là vách đoạn nhiệt, biên đỉnh là áp suất thoát ($p_{out} = 0\text{ Pa gauge}$).
    • Mô hình 2 (3D): Bể chứa chất lỏng tĩnh có đặt ống trơn nằm ngang bằng đồng ($k = 387.6\text{ W/m.K}$) đường kính $D = 19\text{ mm}$ tại vị trí tâm bể.
  • Thiết lập tính chất Nanofluid $\text{SiO}_2$/Nước:
    • Khối lượng riêng hỗn hợp: $\rho_{nf} = (1 - \phi)\rho_{bf} + \phi \rho_{np}$
    • Nhiệt dung riêng: $(\rho C_p){nf} = (1 - \phi)(\rho C_p){bf} + \phi (\rho C_p)_{np}$
    • Độ dẫn nhiệt hiệu dụng (Mô hình Maxwell): $k_{nf} = k_{bf} \left[ \frac{k_{np} + 2k_{bf} - 2\phi(k_{bf} - k_{np})}{k_{np} + 2k_{bf} + \phi(k_{bf} - k_{np})} \right]$

Testing và validation

1. Kiểm nghiệm độc lập lưới (Grid Independence Test)

Khảo sát sự thay đổi profin vận tốc dọc trục trung tâm với 3 cấu hình lưới:

  • Lưới $M_1$: 45,200 phần tử (Sai số vận tốc so với thực nghiệm: $6.8%$)
  • Lưới $M_2$: 98,600 phần tử (Sai số vận tốc: $2.1%$)
  • Lưới $M_3$: 185,400 phần tử (Sai số vận tốc so với $M_2$: $< 0.8%$) $\to$ Cấu hình lưới $M_2$ được lựa chọn nhằm tối ưu hóa giữa độ chính xác vật lý và tài nguyên tính toán.
Tiêu chí hội tụ phần dư (Residual Convergence Criteria):
+----------------------+--------------------+--------------------+
| Phương trình         | Ngưỡng hội tụ      | Trạng thái nghiệm  |
+----------------------+--------------------+--------------------+
| Continuity           | < 1.0e-04          | Hội tụ ổn định     |
| X-Velocity           | < 1.0e-04          | Hội tụ ổn định     |
| Y-Velocity           | < 1.0e-04          | Hội tụ ổn định     |
| Energy               | < 1.0e-06          | Hội tụ tuyệt đối   |
| k - epsilon          | < 1.0e-04          | Hội tụ ổn định     |
| Volume Fraction      | < 1.0e-04          | Hội tụ ổn định     |
+----------------------+--------------------+--------------------+

Kết quả đạt được

1. Ảnh hưởng của độ quá nhiệt $\Delta T$ ($5 \div 30^\circ\text{C}$) đối với nước nguyên chất

Khi tăng độ quá nhiệt vách $\Delta T$, mật độ dòng nhiệt ($q$) và hệ số tỏa nhiệt đối lưu ($\alpha$) tăng phi tuyến tính theo quy luật đường cong Nukiyama:

Độ quá nhiệt $\Delta T$ ($^\circ\text{C}$) Mật độ dòng nhiệt $q$ ($\text{W/m}^2$) Hệ số tỏa nhiệt $\alpha$ ($\text{W/m}^2\text{.K}$) Chế độ dòng chảy quan sát
5 $4.25 \times 10^3$ $850.0$ Đối lưu tự nhiên đơn pha
10 $1.86 \times 10^4$ $1,860.0$ Khởi phát tâm sôi (ONB)
15 $5.42 \times 10^4$ $3,613.3$ Sôi bọt cục bộ
20 $1.28 \times 10^5$ $6,400.0$ Sôi bọt phát triển mạnh
25 $2.65 \times 10^5$ $10,600.0$ Sôi bọt dày đặc
30 $4.82 \times 10^5$ $16,066.7$ Cận kề điểm sôi màng 
Biểu đồ tương quan hệ số tỏa nhiệt alpha theo độ quá nhiệt:
alpha (W/m2.K)
  16000 |                                                 * (30°C, 16066.7)
  12000 |                                       * (25°C, 10600.0)
   8000 |                             * (20°C, 6400.0)
   4000 |                   * (15°C, 3613.3)
      0 +---------*---------*-----------------------------------> Delta T (°C)
        0         5(850)    10(1860)  15        20        25        30

2. Ảnh hưởng của góc tiếp xúc $\theta$ và sức căng bề mặt $\sigma$

  • Góc tiếp xúc $\theta$ ($75^\circ \to 90^\circ$): Khi góc dính ướt tăng lên $90^\circ$ (bề mặt kém dính ướt hơn), đường kính tách ly bọt hơi $D_t$ tăng từ $1.82\text{ mm}$ lên $2.18\text{ mm}$. Bọt hơi bám trải rộng trên bề mặt làm giảm diện tích tiếp xúc giữa pha lỏng và vách nhiệt, dẫn đến hệ số tỏa nhiệt đối lưu $\alpha$ giảm bình quân $14.6%$.
  • Sức căng bề mặt $\sigma$ ($0.059 \to 0.072\text{ N/m}$): Làm tăng bán kính tới hạn $R_o$, tần số sinh bọt giảm nhẹ nhưng biến thiên thông lượng nhiệt tổng thể không vượt quá $3.8%$.

3. Hiệu quả vượt trội của Nanofluid $\text{SiO}_2$/Nước

Bổ sung các hạt nano $\text{SiO}_2$ kích thước $20\text{ nm}$ với nồng độ thể tích từ $0.01%$ đến $0.06%$ mang lại bước nhảy vọt về hiệu suất trao đổi nhiệt:

  • Tại $\phi = 0.01%$: Hệ số tỏa nhiệt $\alpha$ tăng $8.2%$ so với nước nguyên chất.
  • Tại $\phi = 0.06%$: Hệ số $\alpha$ đạt giá trị cực đại, tăng tới $34.5%$ ở $\Delta T = 25^\circ\text{C}$. Hạt nano lơ lửng làm tăng độ dẫn nhiệt hiệu dụng của môi chất, đồng thời các hạt lắng đọng vi mô tạo thêm vô số tâm hóa hơi nhân tạo trên bề mặt gia nhiệt.

4. So sánh cấu hình 2D Tấm phẳng thép vs 3D Ống trơn đồng

  • Ống trơn nằm ngang 3D bằng vật liệu đồng ($k = 387.6\text{ W/m.K}$) cho hiệu suất truyền nhiệt vượt trội hơn tấm phẳng thép đáy bể ($k = 16.27\text{ W/m.K}$) từ $38%$ đến $52%$ ở cùng mức độ quá nhiệt.
  • Dạng hình học cong giúp các bọt hơi hình thành ở nửa dưới của ống trượt dọc theo chu vi lên đỉnh ống, tạo hiệu ứng khuấy trộn cưỡng bức tự nhiên cực mạnh trong lớp biên nhiệt.

Đổi mới và đóng góp

  1. Ứng dụng thành công mô hình Eulerian giải bài toán sôi bể có hạt nano: Khác với các nghiên cứu trước đây chỉ dùng mô hình Mixture đơn giản hoặc thực nghiệm đơn thuần, đề tài đã tích hợp cơ chế tương tác đa pha Eulerian đầy đủ với các hạt nano $\text{SiO}_2$, giúp mô phỏng chính xác sự tăng cường truyền nhiệt mà không bị phân kỳ số.
  2. Khảo sát đa thông số đồng thời: Định lượng hóa mối quan hệ phi tuyến giữa 5 biến độc lập: Độ quá nhiệt ($\Delta T$), Góc tiếp xúc ($\theta$), Sức căng bề mặt ($\sigma$), Nồng độ thể tích nano ($\phi$), và Cấu hình hình học/vật liệu (Thép vs Đồng; Tấm phẳng vs Ống tròn).
  3. Cải thiện hiệu suất truyền nhiệt định lượng:
    • Sử dụng $\text{SiO}_2$ $0.06%$ vol giúp tăng hệ số $\alpha$ thêm $34.5%$.
    • Chuyển đổi cấu hình sang ống đồng nằm ngang 3D giúp tăng thông lượng nhiệt thêm $45.2%$.
  4. Đối chiếu thực nghiệm với các công bố quốc tế: Kết quả mô phỏng đạt độ tương thích cao với dữ liệu thực nghiệm của Mahian et al. (sai số $< 7.5%$) và Abhishek Kumar et al. (sai số $< 6.2%$).

Ứng dụng thực tế và triển khai

+-----------------------------------------------------------------------------------+
|                        LĨNH VỰC ỨNG DỤNG THỰC TIỄN                                |
+-----------------------------------------------------------------------------------+
|  1. Data Center Cooling: Tản nhiệt ngâm pha lỏng (Immersion Phase-Change Cooling).|
|  2. Điện tử công suất: Module tản nhiệt cho Inverter xe điện (EV) và trạm sạc.     |
|  3. Năng lượng hạt nhân: Thiết kế thanh nhiên liệu và hệ thống làm mát khẩn cấp.  |
|  4. Công nghiệp lạnh & HVAC: Tối ưu dàn bay hơi dạng ngập dịch (Flooded Evaporator)|
+-----------------------------------------------------------------------------------+
  • Mô hình triển khai công nghiệp: Các trung tâm dữ liệu AI/HPC có thể ứng dụng trực tiếp kết quả về góc dính ướt và ống đồng bọc vi cấu trúc để thiết kế các bể làm mát ngâm (Immersion Cooling Tanks) sử dụng chất lỏng điện môi hoặc dung dịch hạt nano loãng.
  • Phân tích chi phí - lợi ích (ROI): Việc tăng $34.5%$ hiệu số truyền nhiệt cho phép giảm $25%$ diện tích bề mặt trao đổi nhiệt của dàn bay hơi công nghiệp, tiết kiệm ước tính $18 - 22%$ chi phí nguyên vật liệu kim loại màu (đồng/nhôm) và giảm $15%$ điện năng tiêu thụ cho hệ thống bơm tuần hoàn.
  • Lộ trình thương mại hóa:
    • Giai đoạn 1 (6 tháng): Chế tạo mẫu thử nghiệm (PoC) ống truyền nhiệt tối ưu hóa biên dạng và lớp phủ ưa nước.
    • Giai đoạn 2 (12 tháng): Thử nghiệm độ bền dài hạn của dung dịch nanofluid chống lắng cặn.
    • Giai đoạn 3 (18 tháng): Tích hợp vào dây chuyền sản xuất thiết bị trao đổi nhiệt công nghiệp.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật:

  • Chưa xét đến hiện tượng kết tụ và lắng đọng lâu ngày của các hạt nano $\text{SiO}_2$ trên bề mặt tạo thành lớp cặn nhiệt trở (fouling resistance).
  • Mô hình số giả định góc tiếp xúc $\theta$ đồng nhất trên toàn bộ bề mặt, chưa mô phỏng độ nhám bất đồng nhất vi mô (micro-roughness).
  • Giới hạn mô phỏng ở điều kiện sôi bể tĩnh, chưa kết hợp vận tốc dòng chảy cưỡng bức ngoài.

Hướng nghiên cứu phát triển:

  • Ứng dụng mô hình bọt hơi rời rạc RPI (Rensselaer Polytechnic Institute) Boiling Model để mô phỏng chính xác hơn tần số rời bọt.
  • Nghiên cứu các chất làm lạnh thế hệ mới có chỉ số GWP thấp như R-1234ze, R-513A kết hợp với ống có cánh vi mô (micro-fin tubes).
  • Áp dụng trí tuệ nhân tạo (Physics-Informed Neural Networks - PINNs) để tăng tốc độ dự đoán thông lượng nhiệt tới hạn CHF theo thời gian thực.

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

  • Sinh viên & Học viên cao học: Cung cấp tài liệu tham khảo chuẩn mực về phương pháp thiết lập mô hình đa pha Eulerian và mô hình rối trong ANSYS Fluent 19.0.
  • Kỹ sư mô phỏng & Thiết kế R&D: Nắm vững bộ thông số biên, phương pháp chia lưới và kỹ thuật xử lý hội tụ cho các bài toán truyền nhiệt đổi pha phức tạp.
  • Doanh nghiệp sản xuất thiết bị nhiệt (HVAC/Boiler): Nhận được các khuyến nghị thiết kế định lượng về vật liệu, hình học ống và xử lý bề mặt nhằm giảm chi phí vật liệu.
  • Các nhà nghiên cứu chuyên sâu: Dữ liệu benchmark tin cậy về cơ chế truyền nhiệt của chất lỏng nano $\text{SiO}_2$ phục vụ các công bố khoa học quốc tế.

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

1. Cấu hình phần cứng tối thiểu để chạy mô phỏng CFD sôi đa pha này là gì?

Hệ thống khuyến nghị trang bị CPU tối thiểu 8 nhân/16 luồng (như Intel Xeon hoặc AMD Ryzen 7 trở lên), 32 GB RAM (tốt nhất là 64 GB cho lưới 3D trên 1 triệu phần tử) và card đồ họa chuyên dụng hỗ trợ tăng tốc tính toán song song.

2. Tại sao lại chọn mô hình Eulerian thay vì mô hình VOF?

Mô hình Eulerian xử lý các pha như các môi trường liên tục tương tác lẫn nhau, cho phép mô phỏng ổn định với khối lượng lớn bọt hơi sinh ra đồng thời và tương thích tự nhiên với việc hòa trộn hạt nano phân tán, trong khi VOF đòi hỏi lưới cực mịn và dễ phân kỳ khi số lượng bọt hơi bùng nổ.

3. Làm thế nào để ngăn chặn hiện tượng hạt nano bị lắng đọng trong ứng dụng thực tế?

Trong vận hành công nghiệp, cần bổ sung các chất hoạt động bề mặt (surfactants như SDS, CTAB) kết hợp điều chỉnh độ pH của dung dịch và duy trì chế độ rung siêu âm định kỳ để giữ các hạt nano $\text{SiO}_2$ luôn ở trạng thái phân tán đồng nhất.

4. Độ chính xác của mô hình CFD này so với thực nghiệm thực tế đạt bao nhiêu phần trăm?

Sai số tương đối giữa kết quả mô phỏng số và dữ liệu thực nghiệm dao động trong khoảng $5.8% \div 8.5%$, nằm hoàn toàn trong ngưỡng sai số cho phép của tính toán kỹ thuật nhiệt ($< 10%$).

5. Việc thay đổi vật liệu từ thép sang đồng cải thiện trao đổi nhiệt như thế nào?

Đồng có hệ số dẫn nhiệt $k \approx 387.6\text{ W/m.K}$, cao hơn gấp 23 lần so với thép ($k \approx 16.27\text{ W/m.K}$). Điều này giúp nhiệt lượng phân bổ cực kỳ đồng đều trên toàn bộ thành ống, triệt tiêu các điểm quá nhiệt cục bộ và thúc đẩy các tâm sôi hoạt động đồng thời.

Kết luận

Đồ án tốt nghiệp "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)" đã hoàn thành xuất sắc các mục tiêu đề ra thông qua việc xây dựng thành công mô hình mô phỏng số CFD trên nền tảng ANSYS Fluent 19.0. Nghiên cứu đã chứng minh một cách khoa học rằng:

  • Việc gia tăng độ dính ướt (giảm góc tiếp xúc $\theta$ từ $90^\circ$ xuống $75^\circ$) giúp cải thiện hệ số tỏa nhiệt đối lưu thêm $\sim 14.6%$.
  • Ứng dụng dung dịch Nanofluid $\text{SiO}_2$/nước ở nồng độ thể tích $0.06%$ mang lại hiệu suất tăng cường truyền nhiệt vượt trội tới $34.5%$.
  • Cấu hình bề mặt dạng ống trơn nằm ngang bằng vật liệu đồng trong không gian 3D tối ưu hóa quá trình giải phóng bọt hơi, nâng cao thông lượng nhiệt lên tới $45.2%$ so với tấm phẳng thép.

Những kết quả này không chỉ đóng góp giá trị học thuật quan trọng cho lĩnh vực nhiệt - thủy động lực học mà còn mở ra giải pháp công nghệ mang tính ứng dụng cao cho việc thiết kế các hệ thống tản nhiệt thế hệ mới trong kỷ nguyên chuyển đổi số và phát triển năng lượng bền vững.