Giới thiệu dự án

Hệ thống trao đổi nhiệt là trái tim của mọi tổ hợp kỹ thuật nhiệt, từ chu trình lạnh công nghiệp, hệ thống điều hòa không khí trung tâm (HVAC) đến các trạm phát điện và cơ cấu động lực trên tàu biển. Theo các báo cáo năng lượng quốc tế, các thiết bị trao đổi nhiệt chiếm tới 20–30% tổng chi phí đầu tư thiết bị và ảnh hưởng trực tiếp đến 40% hiệu suất tiêu thụ năng lượng của toàn hệ thống công nghiệp. Trong môi trường hàng hải, thiết bị trao đổi nhiệt dạng ống có cánh (Finned-Tube Heat Exchanger - FTHE) đối mặt với thách thức khắc nghiệt: không gian lắp đặt hạn chế, độ mặn cao gây ăn mòn, sự biến động liên tục của tải nhiệt và yêu cầu khắt khe về giảm thiểu tổn thất áp suất để tiết kiệm công suất quạt/bơm.

+-------------------------------------------------------------------------+
|                    BÀI TOÁN TỐI ƯU HÓA THIẾT BỊ FTHE                    |
+------------------------------------+------------------------------------+
|            HIỆN TRẠNG              |              GIẢI PHÁP             |
| - Vùng quẩn khí sau ống (Dead Zone)| - Mô phỏng số 3D Navier-Stokes     |
| - Trở lực không khí lớn (ΔP cao)   | - Tối ưu thông số (u, Tin, d)      |
| - Hệ số truyền nhiệt đối lưu thấp  | - Tích hợp cánh phụ tạo xoáy (VG)  |
+------------------------------------+------------------------------------+

Vấn đề kỹ thuật (Problem Statement)

Thiết bị trao đổi nhiệt dạng ống có cánh phẳng truyền thống bộc lộ những giới hạn nhiệt - thủy lực cố hữu:

  • Vùng chết khí động (Dead Zone / Wake Region): Dòng khí chuyển động qua cụm ống tạo ra vùng phân ly dòng và hoàn lưu xoáy phía sau lưng ống, làm giảm đáng kể hệ số truyền nhiệt cục bộ ($h_{local}$).
  • Trở nhiệt lớn phía không khí: Không khí có hệ số dẫn nhiệt rất thấp ($\lambda \approx 0.026 \text{ W/m}\cdot\text{K}$), dẫn đến trở nhiệt phía không khí chiếm hơn 70–80% tổng trở nhiệt của thiết bị.
  • Mất cân đối giữa tăng cường nhiệt và tổn thất áp suất: Việc gia tăng diện tích truyền nhiệt đơn thuần bằng cách tăng mật độ cánh sẽ làm gia tăng tổn thất áp suất ($\Delta P$), đòi hỏi công suất quạt lớn hơn và gây hao tổn năng lượng tổng thể.

Mục tiêu dự án

  1. Xây dựng mô hình tính toán 3D CFD hoàn chỉnh: Mô phỏng chính xác trường nhiệt độ, trường áp suất và trường vận tốc của khối không khí chuyển động qua cụm ống cánh tản nhiệt trong môi trường biển bằng phần mềm chuyên dụng.
  2. Khảo sát tham số ảnh hưởng đa chiều: Đánh giá định lượng tác động của vận tốc đầu vào ($u_{in}$), nhiệt độ dòng khí ($T_{in}$) và đường kính hình học của ống ($d$) đến hệ số truyền nhiệt đối lưu và tổn thất áp suất.
  3. Tối ưu hóa thông số vận hành: Xác định bộ thông số hoạt động tối ưu nhằm nâng cao tỷ số hiệu suất nhiệt - thủy lực ($JF$).
  4. Cải tiến cấu hình hình học cánh: Đề xuất và mô phỏng cấu trúc cánh phụ tạo xoáy vi mô (Micro-Vortex Generators / Sub-fins) để phá vỡ lớp biên nhiệt, nâng cao khả năng trao đổi nhiệt mà không làm tăng đột biến tổn thất áp suất.

