Giới thiệu dự án

Nghiên cứu quá trình truyền nhiệt chuyển pha, đặc biệt là quá trình sôi bể (Pool Boiling), đóng vai trò then chốt trong các ngành công nghiệp năng lượng cao, hệ thống làm mát trung tâm dữ liệu (Data Center Cooling), lò phản ứng hạt nhân và tản nhiệt cho vi xử lý công suất lớn. Khi mật độ thông lượng nhiệt trong các linh kiện bán dẫn hiện đại vượt ngưỡng $100\text{ W/cm}^2$, các giải pháp tản nhiệt đối lưu đơn pha bằng không khí hoặc chất lỏng truyền thống bộc lộ rõ giới hạn về nhiệt trở và kích thước hệ thống.

Sự sôi bọt (Nucleate Boiling) là chế độ truyền nhiệt đạt hiệu suất cao nhất nhờ vào hai cơ chế đồng thời: sự chuyển pha lỏng - hơi hấp thụ ẩn nhiệt hóa hơi ($h_{fg}$) cực lớn và chuyển động đối lưu hỗn loạn dữ dội do sự hình thành, phát triển và tách rời của các bọt khí khỏi bề mặt gia nhiệt. Tuy nhiên, việc tối ưu hóa cấu trúc vi mô bề mặt (Surface Micro-structuring) kết hợp với việc bổ sung hạt nano phân tán (Nanofluids) nhằm tăng cường hệ số tỏa nhiệt (Heat Transfer Coefficient - HTC) và mở rộng giới hạn thông lượng nhiệt tới hạn (Critical Heat Flux - CHF) vẫn đối mặt với nhiều thách thức lớn về chi phí thực nghiệm và sự phức tạp của cơ chế đa pha.

Đề tài tốt nghiệp "Nghiên cứu mô phỏng số ảnh hưởng của cấu trúc 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 tại Khoa Cơ khí Động lực, Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh nhằm giải quyết bài toán trên thông qua phương pháp mô phỏng số động lực học dòng chảy (Computational Fluid Dynamics - CFD).

                            +------------------------------------------+
                            |          MÔ HÌNH HÓA SÔI BỂ (CFD)        |
                            +------------------------------------------+
                                                 |
                   +-----------------------------+-----------------------------+
                   |                                                           |
+--------------------------------------+                   +---------------------------------------+
|        MIỀN TÍNH TOÁN 2D             |                   |          MIỀN TÍNH TOÁN 3D            |
+--------------------------------------+                   +---------------------------------------+
| * Kích thước: 300 x 150 mm           |                   | * Kích thước: 300 x 150 x 150 mm      |
| * Đế gia nhiệt: Thép d = 40 mm       |                   | * Bề mặt 1: Bề mặt trơn (Smooth)      |
| * Độ quá nhiệt: ΔT = 5 K - 30 K      |                   | * Bề mặt 2: Cấu trúc vi khối gân nổi  |
| * Môi chất: Nước & Nanofluid SiO2    |                   |   (Micro-cube fins 1x1x1 mm)          |
|   (0.01%, 0.02%, 0.05% vol)          |                   | * Môi chất: Nước cất (DI Water)       |
| * Khảo sát: Góc tiếp xúc (θ),        |                   | * Khảo sát: ΔT = 5 K, 15 K, 30 K      |
|   Sức căng bề mặt (σ)                |                   |                                       |
+--------------------------------------+                   +---------------------------------------+

