Tổng quan nghiên cứu

Trong kỹ thuật thủy khí động lực học và thiết kế phương tiện ngầm, lực cản ma sát giữa bề mặt vỏ kim loại với môi trường nước chiếm từ 70% đến 90% tổng lực cản tác động lên thân phương tiện có chiều dài lớn. Việc ứng dụng dòng chảy có khoang khí hoặc hơi bao quanh thân vật thể mang lại tiềm năng giảm lực cản ma sát tới hơn 90%, tạo bước đột phá trong việc tối ưu hóa hiệu suất năng lượng và gia tăng vận tốc di chuyển. Tuy nhiên, bản chất phức tạp của dòng chảy đa pha, quá trình hóa hơi khi áp suất giảm xuống áp suất hơi bão hòa, cùng hiện tượng đóng kín khoang sinh ra xung áp suất lên tới hàng trăm bar đặt ra bài toán khó cho việc tính toán lý thuyết và thực nghiệm.

Luận văn thạc sĩ chuyên ngành Cơ kỹ thuật mang tên "Nghiên cứu, ứng dụng bộ chương trình OpenFOAM trong tính toán động lực học dòng chảy không có/có chuyển pha" do học viên Nguyễn Quang Thái thực hiện tại Trường Đại học Công nghệ, Đại học Quốc gia Hà Nội phối hợp với Viện Cơ học vào năm 2018. Mục tiêu trọng tâm của nghiên cứu là làm chủ bộ công cụ mã nguồn mở OpenFOAM bao gồm hơn 100 thư viện C++ và khoảng 250 ứng dụng chuyên dụng để thiết lập mô hình tính toán số cho dòng chảy hai pha không có và có chuyển pha.

Phạm vi nghiên cứu tập trung vào hai bài toán cốt lõi: mô phỏng vật thể xâm nhập qua mặt thoáng vào nước ở dải vận tốc từ 1,47 m/s đến 3,0 m/s và chuyển động của quả cầu tạo khoang hơi ở các số khoang từ 1,0 đến 2,5. Công trình mang ý nghĩa thực tiễn to lớn khi giúp các đơn vị nghiên cứu tại Việt Nam làm chủ công nghệ mô phỏng thủy động lực học chuyên sâu, tiết kiệm 100% chi phí mua bản quyền phần mềm thương mại đắt đỏ và rút ngắn thời gian thử nghiệm thực tế.

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

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng cơ học chất lưu liên tục và nhiệt động lực học dòng chảy đa pha với hai hệ thống lý thuyết chủ đạo:

Thứ nhất là lý thuyết dòng chảy nhiều pha và động lực học bọt khí dựa trên phương trình Rayleigh - Lamb. Khi chất lỏng di chuyển với vận tốc lớn, áp suất cục bộ giảm xuống dưới áp suất hơi bão hòa tại nhiệt độ môi trường, kích hoạt quá trình chuyển pha từ lỏng sang hơi. Phương trình rút gọn mô tả tốc độ phát triển bán kính bọt hơi theo thời gian dưới tác động của độ chênh lệch giữa áp suất bão hòa và áp suất cục bộ trong chất lỏng.

Thứ hai là lý thuyết mô hình hóa dòng chảy rối thông qua phương pháp mô phỏng xoáy lớn kết hợp với mô hình thể tích hữu hạn thể tích pha để theo dõi mặt phân cách giữa pha lỏng và pha khí. Dòng chảy được giải quyết thông qua hệ phương trình Navier - Stokes kết hợp thuật toán PIMPLE, liên kết chặt chẽ giữa hai thuật toán kinh điển là PISO và SIMPLE.

Nghiên cứu chuẩn hóa năm khái niệm và tham số không thứ nguyên then chốt bao gồm: số khoang mô tả mức độ phát triển của túi khí, hệ số áp suất đặc trưng cho phân bố áp lực bề mặt, hệ số lực cản tổng hợp, số Reynolds và hệ số cấp khí khi tạo khoang nhân tạo.

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

Nguồn dữ liệu của luận văn được tổng hợp từ hai kênh chính: dữ liệu mô phỏng số 3D trên nền tảng OpenFOAM và dữ liệu thực nghiệm đối chứng từ hệ thống ống thủy động chế tạo tại Viện Cơ học năm 2016 cùng hệ thống camera tốc độ cao ghi hình từ 60 fps đến 1000 fps.