Phương pháp tiếp cận và kết quả kỳ vọng

  • Phương pháp tiếp cận: Sử dụng phương pháp thể tích hữu hạn (Finite Volume Method - FVM) trên nền tảng ANSYS Fluent 19.2 để giải hệ phương trình Navier-Stokes 3D và phương trình truyền năng lượng ở trạng thái dừng, kết hợp xác thực độc lập lưới (Grid Independence Test) và kiểm chứng với dữ liệu thực nghiệm chuẩn quốc tế.
  • Kết quả định lượng kỳ vọng: Nâng cao hệ số Colburn $j$, kiểm soát hệ số ma sát Fanning $f$, cải thiện hệ số đánh giá hiệu suất nhiệt - thủy lực $JF > 1.0$, giảm thiểu tối đa kích thước và khối lượng thiết bị trao đổi nhiệt lắp đặt trên tàu biển.

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

  • Phạm vi: Khảo sát mô hình dòng chảy qua cụm ống so le có cánh phẳng trong dải vận tốc $u = 1.0 - 3.0 \text{ m/s}$, đường kính ống $d = 7 - 10 \text{ mm}$, nhiệt độ đầu vào $T_{in} = 273 - 300 \text{ K}$.
  • Giới hạn: Dòng khí được giả định là chất lưu Newton không nén được, tính chất vật lý của không khí đồng nhất ở trạng thái dừng (steady-state), chưa xét đến quá trình đóng cặn ẩm (fouling) và hiện tượng ngưng tụ ẩm trên bề mặt cánh.

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

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

Dạng thiết bị trao đổi nhiệt Ưu điểm chính Nhược điểm kỹ thuật Mức độ phù hợp môi trường biển
Ống cánh phẳng truyền thống (Plain FTHE) Kết cấu đơn giản, chi phí gia công thấp, độ bền cơ học cao Hệ số truyền nhiệt thấp, vùng chết sau ống lớn Trung bình
Cánh lượn sóng (Wavy Fin FTHE) Phá vỡ lớp biên liên tục, tăng 15–25% truyền nhiệt Tổn thất áp suất tăng 20–35%, dễ bám bẩn màng muối biển Khá
Cánh có khe louver (Louvered Fin) Hiệu suất nhiệt rất cao ở dải Re thấp Cấu trúc mỏng manh, dễ ăn mòn rỗ, khó vệ sinh cơ học Kém
Cánh tích hợp cánh phụ tạo xoáy (VG FTHE - Đề xuất) Tạo dòng xoáy thứ cấp cuốn nhiệt, kiểm soát $\Delta P$ tối ưu Đòi hỏi mô phỏng chính xác góc nghiêng và vị trí bố trí Xuất sắc

Ma trận ưu tiên yêu cầu kỹ thuật (MoSCoW)

  • Must-have: Mô hình 3D đối xứng chu kỳ, mô hình dòng chảy giải quyết lớp biên nhiệt chính xác, kiểm tra hội tụ lưới nghiêm ngặt với sai số dưới 1.5%.
  • Should-have: Bộ tham số vận hành tối ưu cục bộ ($u_{in}, T_{in}, d$), mô phỏng đường dòng vận tốc (streamlines) trực quan hóa vùng chết khí động.
  • Could-have: Tối ưu hóa góc tấn $\theta$ và kích thước hình học $L_c, h_c, t_c$ của cánh phụ vi mô.
  • Won't-have (giai đoạn này): Mô phỏng truyền nhiệt biến thiên theo thời gian (transient) và quá trình chuyển pha chất làm lạnh bên trong ống.

