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

Thiết bị trao đổi nhiệt là thành phần kỹ thuật quan trọng trong các hệ thống điều hòa không khí (HVAC), hệ thống lạnh công nghiệp, các nhà máy nhiệt điện, hóa chất và đặc biệt là hệ thống động lực, làm mát trên tàu thủy thuộc ngành công nghiệp biển. Tại các công trình biển và tàu thuyền, thiết bị trao đổi nhiệt phải hoạt động trong điều kiện môi trường đặc thù với các yếu tố như độ mặn, nhiệt độ, áp suất, chất lượng nước và tác động động lực học của dòng chảy biển. Việc tối ưu hóa khả năng truyền nhiệt kết hợp giảm thiểu tổn thất áp suất có ý nghĩa quyết định đến việc nâng cao hiệu suất làm việc, tiết kiệm năng lượng và giảm chi phí vận hành hệ thống.

Đồ án tốt nghiệp "Nghiên cứu sự ảnh hưởng của thông số hoạt động và cấu hình đến hiệu quả truyền nhiệt của thiết bị trao đổi nhiệt dạng ống có cánh bằng phương pháp mô phỏng số" được thực hiện bởi nhóm sinh viên Nguyễn Trịnh Hoài My, Trương Sỹ Tùng Phương, Huỳnh Thanh Phong, Nguyễn Hùng Sơn dưới sự hướng dẫn của TS. Trần Thanh Tình tại Bộ môn Công nghệ Kỹ thuật Nhiệt – Điện lạnh, Khoa Cơ khí Động lực, Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh.

Mục tiêu tổng quát của đề tài là khảo sát ảnh hưởng của các thông số hoạt động và thông số hình học đến hiệu quả nhiệt – thủy lực, từ đó tìm kiếm bộ thông số tối ưu và đề xuất cải tiến cấu hình cánh cho thiết bị trao đổi nhiệt dạng ống có cánh tản nhiệt phẳng (FTHE) hoạt động trong môi trường biển.

Nhiệm vụ nghiên cứu được văn bản xác định cụ thể qua hai nội dung:

  1. Nghiên cứu, tìm hiểu về thiết bị trao đổi nhiệt dạng ống có cánh, từ đó tiến hành khảo sát các đặc tính nhiệt - thủy lực dựa trên phương pháp mô phỏng số (CFD).
  2. Rút ra nhận xét từ kết quả khảo sát mô phỏng số, từ đó phát triển các giải pháp tối ưu hóa thông số hoạt động và cải tiến cấu hình hình học cho thiết bị trao đổi nhiệt dạng ống có cánh.

Đối tượng nghiên cứu là thiết bị trao đổi nhiệt dạng ống có cánh (Finned-tube heat exchanger - FTHE), tập trung vào miền không khí giải nhiệt chuyển động cắt ngang qua chùm ống có gắn cánh tản nhiệt. Phạm vi nghiên cứu về mặt không gian là một phân đoạn đại diện có tính chất đối xứng và tuần hoàn của bộ trao đổi nhiệt ba chiều (3D), mô phỏng quá trình truyền nhiệt và dòng chảy chất lưu phía không khí trong điều kiện môi trường biển; phạm vi thời gian thực hiện đồ án từ ngày 06 tháng 03 năm 2023 đến ngày 17 tháng 07 năm 2023.

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

Cơ sở lý luận của đề tài dựa trên các nguyên lý truyền nhiệt (dẫn nhiệt qua vách ống và cánh tản nhiệt, tỏa nhiệt đối lưu giữa bề mặt rắn và dòng không khí) cùng hệ phương trình động lực học chất lưu (phương trình Navier-Stokes cho dòng liên tục, bảo toàn động lượng và bảo toàn năng lượng). Hiệu quả làm việc của thiết bị được đánh giá thông qua các chỉ số nhiệt - thủy lực không thứ nguyên bao gồm: số Reynolds ($Re$, $Re_{Dh}$ tính theo đường kính thủy lực $D_h$), số Nusselt ($Nu$), số Prandtl ($Pr$), hệ số truyền nhiệt trung bình tổng thể $\bar{h}$, hệ số truyền nhiệt Colburn $j$, hệ số ma sát Fanning $f$, độ giảm áp $\Delta P$, và đặc biệt là tỷ số hiệu suất $JF$ (tiêu chí đánh giá hiệu năng tổng hợp giữa tăng cường truyền nhiệt và tổn thất áp suất do ma sát).

