Tổng quan nghiên cứu

Quạt ly tâm công suất 5.5 kW là thiết bị cơ học giữ vai trò trọng yếu trong hệ thống thông gió công nghiệp, xử lý khí thải và điều hòa không khí nhà xưởng, nơi hệ thống động cơ quạt thường chiếm từ 60% đến 70% tổng điện năng tiêu thụ toàn trạm. Tại Việt Nam, quy trình thiết kế bánh công tác quạt ly tâm hiện nay phần lớn vẫn phụ thuộc vào kinh nghiệm truyền thống hoặc gia công phỏng theo các mẫu có sẵn của nước ngoài. Cách tiếp cận này khiến thiết bị thường chịu tổn thất khí động học lớn, vận hành dưới điểm làm việc tối ưu và tiêu tốn nhiều năng lượng.

Mục tiêu trọng tâm của nghiên cứu là xây dựng và lựa chọn mô hình mô phỏng số 3D chính xác đặc tính dòng chảy của quạt ly tâm công suất 5.5 kW trên nền tảng phần mềm ANSYS CFX, đánh giá toàn diện ảnh hưởng của cấu trúc lưới và các mô hình dòng rối. Đồng thời, đề tài tiến hành thiết lập chu trình tự động tối ưu hóa biên dạng hình học cánh bằng giải thuật tiến hóa thông qua việc kết nối hai phần mềm OpenFOAM và DAKOTA, hướng đến mục tiêu giảm thiểu hệ số lực cản khí động học.

Phạm vi nghiên cứu tập trung vào mẫu quạt ly tâm công suất 5.5 kW của Công ty TNHH Kỹ thuật Công nghiệp Đồng Tâm kết hợp các bộ dữ liệu khí động học chuẩn quốc tế, được thực hiện trong thời gian từ tháng 09/2017 đến tháng 06/2018 tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh. Việc ứng dụng chu trình tính toán này giúp rút ngắn khoảng 40% thời gian phát triển thiết bị, cắt giảm hơn 30% chi phí thử nghiệm thực tế và cung cấp công cụ tự động hóa đắc lực nhằm nâng cao hiệu suất động lực học cho ngành chế tạo máy thủy khí trong nước.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên hai nền tảng lý thuyết cốt lõi trong cơ học chất lỏng và toán học tối ưu. Trước hết là hệ phương trình vi phân Navier-Stokes 3 chiều mô tả các định luật bảo toàn khối lượng (phương trình liên tục), bảo toàn động lượng và bảo toàn năng lượng cho dòng lưu chất không nén được. Để giải quyết các chuyển động xoáy phức tạp, lý thuyết dòng chảy rối được áp dụng thông qua hai mô hình hai phương trình tiêu biểu: mô hình k - epsilon (tập trung vào động năng rối và tốc độ tiêu tán rối) và mô hình k - omega (đặc biệt là biến thể SST của Menter, cho phép kiểm soát chính xác sự phân tách dòng tại lớp biên sát bề mặt kim loại).

Về mặt tối ưu hóa, đề tài ứng dụng Lý thuyết Thuật toán tiến hóa (Evolutionary Algorithm - EA) kết hợp giải thuật di truyền đơn mục tiêu (SOGA) dựa trên nguyên lý chọn lọc tự nhiên của Darwin. Mô hình toán học kiểm soát chặt chẽ 4 cơ chế tiến hóa sinh học gồm: sinh sản, chọn lọc tự nhiên, lai ghép chéo và đột biến ngẫu nhiên nhằm tìm kiếm lời giải toàn cục. Bên cạnh đó, các khái niệm chuyên ngành then chốt được tham số hóa bao gồm: hệ số lực cản, khoảng cách không thứ nguyên tới bề mặt vách để đánh giá chất lượng lớp biên, phương pháp hệ quy chiếu quay phân tách miền tĩnh và miền xoay, cùng kỹ thuật đường cong B-spline để nội suy độ cong hình học cánh quạt qua các điểm điều khiển linh hoạt.

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

Dữ liệu đầu vào của đề tài bao gồm mô hình 3D thực tế của quạt ly tâm 5.5 kW do Công ty Đồng Tâm cung cấp và các bộ tham số biên dạng khí động học chuẩn từ các công bố quốc tế của Manuel J. Garcia và Lucio Cardillo. Về phương pháp phân tích số, phương pháp thể tích hữu hạn (Finite Volume Method - FVM) được lựa chọn thay vì sai phân hữu hạn hay phần tử hữu hạn nhờ khả năng bảo toàn vật lý tuyệt đối trên từng ô lưới và tối ưu hóa bộ nhớ xử lý khi giải các bài toán dòng chảy rối phức tạp.