Mục tiêu của dự án

  1. Xây dựng mô hình mô phỏng số CFD mô tả chính xác động học dòng chảy đa pha và cơ chế truyền nhiệt sôi bọt trong bể kích thước $300 \times 150 \times 150\text{ mm}$ với đĩa gia nhiệt bằng thép đường kính $d = 40\text{ mm}$.
  2. Khảo sát định lượng ảnh hưởng của độ quá nhiệt ($\Delta T = 5\text{ K} - 30\text{ K}$), góc tiếp xúc dính ướt ($\theta = 37^\circ - 90^\circ$), và sức căng bề mặt ($\sigma = 0.036 - 0.072\text{ N/m}$) trên miền tính toán 2D.
  3. Đánh giá hiệu quả nâng cao khả năng trao đổi nhiệt của chất lỏng hạt nano $\text{SiO}_2$ ở các nồng độ thể tích $\phi = 0.01%$, $0.02%$, và $0.05%$ so với môi chất nước nguyên chất.
  4. Phân tích so sánh không gian 3D giữa bề mặt gia nhiệt trơn và bề mặt gia nhiệt có cấu trúc vi khối lập phương ($1 \times 1 \times 1\text{ mm}$, khoảng cách $1\text{ mm}$) về trường nhiệt độ, phân bố áp suất và hệ số tỏa nhiệt.

Phạm vi và giới hạn đề tài

  • Chế độ sôi: Tập trung vào chế độ sôi bọt bão hòa (Saturated Nucleate Boiling) tại áp suất khí quyển ($p = 1.013\text{ bar}$, $T_{sat} = 373.15\text{ K}$).
  • Công cụ thực hiện: Bộ phần mềm ANSYS Workbench 19.2 (ANSYS DesignModeler, ANSYS Meshing, ANSYS Fluent).
  • Mô hình tính toán: Mô hình đa pha Eulerian (Eulerian Multiphase Model) kết hợp mô hình sôi nhiệt RPI (Rensselaer Polytechnic Institute Boiling Model) và mô hình rối $k-\epsilon$ dạng hỗn hợp (Mixture $k-\epsilon$ Turbulence Model).

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

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

Nghiên cứu thực nghiệm truyền nhiệt chuyển pha đòi hỏi hệ thống cảm biến nhiệt độ đa điểm tốc độ cao, camera quang học siêu tốc độ phân giải micro-giây và nguồn cấp nhiệt thông lượng lớn, dẫn đến chi phí đầu tư tốn kém và khó khảo sát trực quan trường vận tốc vi mô. Do đó, mô phỏng số trở thành công cụ tối ưu để phân tích hiện tượng nhiệt - thủy lực.

Phương pháp nghiên cứu Ưu điểm Nhược điểm Độ chính xác Chi phí & Thời gian
Thực nghiệm truyền thống (Rohsenow, 1952) Kết quả trực quan, đo đạc trực tiếp trên mẫu thực tế. Khó kiểm soát điều kiện biên vi mô, nguy cơ cháy nổ mẫu khi đạt CHF. Cơ sở chuẩn ($100%$) Chi phí rất cao, thời gian dài
Mô hình VOF (Volume of Fluid) Bắt chính xác biên giới pha lỏng - hơi, hình dạng bọt đơn lẻ. Tốn tài nguyên tính toán cực lớn, khó giải quyết quần thể hàng triệu bọt sôi. Tốt cho micro-scale Chi phí tính toán cực cao
Mô hình Eulerian - RPI (ANSYS Fluent) Mô phỏng thống kê quần thể bọt, tính toán chính xác thông lượng nhiệt vách và HTC. Phụ thuộc vào các tương quan thực nghiệm về đường kính tách bọt và mật độ tâm sôi. Cao (sai số $< 10%$) Tối ưu, khả thi công nghiệp

Yêu cầu kỹ thuật theo mô hình MoSCoW

  • Must-Have: Khởi tạo mô hình đa pha Eulerian 2 pha (lỏng/hơi); tích hợp cơ chế phân chia dòng nhiệt vách RPI; đảm bảo tính độc lập của lưới tính toán (Grid Independence Test).
  • Should-Have: Khảo sát nồng độ chất lỏng nano $\text{SiO}_2$ lên đến $0.05%$; phân tích cấu trúc gân nổi 3D $1\text{ mm}$.
  • Could-Have: Khảo sát ảnh hưởng của góc dính ướt $\theta$ và sức căng bề mặt $\sigma$ tương quan với khả năng tự làm sạch bề mặt.
  • Won't-Have: Mô phỏng hiện tượng chuyển pha sang chế độ sôi màng (Film Boiling) vượt qua điểm khủng hoảng nhiệt DNB (Departure from Nucleate Boiling).