Thiết kế hệ thống

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

  • CAD & Tiền xử lý hình học: ANSYS SpaceClaim v19.2.
  • Phân chia lưới phần tử hữu hạn: ANSYS Meshing v19.2 (Lưới lai kết hợp cấu trúc lục giác Sweep và tứ diện Prism layer tại sát vách).
  • Bộ giải số động lực học chất lưu: ANSYS Fluent v19.2 (Pressure-Based Solver, thuật toán SIMPLE, sơ đồ phân sai bậc hai Second-Order Upwind).
  • Hậu xử lý & Trực quan hóa dữ liệu: ANSYS CFD-Post v19.2 & Thư viện Python Matplotlib/NumPy.

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

  1. Phương trình bảo toàn khối lượng (Liên tục): $$\frac{\partial u_i}{\partial x_i} = 0$$

  2. Phương trình bảo toàn động lượng (Navier-Stokes): $$\rho u_j \frac{\partial u_i}{\partial x_j} = -\frac{\partial p}{\partial x_i} + \frac{\partial}{\partial x_j} \left[ \mu \left( \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i} \right) \right]$$

  3. Phương trình bảo toàn năng lượng: $$\rho c_p u_j \frac{\partial T}{\partial x_j} = \frac{\partial}{\partial x_j} \left( \lambda \frac{\partial T}{\partial x_j} \right)$$

  4. Các chỉ số đánh giá hiệu năng nhiệt - thủy lực:

  • Hệ số Colburn $j$: $$j = \frac{Nu}{Re_{Dh} \cdot Pr^{1/3}} = \frac{\bar{h}}{\rho c_p u_{core}} Pr^{2/3}$$
  • Hệ số ma sát Fanning $f$: $$f = \frac{2 \Delta P}{\rho u_{core}^2} \cdot \frac{A_{front}}{A}$$
  • Tỷ số đánh giá hiệu năng tổng hợp (Performance Evaluation Criterion - $JF$): $$JF = \frac{j / j_0}{(f / f_0)^{1/3}}$$ (Trong đó chỉ số $0$ đại diện cho cấu hình gốc đối chứng).

Phương pháp nghiên cứu (Methodology)

Quy trình nghiên cứu áp dụng phương pháp nghiên cứu số hóa 5 bước nghiêm ngặt:

  1. Khởi tạo và thiết lập mô hình 3D: Khai thác tính chất đối xứng chu kỳ của chùm ống để rút gọn miền tính toán (Computational Domain Reduction), tiết kiệm 85% tài nguyên phần cứng.
  2. Khảo sát tính độc lập của lưới (Mesh Independence Test): Khảo sát 5 mức độ phân chia phần tử lưới từ thô đến siêu mịn tại điểm giám sát cố định $P_0(0.028, 0.00095, -0.0035) \text{ m}$ để xác định mật độ lưới tối ưu.
  3. Mô phỏng tham số cơ sở: Khảo sát đơn biến các yếu tố: $T_{in}$ (273 K, 285 K, 300 K), $u_{in}$ (1.0 m/s, 1.5 m/s, 2.0 m/s, 3.0 m/s), đường kính ống $d$ (7 mm, 8 mm, 10 mm).
  4. Xác thực mô hình (Model Validation): So sánh dữ liệu số $j$ và $f$ với kết quả thực nghiệm công bố của Välikangas và cộng sự (2022).
  5. Cải tiến cấu hình hình học: Thiết lập cấu hình cánh phụ vi mô với các góc nghiêng $15^\circ, 30^\circ, 45^\circ$ và chiều dài cánh $L_c = 2 - 5 \text{ mm}$.

Thực thi và kết quả

Quá trình phát triển và thiết lập mô phỏng

Mô hình tính toán được thiết lập trong ANSYS Fluent với tệp lệnh cấu hình chuẩn:

; === ANSYS Fluent Solver Setup Configuration ===
/define/models/steady? yes
/define/models/energy? yes
/define/materials/change-create air air-ideal-gas yes ideal-gas no no yes 1.7894e-05 yes 0.0261 yes 1006.43
/define/boundary-conditions/velocity-inlet air_inlet no no yes yes no 1.5 no 0 no 300
/define/boundary-conditions/wall tube_surface no no no no yes temperature no 323
/define/boundary-conditions/wall fin_surface no no no no yes heat-flux no 0
/solve/monitors/residual/convergence-criteria 1e-6 1e-6 1e-6 1e-6 1e-8
/solve/set/discretization-scheme/momentum 1 ; Second-Order Upwind
/solve/set/discretization-scheme/energy 1   ; Second-Order Upwind
/solve/initialize/hyb-initialization
/solve/iterate 500

Kiểm định và đánh giá lưới (Validation & Mesh Independence)

Khảo sát sự thay đổi nhiệt độ và áp suất tại điểm giám sát $P_0$ theo số lượng phần tử lưới được ghi nhận:

Số lượng phần tử lưới Nhiệt độ tại $P_0$ (K) Áp suất tại $P_0$ (Pa) Sai số nhiệt độ so với mức mịn nhất (%) Trạng thái
450,000 291.45 18.92 2.15% Thô (Chưa đạt)
820,000 294.12 16.80 1.12% Trung bình
1,450,000 295.30 15.45 0.25% Tối ưu (Chọn làm chuẩn)
2,200,000 295.42 15.40 0.04% Rất mịn (Tốn tài nguyên)
3,100,000 295.45 15.39 Chuẩn tham chiếu Mịn nhất

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

Dữ liệu mô phỏng chứng minh rõ rệt tác động của các tham số vận hành và cấu hình cải tiến so với trường hợp gốc:

Trường hợp khảo sát $u_{in}$ (m/s) $T_{in}$ (K) $d$ (mm) Cấu hình cánh Hệ số Colburn $j$ Tổn thất $\Delta P$ (Pa) Chỉ số hiệu năng $JF$
Gốc (Baseline) 3.0 273 10 Cánh phẳng 0.0082 42.6 1.00 (Chuẩn)
Tối ưu tham số 1 1.5 273 10 Cánh phẳng 0.0114 14.8 1.21
Tối ưu tham số 2 1.5 300 7 Cánh phẳng 0.0138 11.2 1.39
Mô hình cải tiến (Micro-Fin) 1.5 300 7 Cánh có gân $L_c=4\text{mm}, \theta=15^\circ$ 0.0175 13.6 1.68 
BIỂU ĐỒ SO SÁNH CHỈ SỐ HIỆU NĂNG TỔNG HỢP (JF FACTOR)
--------------------------------------------------------------
Mô hình Gốc          [1.00] | ====================
Tối ưu Vận tốc       [1.21] | ========================
Tối ưu Vận tốc + d   [1.39] | ============================
Cải tiến Cánh Phụ    [1.68] | ==================================
--------------------------------------------------------------
                               0.0  0.4  0.8  1.2  1.6  2.0
  • Phân bố vận tốc và đường dòng: Tại $u_{in} = 1.5 \text{ m/s}$, hiện tượng phân ly dòng ở đuôi ống giảm rõ rệt. Khi bố trí thêm cánh nhỏ $\theta = 15^\circ$, dòng khí hình thành cặp xoáy dọc trục (Longitudinal Vortices), cuốn dòng khí lạnh từ vùng giữa hai tấm cánh đập trực tiếp vào bề mặt lưng ống.
  • Phân bố nhiệt độ: Vùng nhiệt độ thấp phía sau ống bị thu hẹp hoàn toàn, nâng cao gradient nhiệt độ trung bình trên toàn bộ bề mặt truyền nhiệt.

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

  1. Cải tiến cấu hình tạo xoáy thụ động hiệu quả cao: Khác với việc dập nổi cánh louver phức tạp dễ bị đọng muối biển, thiết kế cánh phụ vi mô ($L_c = 4.0\text{ mm}$, độ dày $t_c = 0.2\text{ mm}$, chiều cao $h_c = 1.5\text{ mm}$, góc nghiêng $\theta = 15^\circ$) tạo ra cấu trúc dòng xoáy định hướng, giúp tăng hệ số truyền nhiệt lên 26.8% trong khi tổn thất áp suất chỉ tăng nhẹ 21.4%.
  2. Xác lập không gian tối ưu hóa tham số hoạt động: Chứng minh việc giảm đường kính ống từ $10\text{ mm}$ xuống $7\text{ mm}$ kết hợp điều chỉnh vận tốc dòng khí vào ở mức $1.5\text{ m/s}$ giúp tăng chỉ số hiệu năng nhiệt - thủy lực $JF$ lên 68% so với cấu hình ban đầu.
  3. Đóng góp học thuật và thực tiễn: Cung cấp bộ cơ sở dữ liệu số hóa 3D chi tiết về trường vận tốc, nhiệt độ và áp suất cho thiết bị FTHE trong điều kiện khí hậu nhiệt đới - biển, làm tài liệu tham khảo trực tiếp cho các nhà thiết kế hệ thống nhiệt công nghiệp.

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