Cỡ mẫu khảo sát mô phỏng gồm 30 trường hợp tính toán với các cấu hình lưới 3D từ 80.000 đến 650.000 phần tử khối. Phương pháp chọn mẫu áp dụng theo kỹ thuật lấy mẫu phân tầng dựa trên các dải vận tốc chuyển động và biên dạng đầu dính ướt của vật thể (đầu phẳng, đầu bán cầu và đầu côn có đường kính 5 mm). Lý do lựa chọn phương pháp phân tích thể tích hữu hạn với thuật toán PIMPLE là khả năng xử lý xuất sắc các bài toán quá độ không ổn định theo thời gian, cho phép bước thời gian tính toán siêu nhỏ ở mức $10^{-5}$ giây mà vẫn đảm bảo độ hội tụ và tính bảo toàn khối lượng nghiêm ngặt.

Quy trình nghiên cứu được triển khai đồng bộ từ việc tiền xử lý lưới với blockMesh và snappyHexMesh, cấu hình điều kiện biên tại thư mục 0, thiết lập các hằng số nhiệt động lực học trong thư mục constant và xuất kết quả phân tích trường dòng chảy qua phần mềm ParaView trong toàn bộ giai đoạn từ năm 2016 đến 2018.

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

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

Quá trình mô phỏng số trên OpenFOAM đã mang lại bốn phát hiện quan trọng có giá trị khoa học và kỹ thuật cao:

Thứ nhất, đối với bài toán kiểm chứng dòng chảy tầng qua ống thay đổi tiết diện chiều dài 1,5 m và đoạn thu hẹp 0,5 m, bộ giải icoFoam tái hiện chính xác trường vận tốc đầu vào $U_x = 1,0$ m/s và áp suất đầu ra 0 bar. Kết quả cho thấy dòng chất lỏng phân tách đều và hội tụ chính xác qua khe hẹp mà không xảy ra sai số tích lũy biên.

Thứ hai, trong bài toán vật thể hình học xâm nhập nước qua mặt thoáng, mô hình tính toán phản ánh chân thực từng giai đoạn từ va chạm ban đầu, hình thành túi khí tự nhiên đến khi khoang đóng kín ở các dải vận tốc 1,47 m/s, 2,0 m/s và 3,0 m/s. Độ sai lệch về kích thước chiều dài khoang khí giữa mô phỏng OpenFOAM và ảnh chụp thực nghiệm 1000 fps được kiểm soát ở mức dưới 4,5%.

Thứ ba, sự hình thành khoang khí/hơi quanh quả cầu ở số khoang $\sigma = 1,0$ cho thấy lớp khí bao bọc làm giảm diện tích tiếp xúc trực tiếp giữa thân vật thể và nước, dẫn đến lực cản ma sát giảm tới 88% so với khi di chuyển trong dòng chảy đơn pha.

Thứ tư, thuật toán PIMPLE phối hợp phương pháp thể tích pha đã mô tả xuất sắc quá trình biến thiên theo thời gian từ $t = 0,002$ giây đến $t = 0,57$ giây, chỉ ra sự xuất hiện của dòng xoáy ngược tại đuôi khoang khí và hiện tượng áp suất đỉnh khi bọt khí sụp đổ.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự hình thành khoang hơi là hiện tượng giảm áp cục bộ tại lớp biên rối của vật thể khi vận tốc tương đối tăng cao, đưa áp suất chất lỏng về mức áp suất hơi bão hòa. So sánh kết quả tính toán giữa OpenFOAM và phần mềm thương mại ANSYS Fluent cho thấy độ tương đồng đạt trên 96% về trường áp suất $C_p$ và phân bố đường dòng. Tuy nhiên, OpenFOAM thể hiện ưu thế vượt trội nhờ khả năng can thiệp trực tiếp vào mã nguồn C++, cho phép người dùng tùy biến các mô hình chuyển pha phức tạp mà không tốn chi phí bản quyền lên đến hàng chục nghìn USD mỗi năm.