Phương pháp nghiên cứu chính là mô phỏng số động lực học chất lưu (Computational Fluid Dynamics - CFD) kết hợp phương pháp thể tích hữu hạn (Finite Volume Method - FVM) được thực thi trên bộ phần mềm ANSYS phiên bản 19.2, bao gồm các mô-đun:

  • ANSYS SpaceClaim: Xây dựng và xử lý hình học miền tính toán 3D của thiết bị trao đổi nhiệt.
  • ANSYS Meshing: Thiết lập cấu trúc lưới tính toán, tạo lớp biên tại các bề mặt trao đổi nhiệt.
  • ANSYS Fluent: Thiết lập phương trình điều khiển, lựa chọn mô hình dòng chảy, gán điều kiện biên và tiến hành giải số lặp.
  • Công cụ Parameter: Thực hiện kiểm tra tính độc lập của lưới (hội tụ lưới) và tự động hóa khảo sát các trường hợp thay đổi thông số.
  • ANSYS CFD-Post: Trích xuất dữ liệu số học, vẽ trường phân bố nhiệt độ, trường áp suất, trường vận tốc và các đường dòng (streamlines).

Nguồn dữ liệu tham chiếu và xác thực mô hình được kế thừa từ công trình nghiên cứu của Turo Välikangas, Mikko Folkersma, Miikka Dal Mas, Tuomo Keskitalo, Petteri Peltonen và Ville Vuorinen công bố năm 2022 trên tạp chí khoa học quốc tế thuộc nhà xuất bản Elsevier với tiêu đề: "Parametric CFD study for finding the optimal tube arrangement of a fin-and-tube heat exchanger with plain fins in a marine environment".

Yếu tố phương pháp Chi tiết thiết lập trong đồ án
Nền tảng phần mềm ANSYS 19.2 (SpaceClaim, Meshing, Fluent, Parameter, CFD-Post)
Phương pháp số Phương pháp thể tích hữu hạn (FVM), mô phỏng 3D dòng ổn định
Chỉ số đánh giá Hệ số Colburn $j$, hệ số ma sát Fanning $f$, tỷ số hiệu suất $JF$, tổn thất áp suất $\Delta P$, số $Re$
Đối tượng hình học Miền không khí giữa hai tấm cánh cắt qua cụm ống tròn có cánh phẳng
Nguồn đối chuẩn Công trình mô phỏng CFD FTHE của Välikangas và cộng sự (Elsevier, 2022)

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

CHƯƠNG I. Lý do chọn đề tài

Chương này trình bày tổng quan về các dạng thiết bị trao đổi nhiệt công nghiệp, phân tích ưu nhược điểm của ba nhóm kết cấu chính: thiết bị trao đổi nhiệt kiểu hồi nhiệt, thiết bị trao đổi nhiệt kiểu hỗn hợp và thiết bị trao đổi nhiệt dạng vách ngăn cách. Nhóm tác giả nhấn mạnh thiết bị kiểu vách ngăn cách chiếm hơn 98% ứng dụng trong thực tế sản xuất nhờ khả năng giữ độ tinh khiết cho các lưu chất và cho phép làm việc ở độ chênh lệch áp suất lớn. Tiếp đó, cấu tạo chi tiết của thiết bị trao đổi nhiệt dạng ống có cánh (FTHE) được phân tích gồm hai phần chính: ống dẫn môi chất mang nhiệt và các tấm cánh mỏng nhằm gia tăng diện tích tiếp xúc nhiệt.