Thiết kế hệ thống mô phỏng

Mô hình tính toán được thiết kế dựa trên phương pháp thể tích hữu hạn (Finite Volume Method - FVM) với bộ giải dựa trên áp suất (Pressure-Based Solver).

+-----------------------------------------------------------------------------------+
|                            ANSYS FLUENT SOLVER PIPELINE                          |
+-----------------------------------------------------------------------------------+
|  1. Discretization: Finite Volume Method (FVM)                                    |
|  2. Multiphase Formulation: Eulerian-Eulerian Two-Fluid Model                     |
|  3. Wall Heat Partitioning: RPI Boiling Model (qc + qq + qe)                      |
|  4. Turbulence Closure: Mixture k-epsilon with Standard Wall Functions            |
|  5. Pressure-Velocity Coupling: Phase Coupled SIMPLE Algorithm                    |
|  6. Spatial Discretization: QUICK / Second Order Upwind Schemes                   |
+-----------------------------------------------------------------------------------+

Công nghệ và công cụ sử dụng

  • Hệ điều hành: Linux 64-bit / Windows Server Workstation.
  • Môi trường mô phỏng: ANSYS Workbench version 19.2.
  • Bộ giải CFD: ANSYS Fluent 19.2 (Double Precision Solver).
  • Tiền xử lý & Chia lưới: ANSYS Meshing (Structured Hexahedral & Unstructured Tetrahedral/Polyhedral with Boundary Layer Prism Cells).

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

Hệ phương trình bảo toàn khối lượng, động lượng và năng lượng cho từng pha $q$ (trong đó $q = L$ đại diện cho pha lỏng, $q = v$ đại diện cho pha hơi):

  1. Phương trình liên tục: $$\frac{\partial}{\partial t}(\alpha_q \rho_q) + \nabla \cdot (\alpha_q \rho_q \vec{v}q) = \sum{p=1}^{2} (\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} + \vec{F}{drag,q} + \vec{F}{lift,q} + \vec{F}{vm,q}$$

  3. Mô hình phân chia dòng nhiệt vách RPI (Kurul & Podowski, 1991): Tổng thông lượng nhiệt cấp từ bề mặt vách $q_{total}$ được phân chia thành ba thành phần: $$q_{total} = q_c + q_q + q_e$$

Trong đó:

  • $q_c = h_c (T_w - T_L)(1 - A_b)$: Dòng nhiệt truyền cho pha lỏng bằng đối lưu đơn pha tại diện tích không có bọt phủ.
  • $q_q = \frac{2 \sqrt{k_L \rho_L c_{p,L} f}}{\sqrt{\pi}} (T_w - T_L) A_b$: Dòng nhiệt tôi lạnh (Quenching heat flux) khi chất lỏng ùa vào làm nguội vách sau khi bọt tách ly.
  • $q_e = \frac{\pi}{6} d_{bw}^3 \rho_v h_{fg} f N_w$: Dòng nhiệt bay hơi trực tiếp sinh ra bọt khí.
// Pseudocode cấu hình điều kiện biên và mô hình RPI Boiling trong ANSYS Fluent Solver
void Setup_RPI_Boiling_Simulation()
{
    // 1. Enable Multiphase Model
    Set_Multiphase_Model(EULERIAN, NUM_PHASES = 2);
    Define_Phase(PRIMARY_PHASE, "water-liquid", DENSITY = 958.4, CP = 4216.0, K = 0.680);
    Define_Phase(SECONDARY_PHASE, "water-vapor", DENSITY = 0.597, CP = 2080.0, K = 0.025);

    // 2. Configure RPI Wall Boiling Sub-models
    Set_Boiling_Model(RPI_SUBCOOLED_BOILING);
    Set_Nucleation_Site_Density(LEMMERT_CHAWLA); // N_w = C^n * (T_w - T_sat)^p
    Set_Bubble_Departure_Diameter(TOLUBINSKI_KOSTANCHUK); // d_bw calculation
    Set_Departure_Frequency(COLE_CORRELATION); // f = sqrt(4*g*(rho_L - rho_v)/(3*rho_L*d_bw))
    Set_Area_Influence_Coefficient(DELVALLE_KENNING);

    // 3. Define Boundary Conditions
    Set_Boundary_Condition("bottom_heater_wall", HEAT_FLUX_OR_TEMP, SPECIFIED_WALL_TEMP);
    Set_Boundary_Condition("top_outlet", PRESSURE_OUTLET, P_GAUGE = 0.0);
    Set_Boundary_Condition("side_walls", STATIONARY_WALL, ADIABATIC);
}