Dữ liệu nghiên cứu được biểu diễn trực quan qua các bảng đối chiếu tính năng kỹ thuật và biểu đồ phân bố hệ số áp suất $C_p$ theo tỉ số tắc nghẽn của ống thủy động. Biểu đồ đường dòng chỉ ra rằng các vật thể có đầu dính ướt vát góc tạo ra khoang khí ổn định hơn và triệt tiêu được 60% dòng chảy ngược ở phần đuôi so với dạng đầu phẳng truyền thống.

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

Dựa trên các kết quả đạt được từ luận văn, bốn giải pháp trọng tâm được khuyến nghị nhằm thúc đẩy ứng dụng công nghệ mô phỏng CFD tại Việt Nam:

Thứ nhất, chuẩn hóa và đóng gói bộ thư viện OpenFOAM chuyên dụng cho bài toán dòng chảy nhiều pha và siêu khoang. Bộ môn Cơ kỹ thuật phối hợp cùng các viện nghiên cứu chuyên ngành cần xây dựng bộ template tính toán tự động trong giai đoạn 2026-2027, đặt mục tiêu rút ngắn 50% thời gian thiết lập lưới và điều kiện biên cho các nghiên cứu sinh mới tiếp cận.

Thứ hai, nâng cấp trang thiết bị thực nghiệm kiểm chứng thủy động lực học. Viện Hàn lâm Khoa học và Công nghệ Việt Nam cần đầu tư bổ sung hệ thống camera ghi hình quang học đạt tốc độ trên 2000 fps cùng cảm biến đo trường vận tốc PIV laser 3D trong vòng 18 tháng tới, nhằm nâng cao độ chính xác đo đạc thực nghiệm lên 99%.

Thứ ba, ứng dụng mô hình khoang khí nhân tạo vào thiết kế phương tiện ngầm và chân vịt tàu thủy. Các nhà máy đóng tàu và cơ sở công nghiệp hàng hải cần áp dụng giải pháp tạo túi khí dưới đáy thân tàu và cánh chân vịt góc vát từ năm 2027, hướng tới chỉ tiêu cắt giảm từ 15% đến 25% lượng tiêu thụ nhiên liệu hành trình.

Thứ tư, mở rộng thuật toán mô phỏng cho dòng nén được siêu âm dưới nước. Nhóm nghiên cứu cần tiếp tục phát triển mã nguồn OpenFOAM để xử lý bài toán chuyển động vượt vận tốc âm trong nước (trên 1400 m/s) trong kế hoạch 24 tháng tiếp theo, phục vụ các đề tài nghiên cứu quốc phòng công nghệ cao.

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

Luận văn là tài liệu tham khảo giá trị cho bốn nhóm độc giả chuyên môn:

Nhóm thứ nhất là giảng viên, học viên cao học và sinh viên ngành Cơ kỹ thuật, Thủy khí động lực học, Kỹ thuật Hàng hải. Tài liệu cung cấp hướng dẫn tường tận về cấu trúc case thư mục, lập trình C++ trong OpenFOAM và phương pháp luận thiết lập mô hình tính toán số chuẩn mực.

Nhóm thứ hai là các kỹ sư R&D trong ngành công nghiệp đóng tàu và cơ khí thủy lực. Luận văn cung cấp phương pháp tính toán lực cản, thiết kế chân vịt siêu khoang và giải pháp bọc túi khí giảm ma sát cho thân tàu vỏ dài trong các dự án đóng mới phương tiện thủy.

Nhóm thứ ba là các nhà khoa học và chuyên gia tại các viện nghiên cứu kỹ thuật quân sự. Công trình mang lại các cơ sở dữ liệu số quan trọng về hiện tượng xâm nhập nước của vật thể và chuyển động tốc độ cao của khí tài ngầm dạng ngư lôi Shkval.

Nhóm thứ tư là các lập trình viên tính toán hiệu năng cao và kỹ sư phát triển phần mềm mô phỏng. Nghiên cứu minh họa chi tiết cách can thiệp mã nguồn hướng đối tượng, chuyển đổi định dạng lưới và tối ưu hóa xử lý song song trên hệ điều hành Linux.

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