Chương I cũng tổng thuật 10 công trình khoa học tiêu biểu trong và ngoài nước liên quan đến FTHE (Välikangas và cộng sự, Lindqvist và cộng sự, Wang và cộng sự, Batista và cộng sự, Wen và cộng sự, Ameel và cộng sự, Bhuiyan và cộng sự, Wu và cộng sự, Xie và cộng sự, Borrajo-Peláez và cộng sự), phân tích các yếu tố như góc tấn của bộ tạo xoáy cánh nhỏ delta, bước ống ngang, bước ống dọc, bước cánh và đường kính ống. Trên cơ sở đó, đồ án lựa chọn nghiên cứu của Välikangas và cộng sự (2022) làm công trình chuẩn để kế thừa mô hình và đối sánh kết quả.

Kiểu thiết bị trao đổi nhiệt Nguyên lý làm việc Ưu điểm chính Nhược điểm chính
Kiểu hồi nhiệt (HRV) Bề mặt rắn tích và nhả nhiệt theo chu kỳ khi lần lượt tiếp xúc dòng nóng và lạnh Ít tốn kém vật tư, kết cấu đơn giản, chi phí thấp hơn kiểu bề mặt Làm việc không ổn định theo chu kỳ, mức độ vệ sinh không cao, áp suất hai dòng phải tương đương
Kiểu hỗn hợp Hai chất lỏng tiếp xúc trực tiếp và hòa trộn với nhau để truyền nhiệt Kết cấu đơn giản, giá thành rẻ, hiệu suất cao, quán tính nhiệt rất nhỏ Hai chất lỏng bị hòa trộn lẫn nhau, hạn chế môi trường làm việc
Kiểu vách ngăn cách (gồm FTHE) Truyền nhiệt qua bề mặt vách ngăn cách cố định bằng dẫn nhiệt và đối lưu Đảm bảo độ tinh khiết của lưu chất, chịu được chênh lệch áp suất lớn Tiêu tốn nhiều vật tư, chế tạo phức tạp, quán tính nhiệt lớn

CHƯƠNG II. XÂY DỰNG MÔ HÌNH VÀ MÔ PHỎNG

Chương II mô tả việc thiết lập mô hình tính toán số học trên không gian 3D. Do tính chất đối xứng và lặp lại tuần hoàn của chùm ống và các tấm cánh phẳng, nhóm tác giả xây dựng một miền tính toán đối xứng đại diện cho dòng không khí chảy qua các ống. Miền tính toán có phương song song với dòng không khí đầu vào và vuông góc với trục các đường ống.

Các điều kiện biên được cài đặt trên phần mềm ANSYS Fluent bao gồm:

  • Vận tốc đầu vào (Velocity inlet): Cố định giá trị vận tốc dòng không khí $u$ và nhiệt độ đầu vào $T_{in}$.
  • Đầu ra dòng chảy (Pressure outlet): Áp suất tĩnh bằng áp suất khí quyển.
  • Bề mặt ống: Gán điều kiện nhiệt độ bề mặt vách không đổi ($T_w$).
  • Bề mặt cánh tản nhiệt: Gán điều kiện thông lượng nhiệt bề mặt phù hợp với trao đổi nhiệt đối lưu và dẫn nhiệt.
  • Các mặt biên bên và mặt trên/dưới: Thiết lập điều kiện biên đối xứng (Symmetry) để đại diện cho toàn bộ chùm ống trong thiết bị thực tế.

Quy trình khảo sát tính độc lập của lưới được thực hiện thông qua công cụ Parameter với nhiều mức độ phân giải lưới khác nhau. Tiêu chí đánh giá hội tụ lưới dựa trên sự biến thiên nhiệt độ đầu ra, áp suất đầu ra và các giá trị nhiệt độ, áp suất cục bộ tại một điểm kiểm tra có tọa độ cụ thể $(x = 0{,}028\text{ m}; y = 0{,}00095\text{ m}; z = -0{,}0035\text{ m})$ cùng nhiệt độ trung bình toàn mô hình.

CHƯƠNG III. KẾT QUẢ VÀ THẢO LUẬN