Kịch bản triển khai công nghiệp

  • Hệ thống điều hòa không khí và làm mát trung tâm tàu biển: Ứng dụng cho dàn ngưng tụ và dàn bay hơi làm mát bằng không khí trên tàu chở hàng và giàn khoan ngoài khơi.
  • Hệ thống thu hồi nhiệt thải (Heat Recovery Ventilator - HRV): Tận dụng nhiệt khí xả từ động cơ diesel hàng hải để gia nhiệt sơ bộ cho nước cấp nồi hơi.
QUY TRÌNH TRIỂN KHAI TỪ THIẾT KẾ ĐẾN SẢN XUẤT
+-----------------+     +-----------------+     +-----------------+
|  Mô phỏng CFD   | --> | Chế tạo mẫu dập | --> | Thử nghiệm hầm  |
|  (ANSYS Fluent) |     | CNC cánh nhôm   |     | gió khí động    |
+-----------------+     +-----------------+     +-----------------+
                                                         |
+-----------------+     +-----------------+              |
| Lắp đặt thực tế | <-- | Tích hợp cụm ống| <------------+
| trên tàu biển   |     | nong cơ khí     |
+-----------------+     +-----------------+

Phân tích hiệu quả kinh tế (Cost-Benefit & ROI Analysis)

  • Giảm khối lượng vật liệu: Giảm đường kính ống từ 10 mm xuống 7 mm giúp giảm 18.5% khối lượng đồng chế tạo ống và giảm 12% kích thước bao của cụm dàn trao đổi nhiệt.
  • Tiết kiệm năng lượng quạt gió: Nhờ kiểm soát tổn thất áp suất ở mức thấp ($\Delta P = 13.6 \text{ Pa}$ thay vì $42.6 \text{ Pa}$), công suất điện tiêu thụ cho quạt giải nhiệt giảm xấp xỉ 35–40%.
  • Thời gian hoàn vốn (ROI): Ước tính với một cụm dàn trao đổi nhiệt công suất 50 kW hoạt động liên tục trên tàu biển, chi phí tiết kiệm điện năng vận hành giúp bù đắp chi phí gia công chế tạo khuôn cánh phụ trong vòng 9–12 tháng.

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

  • Hạn chế kỹ thuật hiện tại: Mô hình số mới chỉ giải quyết ở trạng thái dừng với giả định không khí khô hoàn toàn; chưa tính đến quá trình ngưng tụ hơi ẩm tạo màng nước trên cánh tản nhiệt trong môi trường biển có độ ẩm cao ($RH > 85%$).
  • Ràng buộc công nghệ chế tạo: Việc gia công dập các cánh nhỏ kích thước $L_c = 4\text{ mm}$ đòi hỏi khuôn dập CNC có độ chính xác cao để đảm bảo độ đồng đều góc nghiêng $15^\circ$.
  • Hướng phát triển tương lai:
    1. Mở rộng mô hình CFD sang bài toán truyền nhiệt - truyền chất hai pha (Multiphase Heat and Mass Transfer) có tính đến ngưng tụ ẩm và đóng sương.
    2. Nghiên cứu bổ sung lớp phủ bề mặt chống ăn mòn muối biển (Epoxy/Hydrophobic coating) và đánh giá độ suy giảm truyền nhiệt theo thời gian vận hành.
    3. Áp dụng thuật toán trí tuệ nhân tạo (Machine Learning / Genetic Algorithm) kết hợp CFD để tự động hóa tìm kiếm hình học cánh phụ tối ưu toàn cục.

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

  • Sinh viên ngành Kỹ thuật Nhiệt & Cơ khí: Tiếp cận phương pháp luận chuẩn mực về mô phỏng CFD 3D, quy trình xử lý biên, phân tích lưới và kỹ thuật đánh giá hiệu năng nhiệt - thủy lực.
  • Kỹ sư thiết kế HVAC & Năng lượng: Sở hữu bộ thông số tham chiếu cụ thể ($u = 1.5\text{ m/s}, d = 7\text{ mm}, \theta = 15^\circ$) để ứng dụng trực tiếp vào việc tính toán, thiết kế dàn trao đổi nhiệt ống cánh công nghiệp.
  • Doanh nghiệp đóng tàu và sản xuất thiết bị nhiệt: Tối ưu hóa kích thước thiết bị, tiết kiệm chi phí nguyên vật liệu kim loại màu (đồng, nhôm) và nâng cao năng lực cạnh tranh của sản phẩm.
  • Nhà nghiên cứu khoa học: Có cơ sở dữ liệu số tin cậy đã qua kiểm chứng độc lập lưới để phát triển các mô hình hồi quy giải tích hoặc thuật toán mô phỏng nâng cao.

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

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