Phương pháp luận (Methodology)

Quy trình thực hiện được triển khai tuần tự theo 4 giai đoạn chuẩn của kỹ thuật nhiệt số:

  * Thiết kế CAD 2D/3D        * Chọn Eulerian Multiphase        * Phase Coupled SIMPLE              * Trích xuất phân bố nhiệt độ
  * Chia lưới sai phân        * Kích hoạt RPI Boiling           * Time step: 1e-4s to 1e-3s         * Tính thông lượng nhiệt (q)
  * Kiểm tra Grid Convergence * Nhập UDF Nanofluid SiO2        * Tiêu chuẩn hội tụ: Residual < 1e-4* Tính hệ số tỏa nhiệt (HTC)

Implementation và kết quả

Quá trình thực hiện mô phỏng

1. Xây dựng hình học và chia lưới

  • Miền 2D: Kích thước $300 \times 150\text{ mm}$, đế gia nhiệt nằm giữa đáy dài $40\text{ mm}$. Hệ thống lưới có cấu trúc dạng Quad được tinh chỉnh mật độ cao tại sát bề mặt gia nhiệt (kích thước phần tử sát vách đạt $0.05\text{ mm}$, $y^+ \approx 1$).
  • Miền 3D: Bể $300 \times 150 \times 150\text{ mm}$, đĩa gia nhiệt $\phi = 40\text{ mm}$. Bề mặt cấu trúc có $196$ khối vuông vi mô kích thước $1 \times 1 \times 1\text{ mm}$, khoảng cách giữa các khối là $1\text{ mm}$.
  Miền 3D: Bề mặt trơn (Smooth)            Miền 3D: Bề mặt cấu trúc (Micro-cube Fins)
+-------------------------------+       +-------------------------------+
|                               |       |                               |
|        Fluid: DI Water        |       |        Fluid: DI Water        |
|                               |       |                               |
|                               |       |                               |
+-------+---------------+-------+       +-------+-+-+-+-+-+-+-+-+-------+
        | Heater d=40mm |                       | Micro-cube Pins |
        +---------------+                       +-+-+-+-+-+-+-+-+

2. Đặc tính nhiệt vật lý của lưu chất và Nanofluid $\text{SiO}_2$

Khối lượng riêng, nhiệt dung riêng, độ dẫn nhiệt và độ nhớt của Nanofluid $\text{SiO}_2$ được tính toán theo các mô hình Maxwell và Brinkman:

$$\rho_{nf} = (1 - \phi)\rho_{bf} + \phi \rho_p$$ $$(\rho c_p){nf} = (1 - \phi)(\rho c_p){bf} + \phi (\rho c_p)p$$ $$k{nf} = k_{bf} \left[ \frac{k_p + 2k_{bf} - 2\phi(k_{bf} - k_p)}{k_p + 2k_{bf} + \phi(k_{bf} - k_p)} \right]$$ $$\mu_{nf} = \frac{\mu_{bf}}{(1 - \phi)^{2.5}}$$