Nghiên cứu tiến hành khảo sát cỡ mẫu mô phỏng trên 2 bộ cấu trúc lưới 3D chính: lưới tam giác có cấu trúc và lưới không cấu trúc tứ diện với mật độ biến thiên từ 50.000 phần tử ở mức thô đến hơn 3.200.000 phần tử ở mức mịn để kiểm tra tính độc lập của lưới. Phương pháp lấy mẫu trong bài toán tối ưu sử dụng quần thể khởi tạo gồm 10 đến 20 cá thể điểm điều khiển và theo dõi qua 15 đến 50 thế hệ tiến hóa. Bộ giải kết hợp của ANSYS CFX được sử dụng cho bài toán 3D quạt ly tâm nhờ giải thuật đa lưới đại số tốc độ cao, trong khi bộ đôi mã nguồn mở OpenFOAM và DAKOTA được chọn cho bài toán tối ưu biên dạng nhờ khả năng lập trình tự động hóa linh hoạt qua dòng lệnh shell script.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, kết quả đánh giá lưới cho thấy mô hình lưới tam giác có cấu trúc mang lại độ ổn định tính toán vượt trội, kiểm soát khoảng cách không thứ nguyên sát thành trong ngưỡng lý thuyết từ 30 đến 100, đồng thời giảm sai lệch áp suất cục bộ xuống dưới 5% so với mô hình lưới không cấu trúc.

Thứ hai, khi so sánh các mô hình rối, mô hình k - omega thể hiện khả năng nắm bắt chính xác hiện tượng phân tách dòng chảy và các xoáy thứ cấp tại mép hút của cánh bánh công tác, mang lại kết quả lưu lượng đầu ra với sai số chỉ 8.5% so với thực nghiệm, trong khi mô hình k - epsilon có xu hướng đánh giá cao hơn thực tế khoảng 10% đến 12%.

Thứ ba, việc tích hợp OpenFOAM và DAKOTA qua thuật toán di truyền đã tối ưu hóa thành công biên dạng hình học cánh, giúp giảm hệ số lực cản 14.6% so với biên dạng ban đầu, trong khi hệ số lực nâng vẫn duy trì ổn định với độ biến thiên không quá 2%.

Thứ tư, phân bố áp suất và vận tốc qua cánh quạt sau tối ưu thể hiện trường dòng chảy trơn mượt, loại bỏ các điểm xoáy cục bộ tại đầu ra, đảm bảo lưu lượng định mức đạt xấp xỉ 4.88 m³/s ở dải tốc độ quay từ 1329 đến 2900 vòng/phút.

Thảo luận kết quả

Sự suy giảm đáng kể của hệ số lực cản bắt nguồn từ việc thuật toán tiến hóa đã tự động điều chỉnh tọa độ các điểm điều khiển B-spline ở mép dẫn và mép thoát của cánh. Sự thay đổi này làm mượt độ cong của cả bề mặt hút lẫn bề mặt đẩy, ngăn chặn sự hình thành sớm của vùng áp suất nghịch khi dòng khí bị chuyển hướng đột ngột 90 độ từ phương xuyên tâm sang phương hướng kính.

Kết quả này hoàn toàn tương đồng với các nghiên cứu của Manuel J. Garcia về tối ưu cánh khí động và khẳng định tính đúng đắn trong các công bố của Lucio Cardillo trên quạt công nghiệp (với sai số kiểm định đều dưới mức 10%). Hơn thế nữa, quy trình tự động hóa được xây dựng đã giúp tiết kiệm hơn 500 giờ tính toán của bộ xử lý so với các quy trình thử sai thủ công truyền thống.

Để phục vụ công tác thiết kế thực tế, toàn bộ dữ liệu phân tích có thể được biểu diễn trực quan qua ba dạng đồ thị chuẩn:

  • Giản đồ đường đặc tính quạt biểu diễn mối quan hệ phi tuyến giữa lưu lượng thể tích và áp suất toàn phần theo từng cấp tốc độ quay.
  • Đồ thị lịch sử hội tụ của hệ số lực cản thể hiện mức độ tối ưu giảm dần sau mỗi chu kỳ thế hệ trong DAKOTA.
  • Bảng ma trận so sánh các chỉ số khí động học giữa các mô hình rối nhằm hỗ trợ kỹ sư lựa chọn phương pháp phân tích phù hợp cho từng dải công suất.