Chương III trình bày toàn bộ kết quả mô phỏng số, đồ thị khảo sát tham số và phân tích các trường vật lý:

  1. Mô hình cơ sở ban đầu: Thiết lập với vận tốc đầu vào $u = 3\text{ m/s}$, nhiệt độ đầu vào $T_{in} = 273\text{ K}$ và đường kính ống $d = 10\text{ mm}$. Kết quả cho thấy tại vận tốc này, trở lực dòng chảy và tổn thất áp suất qua chùm ống còn cao, hiệu quả truyền nhiệt tổng thể chưa đạt trạng thái tối ưu nhiệt - thủy lực.
  2. Khảo sát ảnh hưởng của các thông số vận hành và cấu hình:
    • Ảnh hưởng của nhiệt độ đầu vào $T_{in}$: Khảo sát sự biến thiên của nhiệt độ dòng khí vào đến tỷ số hiệu suất $JF$, xác định trường phân bố nhiệt và giá trị $T_{core}$ tại mặt cắt nhỏ nhất.
    • Ảnh hưởng của vận tốc đầu vào $u$: Khảo sát dải vận tốc dòng chảy, đánh giá sự thay đổi của độ giảm áp $\Delta P$, hệ số Fanning $f$, hệ số Colburn $j$ và vận tốc $u_{core}$. Kết quả chỉ ra khi giảm vận tốc xuống mức thích hợp, tổn thất áp suất giảm đáng kể giúp tỷ số $JF$ tăng cao.
    • Ảnh hưởng của đường kính ống $d$: Khảo sát sự thay đổi kích thước ống đến sự cản trở dòng chảy và hiệu quả trao đổi nhiệt.
  3. Kết quả tối ưu hóa cục bộ: Xác định bộ thông số làm việc tối ưu của mô hình là vận tốc đầu vào $u = 1{,}5\text{ m/s}$, nhiệt độ đầu vào $T_{in} = 300\text{ K}$ và đường kính ống $d = 7\text{ mm}$.
  4. So sánh với bài báo tham khảo: Kết quả mô phỏng trường nhiệt độ, trường áp suất và đường dòng của mô hình tối ưu được đối chiếu với mô hình hồi quy và dữ liệu từ công trình của Välikangas và cộng sự (2022), cho thấy sự tương đồng về xu hướng khí động và phân bố nhiệt lượng.
  5. Cải tiến cấu hình hình học tấm cánh: Để nâng cao hơn nữa hiệu suất, nhóm tác giả đề xuất giải pháp bổ sung các cánh nhỏ (vortex generator/cánh phụ) bố trí giữa miền không khí và hai cánh tản nhiệt chính. Kết quả khảo sát kích thước hình học cho thấy cấu hình đạt hiệu quả cao nhất có chiều dài cánh nhỏ $L_c = 4\text{ mm}$ (so sánh với trường hợp $L_c = 3\text{ mm}$), độ dày cánh $0{,}2\text{ mm}$, chiều cao cánh $1{,}5\text{ mm}$ và góc nghiêng hợp bởi cánh là $15^\circ$. Cấu hình này tạo ra các xoáy lốc cục bộ làm mỏng lớp biên nhiệt nhưng không gây tăng quá mức tổn thất áp suất.

CHƯƠNG IV. KẾT LUẬN VÀ KIẾN NGHỊ

Chương IV tổng kết các kết quả đạt được về việc tối ưu hóa thông số vận hành và cấu hình hình học của thiết bị trao đổi nhiệt FTHE, đồng thời nêu các kiến nghị nhằm mở rộng phạm vi nghiên cứu thực nghiệm trong các giai đoạn tiếp theo.

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

Nghiên cứu mang lại các kết quả định lượng cụ thể thông qua mô phỏng số:

  • Xác định bộ thông số vận hành và hình học ban đầu ($u = 3\text{ m/s}$, $T_{in} = 273\text{ K}$, $d = 10\text{ mm}$) chưa đạt hiệu suất tối ưu do tổn thất áp suất lớn.
  • Tìm ra bộ thông số tối ưu hóa cục bộ cho thiết bị: vận tốc đầu vào $u = 1{,}5\text{ m/s}$, nhiệt độ đầu vào $T_{in} = 300\text{ K}$ và đường kính ống giảm xuống $d = 7\text{ mm}$, giúp cải thiện rõ rệt tỷ số hiệu suất nhiệt - thủy lực $JF$.
  • Đề xuất mô hình cải tiến cấu hình cánh tản nhiệt bằng việc tích hợp các cánh nhỏ phụ giữa hai tấm cánh chính với thông số hình học xác định: chiều dài cánh $L_c = 4\text{ mm}$, độ dày $0{,}2\text{ mm}$, chiều cao $1{,}5\text{ mm}$ và góc đặt $15^\circ$.