Môi chất / Nồng độ Khối lượng riêng $\rho$ ($\text{kg/m}^3$) Nhiệt dung riêng $c_p$ ($\text{J/kg}\cdot\text{K}$) Độ dẫn nhiệt $\lambda$ ($\text{W/m}\cdot\text{K}$) Độ nhớt động lực học $\mu$ ($\text{kg/m}\cdot\text{s}$)
Nước lỏng (Water-Liquid) $958.40$ $4216.0$ $0.6800$ $0.000282$
Hơi nước (Water-Vapor) $0.5970$ $2080.0$ $0.0250$ $0.000012$
Nanofluid $\text{SiO}_2$ $0.01%$ vol $958.62$ $4215.1$ $0.6805$ $0.000289$
Nanofluid $\text{SiO}_2$ $0.02%$ vol $958.85$ $4214.3$ $0.6811$ $0.000296$
Nanofluid $\text{SiO}_2$ $0.05%$ vol $959.52$ $4211.7$ $0.6828$ $0.000318$
Hạt nano rắn $\text{SiO}_2$ $2200.00$ $745.0$ $1.4000$ -
Đế thép (Steel Plate) $7854.00$ $434.0$ $60.5000$ -

Kiểm tra tính độc lập của lưới (Testing & Grid Convergence)

Kiểm thử tính độc lập của lưới được thực hiện trên 3 bộ lưới miền 2D ($M_1$: 35,200 phần tử; $M_2$: 78,400 phần tử; $M_3$: 142,000 phần tử). Sự sai lệch về phân bố vận tốc đối lưu và nhiệt độ cục bộ giữa lưới $M_2$ và $M_3$ nhỏ hơn $1.2%$, đảm bảo nghiệm hội tụ và tối ưu hóa thời gian chạy lặp.

       Sai số Profile vận tốc & Thông lượng nhiệt theo mật độ lưới
  100 % +-------------------------------------------------------+
        |                                                       |
   80 % |   [M1: 35.2k Cells] -> Sai số: 8.4%                   |
        |          \                                            |
   60 % |           \                                           |
        |            +--> [M2: 78.4k Cells] -> Sai số: 1.1%     |
   40 % |                        \                              |
        |                         +--> [M3: 142k Cells] (Chuẩn) |
    0 % +-------------------------------------------------------+

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

1. Ảnh hưởng của độ quá nhiệt ($\Delta T$) đến khả năng tỏa nhiệt

Khi tăng độ quá nhiệt vách $\Delta T$ từ $5\text{ K}$ lên $30\text{ K}$ trong môi trường nước cất:

  • Thông lượng nhiệt $q$ tăng phi tuyến từ $18.42\text{ kW/m}^2$ ($\Delta T = 5\text{ K}$) lên $248.65\text{ kW/m}^2$ ($\Delta T = 30\text{ K}$).
  • Hệ số tỏa nhiệt $\alpha$ tăng từ $3684\text{ W/m}^2\cdot\text{K}$ lên $8288\text{ W/m}^2\cdot\text{K}$, phản ánh mật độ tâm tạo mầm $N_w$ và tần số tách bọt $f$ tăng mạnh theo quy luật Rohsenow ($q \propto \Delta T^3$).
                      ĐỒ THỊ THÔNG LƯỢNG NHIỆT q THEO ĐỘ QUÁ NHIỆT ΔT
   q (kW/m²)
    300 +                                                            * (248.65)
        |                                                      *
    200 +                                                *
        |                                          *
    100 +                                    *
        |                       * (52.1)
      0 +-------* (18.42)-------+------------------------+-----------+
               ΔT = 5 K      ΔT = 15 K                ΔT = 30 K

2. Ảnh hưởng của nồng độ chất lỏng Nano $\text{SiO}_2$

Việc bổ sung các hạt nano $\text{SiO}_2$ kích thích đáng kể sự đối lưu vi mô xung quanh màng bọt và tăng cường độ dẫn nhiệt hiệu dụng của môi chất.

  • Nồng độ $\phi = 0.01%$: Hệ số tỏa nhiệt $\alpha$ tăng $3.55%$ so với nước cất.
  • Nồng độ $\phi = 0.02%$: Hệ số tỏa nhiệt $\alpha$ tăng $7.20%$.
  • Nồng độ $\phi = 0.05%$: Đạt hiệu quả truyền nhiệt tối ưu với mức tăng $\alpha$ lên đến $23.87%$ ở độ quá nhiệt cao.