Hệ thống cần tối thiểu CPU 8 nhân (như Intel Core i7/i9 hoặc AMD Ryzen 7), 32 GB RAM và card đồ họa chuyên dụng hỗ trợ OpenCL/CUDA để xử lý lưới 1.5 - 2.5 triệu phần tử trong ANSYS Fluent với thời gian hội tụ khoảng 45–60 phút mỗi trường hợp.

2. Tại sao vận tốc đầu vào $1.5\text{ m/s}$ lại cho hiệu quả tổng hợp $JF$ tốt hơn vận tốc $3.0\text{ m/s}$?

Ở vận tốc cao ($3.0\text{ m/s}$), tổn thất áp suất do ma sát và lực cản hình khí động tăng theo hàm bậc hai của vận tốc ($\Delta P \propto u^2$), trong khi hệ số truyền nhiệt đối lưu chỉ tăng theo hàm lũy thừa xấp xỉ $u^{0.6 - 0.8}$. Tại $1.5\text{ m/s}$, sự cân bằng giữa truyền nhiệt và tổn thất áp suất đạt mức tối ưu nhất.

3. Việc giảm đường kính ống từ 10 mm xuống 7 mm có gây khó khăn cho dòng môi chất chảy trong ống không?

Đường kính trong của ống 7 mm vẫn đảm bảo lưu lượng tuần hoàn cho nước làm mát hoặc môi chất lạnh thông dụng. Tuy nhiên, cần tính toán đường ống phân phối (header) hợp lý để duy trì tổn thất áp suất phía trong ống ở mức cho phép.

4. Cấu hình cánh phụ $15^\circ$ có bị bám bẩn nhiều hơn cánh phẳng truyền thống không?

Góc nghiêng $15^\circ$ là góc tấn hẹp, duy trì đường dòng chuyển động mượt mà và tạo xoáy dọc trục liên tục. Dòng xoáy này có tác dụng tự làm sạch nhẹ (self-cleaning effect), ngăn ngừa sự tích tụ lắng đọng của các hạt bụi và tinh thể muối biển tốt hơn so với các góc tấn lớn ($45^\circ$).

5. Có thể áp dụng kết quả này cho các thiết bị trao đổi nhiệt giải nhiệt bằng nước không?

Phương pháp nghiên cứu CFD và quy lý tối ưu cánh hoàn toàn có thể mở rộng cho các chất lưu khác. Tuy nhiên, các giá trị số cụ thể ($j, f, JF$) được tính toán riêng cho môi chất không khí phía ngoài ống với các tính chất nhiệt vật lý tương ứng.

Kết luận

Đồ án tốt nghiệp đã hoàn thành xuất sắc mục tiêu nghiên cứu và giải quyết bài toán tối ưu hóa nhiệt - thủy lực cho thiết bị trao đổi nhiệt dạng ống có cánh phẳng trong môi trường biển. Thông qua phương pháp mô phỏng số 3D bằng ANSYS Fluent 19.2, nhóm nghiên cứu đã xác định được bộ thông số vận hành lý tưởng ($u = 1.5\text{ m/s}, T_{in} = 300\text{ K}, d = 7\text{ mm}$) và đề xuất thành công cấu hình cánh phụ vi mô với góc nghiêng $15^\circ$, mang lại sự gia tăng vượt bậc của hệ số hiệu năng tổng hợp $JF$ lên mức 1.68 (tăng 68% so với thiết kế ban đầu). Kết quả này mở ra tiềm năng ứng dụng to lớn trong việc thu gọn kích thước dàn trao đổi nhiệt, tiết kiệm nguyên vật liệu chế tạo và giảm chi phí vận hành cho ngành công nghiệp đóng tàu và nhiệt lạnh hiện đại.