Thông số Mô hình ban đầu Mô hình tối ưu cục bộ Mô hình cải tiến cánh
Vận tốc đầu vào ($u$) $3\text{ m/s}$ $1{,}5\text{ m/s}$ $1{,}5\text{ m/s}$
Nhiệt độ đầu vào ($T_{in}$) $273\text{ K}$ $300\text{ K}$ $300\text{ K}$
Đường kính ống ($d$) $10\text{ mm}$ $7\text{ mm}$ $7\text{ mm}$
Kết cấu cánh tản nhiệt Cánh phẳng tiêu chuẩn Cánh phẳng tiêu chuẩn Bổ sung cánh nhỏ ($L_c = 4\text{ mm}$, dày $0{,}2\text{ mm}$, cao $1{,}5\text{ mm}$, góc $15^\circ$)
Đánh giá hiệu năng Tổn thất áp suất cao, $JF$ chưa tối ưu Giảm độ giảm áp $\Delta P$, tăng tỷ số $JF$ Tăng cường hòa trộn dòng, tối ưu hóa trao đổi nhiệt

Đóng góp của đề tài thể hiện ở việc ứng dụng thành công chu trình mô phỏng CFD trên phần mềm ANSYS để giải quyết bài toán nhiệt - thủy lực của thiết bị trao đổi nhiệt ống có cánh định hướng môi trường biển. Đồ án đưa ra dữ liệu đối chuẩn tin cậy với công trình quốc tế của Välikangas và cộng sự (2022), đồng thời chứng minh giải pháp bổ sung cánh nhỏ góc $15^\circ$ là phương án hình học khả thi để nâng cao hiệu quả trao đổi nhiệt.

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

Hạn chế của đề tài xuất phát từ phạm vi mô phỏng:

  • Toàn bộ kết quả dựa trên phương pháp tính toán số học trên miền tính toán đối xứng lý tưởng, chưa bao quát toàn bộ cụm thiết bị trao đổi nhiệt quy mô công nghiệp đầy đủ.
  • Nghiên cứu chưa thực hiện chế tạo mẫu thực nghiệm vật lý để đo đạc kiểm chứng trực tiếp các thông số nhiệt và tổn thất áp suất trong điều kiện phòng thí nghiệm.
  • Mô hình chưa xét đến hiện tượng bám bẩn (fouling), đóng cặn muối hoặc ăn mòn hóa học bề mặt kim loại trong quá trình vận hành lâu dài ở môi trường biển thực tế.

Hướng nghiên cứu tiếp theo được tác giả định hướng:

  • Chế tạo mô hình thực nghiệm thực tế để kiểm chứng đối sánh dữ liệu đo đạc với kết quả mô phỏng CFD.
  • Mở rộng khảo sát với các dạng hình học cánh phức hợp khác và nghiên cứu độ bền vật liệu trong điều kiện môi trường biển có độ mặn cao.

Giá trị tham khảo

Đồ án là tài liệu tham khảo hữu ích cho 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 Lạnh - Điều hòa không khí, Kỹ thuật Cơ khí Động lực và Kỹ thuật Tàu thủy. Các nội dung có giá trị tham khảo chuyên môn cao bao gồm:

  • Quy trình chuẩn xây dựng miền tính toán 3D đối xứng và phương pháp phân tích hội tụ lưới (Mesh Independence Test) trên ANSYS Meshing và Parameter.
  • Phương pháp thiết lập điều kiện biên nhiệt - thủy lực và xử lý số liệu trường dòng chảy bằng ANSYS Fluent và CFD-Post.
  • Bộ số liệu tham chiếu về kích thước cánh nhỏ tạo xoáy ($L_c = 4\text{ mm}$, cao $1{,}5\text{ mm}$, góc nghiêng $15^\circ$) phục vụ thiết kế các bộ trao đổi nhiệt ống cánh nhỏ gọn.

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