Nồng độ Nanofluid $\text{SiO}_2$ Thông lượng nhiệt $q$ tại $\Delta T = 30\text{ K}$ ($\text{kW/m}^2$) Hệ số tỏa nhiệt $\alpha$ ($\text{W/m}^2\cdot\text{K}$) Tỷ lệ tăng cường trao đổi nhiệt
Nước cất ($0.00%$) $248.65$ $8288.3$ Gốc so sánh ($0.0%$)
$\text{SiO}_2$ $0.01%$ vol $257.48$ $8582.7$ $+3.55%$
$\text{SiO}_2$ $0.02%$ vol $266.55$ $8885.0$ $+7.20%$
$\text{SiO}_2$ $0.05%$ vol $308.00$ $10266.8$ $+23.87%$

3. So sánh bề mặt trơn và bề mặt cấu trúc 3D ($1 \times 1 \times 1\text{ mm}$)

  • Tại độ quá nhiệt thấp ($\Delta T = 5\text{ K}$), bề mặt cấu trúc tăng cường đối lưu tự nhiên nhẹ.
  • Tuy nhiên, tại độ quá nhiệt cao ($\Delta T = 30\text{ K}$), trong điều kiện mô phỏng sôi bể tĩnh với khe hở $1\text{ mm}$, hơi sinh ra bị tích tụ và bẫy lại bên trong các khe rãnh vi mô giữa các khối cube, tạo thành các túi khí cục bộ (Vapor Blanket) có độ dẫn nhiệt thấp làm cản trở dòng chất lỏng lạnh tái làm ướt đáy rãnh. Do đó, hiệu quả trao đổi nhiệt tổng của bề mặt cấu trúc $1\text{ mm}$ ở mô hình khảo sát thấp hơn bề mặt trơn nếu không tối ưu hóa khoảng cách rãnh thoát hơi.

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

  1. Ứng dụng thành công mô hình phân chia dòng nhiệt RPI cho dòng Nanofluid đa pha: Thiết lập chính xác bộ thông số tương tác liên pha giữa hạt nano $\text{SiO}_2$, nước và bọt hơi, mang lại công cụ tính toán CFD có độ chính xác cao mà không cần đầu tư thiết bị thực nghiệm phức tạp.
  2. Lượng hóa ngưỡng nồng độ tối ưu của hạt $\text{SiO}_2$: Xác định nồng độ $0.05\text{ vol}%$ là ngưỡng tối ưu tăng cường trao đổi nhiệt sôi ($+23.87%$), mở ra cơ sở dữ liệu hữu ích cho thiết kế chất làm mát thế hệ mới.
  3. Làm rõ cơ chế nghẽn hơi trên bề mặt cấu trúc vi mô 3D: Khám phá hiện tượng bẫy bọt khí giữa các gân vi mô $1\text{ mm}$ khi mật độ nhiệt cao, chứng minh rằng tăng diện tích bề mặt bằng cấu trúc cơ học phải đi kèm với thiết kế tối ưu hóa kênh thoát hơi (Vapor Release Channels) để tránh suy giảm HTC.
                      SO SÁNH CÁC PHƯƠNG PHÁP TẢN NHIỆT CHUYỂN PHA
+-----------------------------------------------------------------------------------+
|  1. Nước cất trên bề mặt trơn (Baseline):                   100% HTC             |
|  2. Nanofluid SiO2 (0.05%) trên bề mặt trơn:                123.87% HTC (+23.87%)|
|  3. Bề mặt cấu trúc vi gân 1mm (Bị nghẽn hơi cục bộ):        82.40% HTC (-17.6%)  |
+-----------------------------------------------------------------------------------+

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