Đề xuất và khuyến nghị

Thứ nhất, chuẩn hóa quy trình mô phỏng số khí động học: Các doanh nghiệp sản xuất cơ khí máy thủy khí cần thiết lập quy trình kiểm tra lưới định kỳ với giá trị khoảng cách sát thành khống chế dưới 50 và ứng dụng mô hình rối k - omega SST, nhằm kiểm soát sai số thiết kế thủy khí dưới 5% trong vòng 6 tháng tới. Chủ thể thực hiện: Bộ phận Nghiên cứu và Phát triển (R&D) tại các nhà máy chế tạo quạt công nghiệp.

Thứ hai, tích hợp chuỗi công cụ tự động mã nguồn mở OpenFOAM - DAKOTA: Đội ngũ kỹ sư cần xây dựng các tệp kịch bản tự động hóa liên kết giữa phần mềm chia lưới, bộ giải CFD và thuật toán tối ưu hóa di truyền, nhằm mục tiêu giảm tối thiểu 15% hệ số lực cản cho các dòng sản phẩm quạt từ 5.5 kW đến 22 kW trong giai đoạn 2024–2025. Chủ thể thực hiện: Kỹ sư mô phỏng số và chuyên viên thiết kế khí động học.

Thứ ba, đầu tư nâng cấp hạ tầng tính toán song song hiệu năng cao (HPC): Các trường đại học và viện nghiên cứu cần bổ sung hệ thống máy chủ tính toán đa lõi từ 64 đến 128 luồng xử lý nhằm rút ngắn thời gian chạy các vòng lặp tối ưu hóa 3D toàn phần xuống dưới 48 giờ cho mỗi chu kỳ sản phẩm trong quý 3 năm tới. Chủ thể thực hiện: Ban quản lý phòng thí nghiệm tính toán hiệu năng cao và các cơ sở đào tạo sau đại học.

Thứ tư, mở rộng thuật toán sang bài toán tối ưu hóa đa mục tiêu: Nhóm nghiên cứu cần phát triển thuật toán di truyền đa mục tiêu kết hợp mô hình xấp xỉ nhằm đồng thời nâng cao hiệu suất quạt thêm 5% đến 8% và giảm mức độ phát thải tiếng ồn công nghiệp ít nhất 3 dB trong lộ trình 12 tháng tiếp theo. Chủ thể thực hiện: Các nhà khoa học, học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Hàng không.

Đối tượng nên tham khảo luận văn

Thứ nhất, kỹ sư thiết kế và phát triển sản phẩm tại các doanh nghiệp cơ khí: Khai thác phương pháp tham số hóa hình học cánh quạt bằng B-spline và quy trình mô phỏng CFD 3D để thiết kế các dòng quạt ly tâm 5.5 kW tiết kiệm điện năng và đạt chuẩn công nghiệp cao.

Thứ hai, học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Hàng không, Kỹ thuật Cơ khí: Sử dụng tài liệu như một cẩm nang mẫu mực về cách thức lập trình kết nối giữa phần mềm CFD mã nguồn mở OpenFOAM và bộ công cụ tối ưu DAKOTA thông qua giải thuật tiến hóa.

Thứ ba, giảng viên và các nhà nghiên cứu học thuật: Dùng làm tài liệu tham khảo chuyên sâu phục vụ công tác giảng dạy các học phần Cơ học lưu chất tính toán, Động lực học dòng chảy trong máy thủy khí và Các phương pháp tối ưu hóa kỹ thuật.

Thứ tư, chuyên gia tư vấn năng lượng và kỹ sư vận hành hệ thống thông gió công nghiệp: Tham khảo phương pháp phân tích giản đồ đặc tính vận hành để lựa chọn điểm làm việc tối ưu cho quạt ly tâm trong hệ thống hút lọc bụi nhà xưởng, thông gió hầm lò và tăng áp cầu thang tòa nhà cao tầng.

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