Hiện tượng tạo khoang khí hoặc hơi giúp giảm lực cản cho vật thể di chuyển dưới nước như thế nào? Khi khoang khí bao bọc xung quanh thân vật thể, diện tích tiếp xúc trực tiếp giữa vỏ kim loại và chất lỏng giảm đi đáng kể. Do độ nhớt của khí và hơi nhỏ hơn hàng trăm lần so với nước, lực cản ma sát bề mặt có thể giảm tới 90%, giúp phương tiện đạt vận tốc vượt trội với cùng công suất động cơ.

Ưu điểm nổi bật nhất của OpenFOAM so với phần mềm thương mại ANSYS Fluent trong nghiên cứu là gì? OpenFOAM là phần mềm mã nguồn mở hoàn toàn miễn phí, cho phép người dùng tự do can thiệp và chỉnh sửa mã nguồn C++ để phát triển các mô hình vật lý mới. Ngược lại, Fluent là phần mềm mã nguồn đóng, chi phí bản quyền rất cao và hạn chế tối đa khả năng tùy biến thuật toán giải phương trình vi phân.

Luận văn đã sử dụng phương pháp nào để mô phỏng dòng chảy chuyển pha và dòng rối? Nghiên cứu kết hợp phương pháp thể tích pha để theo dõi chính xác ranh giới mặt phân cách lỏng - hơi, phương pháp mô phỏng xoáy lớn để nắm bắt các cấu trúc xoáy rối phức tạp và áp dụng phương trình Rayleigh - Lamb rút gọn để mô tả quá trình chuyển pha nhiệt động học.

Làm thế nào để xác thực độ chính xác của kết quả mô phỏng số trong luận văn? Tác giả đã tiến hành đối chứng trực tiếp các kết quả mô phỏng với dữ liệu thực nghiệm thu được từ hệ thống ống thủy động tại Viện Cơ học và các hình ảnh ghi nhận từ camera tốc độ cao 1000 fps, ghi nhận độ sai lệch về kích thước khoang chỉ dưới 5%.

Khoang khí tự nhiên và khoang khí nhân tạo khác nhau ở điểm nào? Khoang khí tự nhiên hình thành do áp suất chất lỏng tự giảm xuống áp suất hơi bão hòa khi vật thể di chuyển ở vận tốc rất cao. Khoang khí nhân tạo được tạo ra bằng cách chủ động bơm khí không ngưng tụ từ các lỗ nhỏ trên đầu vật thể, giúp hình thành túi khí ngay ở dải vận tốc thấp hơn.

Kết luận

Luận văn thạc sĩ của tác giả Nguyễn Quang Thái đã hoàn thành xuất sắc các mục tiêu nghiên cứu với những dấu ấn chuyên môn nổi bật:

  • Làm chủ hoàn toàn cấu trúc phần mềm mã nguồn mở OpenFOAM với hơn 100 thư viện lập trình C++ và 250 công cụ giải số phục vụ thủy động lực học.
  • Xây dựng thành công mô hình mô phỏng dòng chảy hai pha không có chuyển pha cho bài toán vật thể xâm nhập nước ở vận tốc lên tới 3,0 m/s với độ chính xác cao.
  • Thiết lập hoàn chỉnh mô hình dòng chảy có chuyển pha VOF kết hợp mô phỏng xoáy lớn, chứng minh khả năng giảm 88% lực cản ma sát ở số khoang 1,0.
  • Đối chứng và kiểm chứng tin cậy kết quả mô phỏng số với dữ liệu thực nghiệm từ hệ thống ống thủy động và camera tốc độ cao 1000 fps tại Viện Cơ học.
  • Đặt nền móng kỹ thuật vững chắc cho việc phát triển các bộ giải CFD mã nguồn mở nội địa, loại bỏ hoàn toàn sự phụ thuộc vào các phần mềm thương mại đắt tiền.

Trong giai đoạn tiếp theo, hướng phát triển trọng tâm là mở rộng thuật toán sang bài toán dòng nén được siêu âm và tích hợp tối ưu hóa tự động hình học thân tàu. Quý độc giả, chuyên gia và học viên quan tâm có thể khai thác các mô hình mẫu trong luận văn để ứng dụng vào công tác giảng dạy, nghiên cứu khoa học và phát triển các sản phẩm công nghiệp thực tiễn ngay hôm nay.