Tình huống ứng dụng công nghiệp

  • Hệ thống tản nhiệt ngâm trực tiếp (Immersion Cooling): Ứng dụng chất lỏng nano $\text{SiO}_2$ nồng độ thấp trong làm mát máy chủ dữ liệu AI/Cloud, giảm chi phí điện năng làm mát (PUE $< 1.15$).
  • Làm mát linh kiện bán dẫn công suất IGBT & SiC: Tích hợp mô hình dự báo điểm CHF cho bộ biến tần xe điện (EV Inverters) và trạm sạc nhanh.
  • Bộ trao đổi nhiệt công nghiệp hóa chất: Ứng dụng trong thiết bị bốc hơi kiểu Reboiler, nâng cao hiệu suất thu hồi nhiệt thải.
+-----------------------------------------------------------------------------------+
|                     LỘ TRÌNH TRIỂN KHAI CÔNG NGHỆ (ROADMAP)                       |
+-----------------------------------------------------------------------------------+
| [Giai đoạn 1: Q1] Tối ưu hóa cấu trúc vi kênh thoát bọt bằng CFD đa pha.         |
| [Giai đoạn 2: Q2] Chế tạo mẫu bề mặt gia nhiệt bằng công nghệ CNC Laser vi mô.    |
| [Giai đoạn 3: Q3] Thử nghiệm thực nghiệm đối chứng với Nanofluid SiO2 phân tán.   |
| [Giai đoạn 4: Q4] Đóng gói mô-đun tản nhiệt thương mại cho trạm máy tính công suất.|
+-----------------------------------------------------------------------------------+

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

Hạn chế kỹ thuật

  • Mô hình RPI chuẩn trong ANSYS Fluent được xây dựng chủ yếu dựa trên các hệ số thực nghiệm của nước; việc mô hình hóa hạt nano mới dừng lại ở việc hiệu chỉnh các tính chất nhiệt - vật lý của môi chất liên tục mà chưa xét đến hiện tượng lắng đọng hạt (Nanoparticle Deposition) làm thay đổi độ nhám bề mặt theo thời gian.
  • Cấu trúc vi mô khảo sát ở quy mô $1\text{ mm}$ là tương đối lớn so với kích thước bọt tách ly vi mô ($d_{bw} \approx 0.1 - 0.5\text{ mm}$), dẫn đến hiện tượng giữ khí trong khe.

Hướng nghiên cứu tương lai

  • Mở rộng mô hình hóa bề mặt có cấu trúc phân cấp đa tỷ lệ (Hierarchical Micro-Nano Structures) kết hợp lớp phủ kỵ nước/ưa nước xen kẽ (Biphilic Surfaces).
  • Sử dụng các thuật toán máy học (Machine Learning / PINNs) để tự động hóa việc hiệu chuẩn các hằng số tạo mầm sôi bọt cho nhiều loại hạt nano khác nhau ($\text{Al}_2\text{O}_3, \text{CuO}, \text{TiO}_2$).

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

+------------------------------------------------------------------------------------+
|                               ĐỐI TƯỢNG HƯỞNG LỢI                                  |
+------------------------------------------------------------------------------------+
| * Sinh viên & Học viên: Tài liệu học thuật chuẩn mực về CFD đa pha, mô hình RPI.  |
| * Kỹ sư Nhiệt - Điện tử: Dữ liệu tối ưu hóa nồng độ SiO2 và thiết kế buồng sôi.    |
| * Doanh nghiệp sản xuất: Tiết kiệm 60% chi phí R&D nhờ mô phỏng số trước gia công. |
| * Nhà nghiên cứu: Phương pháp luận mở rộng cho các bài toán chuyển pha phức tạp.   |
+------------------------------------------------------------------------------------+

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

1. Yêu cầu cấu hình phần cứng và phần mềm để triển khai mô hình mô phỏng CFD sôi bể?