Thứ nhất: Tại sao nghiên cứu lại chọn quạt ly tâm công suất 5.5 kW làm đối tượng mô phỏng chính? Mẫu quạt 5.5 kW của Công ty Đồng Tâm là dòng quạt công nghiệp phổ biến trong các hệ thống thông gió và hút lọc bụi. Quy mô công suất này có tính ứng dụng thực tiễn rất cao, cung cấp đầy đủ các thông số kích thước và dữ liệu đo đạc thực nghiệm giúp đối chứng kết quả CFD với độ tin cậy sai số dưới 10%.

Thứ hai: Lợi ích lớn nhất của việc liên kết OpenFOAM và DAKOTA trong tối ưu hóa hình học là gì? OpenFOAM là bộ giải CFD mã nguồn mở mạnh mẽ, còn DAKOTA cung cấp các thuật toán tiến hóa đa dạng. Sự kết hợp này loại bỏ hoàn toàn chi phí bản quyền phần mềm thương mại đắt đỏ, đồng thời tự động hóa hoàn toàn chu trình tính toán, giúp tiết kiệm hàng trăm giờ thao tác thủ công khi tinh chỉnh biên dạng cánh.

Thứ ba: Vì sao mô hình rối k - omega lại cho kết quả vượt trội hơn mô hình k - epsilon trong bài toán này? Mô hình k - omega có độ nhạy rất cao với các hiện tượng biến đổi vận tốc sát vách và dự đoán chính xác điểm phân tách dòng khi có gradient áp suất ngược. Ngược lại, mô hình k - epsilon thường dự đoán điểm tách dòng muộn hơn thực tế, dẫn đến việc ước tính hiệu suất quạt cao hơn bình thường từ 10% trở lên.

Thứ tư: Vai trò cốt lõi của đường cong B-spline trong việc thay đổi biên dạng hình học là gì? Đường cong B-spline cho phép định hình toàn bộ biên dạng cánh mượt mà chỉ thông qua một số lượng nhỏ các điểm điều khiển từ 6 đến 12 điểm. Kỹ thuật này giúp thu hẹp không gian tìm kiếm của thuật toán di truyền, rút ngắn thời gian hội tụ mà vẫn duy trì được độ trơn khí động học cần thiết.

Thứ năm: Sai số giữa mô phỏng số 3D và kết quả thực tế được kiểm chứng ở mức độ nào? Bằng cách xây dựng mô hình lưới chất lượng cao với hơn 3 triệu phần tử và thiết lập điều kiện biên dòng xoay chuẩn xác, sai số về lưu lượng khối lượng và áp suất toàn phần giữa mô phỏng ANSYS CFX và thực nghiệm duy trì ở mức 8.5% đến 9.5%, hoàn toàn đáp ứng tiêu chuẩn kỹ thuật công nghiệp.

Kết luận

  • Xây dựng thành công mô hình mô phỏng số 3D hoàn chỉnh cho quạt ly tâm công nghiệp 5.5 kW trên phần mềm ANSYS CFX với độ tin cậy cao.
  • Chứng minh tính vượt trội của cấu trúc lưới có trật tự và mô hình rối k - omega trong việc tái hiện trường dòng chảy thực tế của máy thủy khí.
  • Thiết lập chu trình tối ưu hóa tự động hóa hoàn chỉnh dựa trên sự kết hợp giữa OpenFOAM và DAKOTA thông qua giải thuật tiến hóa.
  • Giảm thiểu thành công 14.6% hệ số lực cản của biên dạng cánh quạt nhờ kỹ thuật tham số hóa B-spline, nâng cao hiệu quả khí động học.
  • Định hình phương pháp luận tính toán số chuẩn xác, làm nền tảng mở rộng cho các bài toán tối ưu hóa hình học 3D phức tạp tại Việt Nam.

Đóng góp cốt lõi của luận văn là đã làm chủ quy trình mô phỏng số và tối ưu hóa tự động, mở ra hướng đi mới giúp các doanh nghiệp cơ khí trong nước nâng cao chất lượng sản phẩm mà không cần phụ thuộc vào mẫu ngoại nhập. Trong giai đoạn 6 đến 12 tháng tới, hướng nghiên cứu sẽ tiếp tục hoàn thiện mô hình tối ưu 3D đa mục tiêu kết hợp âm học cho bánh công tác. Quý độc giả, kỹ sư và các nhà nghiên cứu quan tâm có thể khai thác toàn văn tài liệu hoặc liên hệ trực tiếp để chuyển giao và ứng dụng quy trình công nghệ này vào thực tiễn sản xuất.