1. Mục tiêu cốt lõi của đề tài là gì?
Nghiên cứu tập trung khảo sát sự ảnh hưởng của các thông số vận hành (vận tốc, nhiệt độ đầu vào) và thông số hình học (đường kính ống, dạng cánh) đến hiệu quả truyền nhiệt và tổn thất áp suất của thiết bị trao đổi nhiệt dạng ống có cánh phẳng trong môi trường biển, từ đó tìm kiếm thông số tối ưu và cải tiến cấu hình cánh.

2. Đồ án sử dụng công trình khoa học nào làm cơ sở đối chuẩn và kiểm chứng mô hình CFD?
Đề tài kế thừa và đối chuẩn dữ liệu với bài báo khoa học của Turo Välikangas và cộng sự xuất bản năm 2022 trên tạp chí thuộc nhà xuất bản Elsevier, nghiên cứu về tối ưu hóa bố trí ống cho bộ trao đổi nhiệt FTHE có cánh phẳng trong môi trường biển.

3. Các thông số làm việc ban đầu và bộ thông số tối ưu cục bộ tìm được trong nghiên cứu là bao nhiêu?
Mô hình ban đầu được thiết lập tại $u = 3\text{ m/s}$, $T_{in} = 273\text{ K}$ và $d = 10\text{ mm}$. Sau quá trình khảo sát tham số, bộ thông số tối ưu cục bộ được xác định là vận tốc đầu vào $u = 1{,}5\text{ m/s}$, nhiệt độ đầu vào $T_{in} = 300\text{ K}$ và đường kính ống $d = 7\text{ mm}$.

4. Cấu hình cánh tản nhiệt cải tiến được nhóm tác giả đề xuất có thông số hình học cụ thể như thế nào?
Nhóm tác giả đề xuất bổ sung các cánh nhỏ bố trí giữa miền không khí và hai cánh chính với chiều dài cánh $L_c = 4\text{ mm}$, độ dày $0{,}2\text{ mm}$, chiều cao $1{,}5\text{ mm}$ và góc nghiêng hợp bởi cánh là $15^\circ$.

5. Việc kiểm tra tính độc lập của lưới trong đồ án được căn cứ trên những tiêu chí nào?
Kiểm tra hội tụ lưới được thực hiện qua công cụ Parameter, theo dõi độ biến thiên của nhiệt độ và áp suất tại đầu ra hệ thống, nhiệt độ trung bình toàn mô hình và giá trị nhiệt độ, áp suất cục bộ tại điểm kiểm tra có tọa độ $(x = 0{,}028\text{ m}; y = 0{,}00095\text{ m}; z = -0{,}0035\text{ m})$.

Kết luận

Đồ án tốt nghiệp của nhóm tác giả đã xây dựng hoàn chỉnh mô hình số 3D khảo sát đặc tính nhiệt - thủy lực của thiết bị trao đổi nhiệt dạng ống có cánh ứng dụng trong môi trường biển bằng phần mềm ANSYS 19.2. Nghiên cứu đã xác định được bộ thông số vận hành và hình học tối ưu cục bộ gồm vận tốc đầu vào $1{,}5\text{ m/s}$, nhiệt độ đầu vào $300\text{ K}$ và đường kính ống $7\text{ mm}$, đồng thời chứng minh hiệu quả của giải pháp cải tiến cánh tản nhiệt bằng các cánh nhỏ phụ dài $4\text{ mm}$ nghiêng một góc $15^\circ$. Kết quả nghiên cứu cung cấp cơ sở dữ liệu mô phỏng tin cậy, đóng góp giải pháp thiết kế phục vụ tối ưu hóa các thiết bị trao đổi nhiệt trong kỹ thuật nhiệt và công nghiệp tàu thủy.