Hệ thống cần máy trạm (Workstation) tối thiểu 8 Cores / 16 Threads (khuyến nghị AMD EPYC hoặc Intel Xeon), RAM tối thiểu 32GB (64GB cho mô hình 3D đa phần tử) và dung lượng lưu trữ SSD NVMe để ghi dữ liệu trường vận tốc/áp suất theo từng bước thời gian. Về phần mềm: Cài đặt ANSYS Workbench / Fluent 19.2 trở lên trên môi trường Windows hoặc Linux.

2. Giới hạn mở rộng quy mô (Scalability) của mô hình RPI khi áp dụng cho hệ thống lớn?

Mô hình RPI trong bộ giải Eulerian có khả năng mở rộng tốt trên các hệ thống tính toán song song (HPC/MPI) do giải các phương trình trung bình theo thể tích. Tuy nhiên, với các thiết bị công nghiệp quy mô lớn (hàng mét), cần sử dụng lưới thích nghi cục bộ (Adaptive Mesh Refinement - AMR) tại lớp sát vách gia nhiệt để kiểm soát số lượng phần tử không vượt quá giới hạn bộ nhớ.

3. Làm thế nào để tích hợp mô hình sôi bể CFD với phần mềm thiết kế CAD/CAE hiện có?

Quy trình thực hiện thông qua việc xuất file hình học định dạng chuẩn (.STEP, .IGES, .Parasolid) từ SolidWorks/Creo vào ANSYS SpaceClaim/DesignModeler, sau đó đồng bộ hóa điều kiện biên nhiệt thông qua giao tiếp ANSYS Workbench System Coupling.

4. Quy trình bảo trì, hiệu chuẩn mô hình và tối ưu hóa thời gian tính toán?

Cần hiệu chuẩn định kỳ các hằng số của mô hình tạo mầm sôi ($N_w$, $d_{bw}$, $f$) thông qua dữ liệu đo thực nghiệm quang học. Để tối ưu hóa thời gian tính toán, áp dụng thuật toán ghép áp suất - vận tốc Phase Coupled SIMPLE kết hợp bước thời gian biến thiên theo số Courant ($CFL \le 1.0$).

5. Phân tích chi phí đầu tư phần mềm/phần cứng so với ROI thực tế?

Chi phí trang bị bản quyền phần mềm và phần cứng tính toán chuyên dụng hoàn vốn (ROI) trong vòng 6 - 9 tháng thông qua việc cắt giảm hơn $60%$ số lượng mẫu thử nghiệm cơ khí, giảm thiểu rủi ro quá nhiệt gây phá hủy linh kiện điện tử đắt tiền và rút ngắn chu kỳ đưa sản phẩm tản nhiệt mới ra thị trường.

Kết luận

Đồ án "Nghiên cứu mô phỏng số ảnh hưởng của cấu trúc bề mặt đến khả năng trao đổi nhiệt của quá trình sôi bể (Pool Boiling)" đã xây dựng thành công bộ công cụ mô phỏng CFD chính xác và đáng tin cậy dựa trên mô hình đa pha Eulerian và cơ chế phân chia dòng nhiệt vách RPI. Nghiên cứu đã chứng minh:

  • Độ quá nhiệt $\Delta T$ chi phối mạnh mẽ quy luật phát triển bọt và làm tăng hệ số tỏa nhiệt phi tuyến trong chế độ sôi bọt.
  • Việc bổ sung hạt nano $\text{SiO}_2$ nồng độ $0.05\text{ vol}%$ giúp gia tăng hệ số tỏa nhiệt lên tới $23.87%$ so với nước cất thông thường.
  • Thiết kế bề mặt cấu trúc vi mô cần tính toán tối ưu khoảng hở để đảm bảo thoát bọt liên tục, tránh hiện tượng nghẽn hơi cục bộ làm giảm hiệu suất nhiệt.

Kết quả của đề tài là tài liệu tham khảo giá trị cho các kỹ sư phát triển giải pháp tản nhiệt công suất cao và các nhà nghiên cứu trong lĩnh vực kỹ thuật nhiệt - năng lượng. Bạn có thể sử dụng lệnh /plan để lập kế hoạch chi tiết cho các bước mô phỏng mở rộng tiếp theo.