Tổng quan nghiên cứu

Trong kỹ thuật thủy động lực học và công nghệ hàng hải, lực cản ma sát nhớt chiếm tới 90% tổng lực cản tác dụng lên các vật thể thon dài di chuyển trong môi trường chất lỏng. Sự tổn hao năng lượng lớn này làm suy giảm đáng kể vận tốc hành trình và gia tăng chi phí nhiên liệu của tàu ngầm, thân tàu thủy và chân vịt. Nhằm giải quyết triệt để rào cản vật lý trên, công nghệ tạo dòng chảy có khoang khí hoặc hơi (supercavitation) nổi lên như một giải pháp đột phá, cho phép giảm tới 90% lực ma sát bề mặt thông qua việc hình thành lớp đệm khí bao bọc vật thể. Tuy nhiên, cơ chế thủy động của dòng chảy đa pha có chuyển pha vô cùng phức tạp, liên quan đến tương tác xoáy rối, tách pha và hiện tượng sụp đổ khoang khí sinh xung áp suất cục bộ lên tới hàng trăm bar.

Luận văn thạc sĩ chuyên ngành Cơ kỹ thuật của tác giả Nguyễn Quang Thái, thực hiện dưới sự hướng dẫn của Giáo sư Nguyễn Tất Thắng tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội phối hợp cùng Viện Cơ học giai đoạn 2014–2018, tập trung làm chủ và ứng dụng bộ phần mềm mã nguồn mở OpenFOAM để tính toán động lực học dòng chảy không có và có chuyển pha. Nghiên cứu thực hiện mô phỏng hai bài toán trọng tâm: dòng chảy đa pha nước - không khí quanh vật thể xâm nhập mặt nước ở dải vận tốc 1,0 m/s đến 3,0 m/s và dòng chuyển động nhanh trong lòng chất lỏng sinh khoang khí tự nhiên ở chuẩn số khoang bằng 1,0. Công trình mang lại giải pháp mô phỏng số chi phí bản quyền 0 đồng, đạt độ tin cậy cao khi đối sánh với dữ liệu thực nghiệm quang học, tạo nền tảng vững chắc cho việc thiết kế các phương tiện thủy động lực học hiệu năng cao tại Việt Nam.

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 ứng dụng hệ phương trình Navier-Stokes trung bình cho dòng chảy đa pha chịu tác động của trọng trường 9,81 m/s² và khối lượng riêng chất lỏng xấp xỉ 1000 kg/m³. Mặt phân cách lỏng - khí được theo dõi chính xác bằng phương pháp thể tích hữu thể Volume of Fluid (VOF), trong đó tỷ phần thể tích pha được bảo toàn qua phương trình truyền tải. Đối với các cấu trúc xoáy rối ở vận tốc cao, mô hình mô phỏng xoáy lớn Large Eddy Simulation (LES) được áp dụng để phân giải trực tiếp các xoáy lớn mang năng lượng và mô hình hóa các xoáy quy mô nhỏ dưới lưới. Quá trình bốc hơi và sụp đổ bọt khí tuân theo phương trình động lực học bọt Rayleigh-Lamb rút gọn từ phương trình Rayleigh-Plesset.

Các chuẩn số thủy động lực học cốt lõi được áp dụng gồm: số khoang sigma xác định mức độ bốc hơi của dòng chảy; hệ số áp suất Cp với giá trị áp suất tối thiểu tại vùng bao phủ khoang đạt mức trừ sigma; chuẩn số Reynolds Re thể hiện tỷ số lực quán tính trên lực nhớt; chuẩn số Froude Fr phản ánh tương tác trọng trường; và hệ số cấp khí Cq điều khiển kích thước khoang khí nhân tạo. Tổng lực cản tác dụng lên vật thể được phân tách thành lực cản chênh áp và lực cản ma sát, trong đó lớp màng khí giúp triệt tiêu phần lớn lực ma sát tiếp xúc.

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

Nghiên cứu kết hợp phương pháp tổng hợp phân tích tài liệu và phương pháp thí nghiệm số chuyên sâu. Nguồn dữ liệu kiểm chứng thực nghiệm được kế thừa từ các công trình đo trường vận tốc hạt PIV (Particle Image Velocimetry) của Wornik công bố năm 2013 và hệ thống quan sát quang học sử dụng camera tốc độ cao từ 60 fps đến 1000 fps tại Viện Cơ học vận hành giai đoạn 2014–2016.

Lý do lựa chọn OpenFOAM viết trên nền tảng C++ hướng đối tượng là tính mở hoàn toàn của mã nguồn, cho phép tùy biến các phương trình vi phân và điều kiện biên phức tạp mà không bị khóa tính năng như các phần mềm thương mại. Quá trình tính toán sử dụng thuật toán liên kết vận tốc - áp suất PIMPLE kết hợp giữa PISO và SIMPLE. Miền tính toán 3D được số hóa với cỡ mẫu hơn 1.200.000 ô lưới lục giác có cấu trúc bằng công cụ blockMesh và snappyHexMesh, tập trung làm mịn cục bộ tại lớp biên đầu dính ướt đường kính 5 mm và vùng đuôi khoang khí nhằm kiểm soát sai số số học trong suốt chu kỳ tính toán.

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

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

Thứ nhất, mô phỏng quá trình vật thể dạng hình trụ và thanh dài đầu phẳng xâm nhập mặt nước ở các mức vận tốc 1,0 m/s, 2,0 m/s và 3,0 m/s đã tái hiện trọn vẹn các giai đoạn: va chạm mặt thoáng, giãn nở túi khí, bóp nghẹt tại bề mặt và đóng kín hoàn toàn trong nước. Kích thước hình học khoang khí mô phỏng tương thích hoàn hảo với ảnh chụp thực nghiệm ở tốc độ 1000 fps với sai lệch thể tích dưới 6%.

Thứ hai, trong bài toán vật thể hình cầu chuyển động ngập hoàn toàn ở số khoang sigma bằng 1,0 tại thời điểm t = 0,36 giây, OpenFOAM mô tả chính xác sự hình thành màng hơi tự nhiên khi áp suất cục bộ hạ xuống dưới áp suất hơi bão hòa. Diện tích tiếp xúc chất lỏng giảm hơn 85%, giúp lực cản ma sát giảm tới 90%, đồng thời làm thay đổi tỷ trọng lực cản áp suất chiếm từ 90% đến 95% tổng lực cản tác dụng.

Thứ ba, nghiên cứu ghi nhận sự xuất hiện của dòng chảy ngược re-entrant jet tại đuôi khoang khí, làm mất ổn định màng hơi và sinh ra các xung áp suất cực đại vượt 100 bar khi bọt khí sụp đổ.

Thứ tư, đối sánh số liệu giữa OpenFOAM và các phần mềm ANSYS Fluent hay UNCEL Code chứng minh OpenFOAM đạt độ chính xác tương đồng với sai số phân bố áp suất bề mặt dưới 5%, đồng thời tiết kiệm toàn bộ chi phí mua bản quyền thương mại.

Thảo luận kết quả

Nguyên nhân vật lý của việc giảm lực cản bắt nguồn từ việc độ nhớt động lực học của pha hơi và khí nhỏ hơn hàng chục lần so với độ nhớt động lực học của nước (khoảng 0,001 kg/ms), làm tiêu giảm triệt để gradient vận tốc tại lớp biên sát tường. Hình dạng đầu dính ướt dạng phẳng, côn hay thuôn quyết định trực tiếp đến đường kính cực đại Dmax và chiều dài cực đại Lmax của khoang khí, hoàn toàn phù hợp với công thức giải tích Garabedian công bố năm 1956. Khi thể hiện kết quả qua đồ thị hệ số áp suất Cp theo chiều dài vật thể, đường cong áp suất duy trì một đoạn phẳng nằm ngang tương ứng với vùng khoang khí có giá trị bằng trừ sigma trước khi tăng vọt tại vùng tái bám đuôi.

Các dữ liệu phân bố pha và trường vận tốc dòng chảy được hiển thị trực quan thông qua biểu đồ vector vận tốc PIV và các lát cắt đẳng trị tỷ phần thể tích. So với các thí nghiệm ống thủy động tại Ukraine và Trung Quốc, việc ứng dụng mô hình LES trên OpenFOAM đã vượt qua rào cản quan sát phẳng 2D của camera thông thường, cung cấp bức tranh không gian 3 chiều đầy đủ về cấu trúc xoáy rối và cơ chế bốc hơi cục bộ quanh thân vật thể.

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

  • Chuẩn hóa thư viện mã nguồn mở mô phỏng siêu khoang: Đề nghị Viện Cơ học phối hợp các trường đại học kỹ thuật xây dựng và đóng gói bộ solver chuyên dụng trên OpenFOAM trước quý 4 năm 2026, tích hợp mô hình nén được cho dải vận tốc trên 1400 m/s nhằm rút ngắn 20% thời gian thiết lập lưới và hội tụ bài toán.
  • Nâng cấp trang thiết bị hầm thủy động thực nghiệm: Đề xuất các cơ quan quản lý khoa học đầu tư hệ thống camera quang học tốc độ cao trên 5000 fps và hệ thống cảm biến áp suất dải đo đến 150 bar trong giai đoạn 2026–2027, nâng độ chính xác kiểm chứng mô hình số lên trên 98%.
  • Phát triển giao diện đồ họa người dùng chuyên biệt: Khuyến nghị các nhóm nghiên cứu công nghệ tính toán phát triển giao diện đồ họa GUI tích hợp trên nền Helyx-OS hoặc SimFlow trong vòng 18 tháng tới, giúp gia tăng 35% năng suất làm việc cho các kỹ sư phân tích thủy động lực học.
  • Ứng dụng công nghệ cấp khí giảm cản thân tàu: Đề xuất các viện thiết kế tàu thủy và tập đoàn đóng tàu triển khai thử nghiệm hệ thống thổi khí tạo đệm khoang khí nhân tạo dưới đáy tàu vận tải thương mại, hướng tới mục tiêu tiết kiệm 8% đến 12% mức tiêu hao nhiên liệu hành trình vào năm 2028.

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

  • Kỹ sư thiết kế khí động và thủy động lực học tàu thủy: Nắm vững quy luật hình thành dòng siêu khoang để tối ưu hóa biên dạng cánh chân vịt tốc độ cao và thiết kế thân vỏ tàu ngầm, giảm thiểu tổn thất ma sát và ngăn ngừa xâm thực phá hủy cơ học.
  • Giảng viên và học viên cao học chuyên ngành Cơ kỹ thuật: Sử dụng công trình như tài liệu học thuật tiêu chuẩn về kỹ thuật giải phương trình đạo hàm riêng bằng C++ trên OpenFOAM, nắm bắt phương pháp chia lưới blockMesh và snappyHexMesh cho dòng đa pha.
  • Chuyên gia phân tích mô phỏng CFD trong công nghiệp: Làm chủ quy trình tiền xử lý, thiết lập thuật toán PIMPLE và khai báo điều kiện biên phức tạp trên nền tảng mã nguồn mở, giúp doanh nghiệp tiết kiệm hàng trăm triệu đồng chi phí bản quyền mỗi năm.
  • Nhà phát triển thiết bị lặn tự hành dưới nước: Ứng dụng nguyên lý tạo khoang khí nhân tạo qua đầu cấp khí để đột phá giới hạn vận tốc di chuyển ngầm của các phương tiện tự hành khảo sát biển sâu.

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

  • Cơ chế nào giúp dòng chảy siêu khoang làm giảm lực cản chuyển động dưới nước? Khi vật thể di chuyển ở vận tốc cao hoặc được cấp khí chủ động, một khoang khí/hơi bao bọc toàn bộ thân vật thể được hình thành. Do độ nhớt chất khí thấp hơn nước hàng chục lần, diện tích tiếp xúc ướt giảm mạnh giúp lực cản ma sát giảm tới 90%, nâng cao vượt bậc tốc độ hành trình.

  • Tại sao OpenFOAM được ưu tiên sử dụng thay thế các phần mềm thương mại đắt tiền? OpenFOAM là phần mềm mã nguồn mở hoàn toàn miễn phí theo giấy phép GPL, giúp tiết kiệm 100% kinh phí bản quyền. Cấu trúc C++ linh hoạt cho phép người dùng can thiệp sâu vào mã nguồn, tùy biến phương trình chuyển pha và chạy song song hiệu năng cao không giới hạn số lõi vi xử lý.

  • Sự khác biệt giữa khoang khí tự nhiên và khoang khí nhân tạo là gì? Khoang khí tự nhiên hình thành do vận tốc chuyển động lớn làm áp suất chất lỏng giảm sâu xuống dưới áp suất bão hòa gây sôi hóa hơi ở nhiệt độ thường. Ngược lại, khoang khí nhân tạo được tạo ra bằng cách chủ động bơm khí không ngưng từ đầu vật thể, duy trì lớp đệm khí ngay cả khi vận tốc dòng chảy chưa đạt ngưỡng hóa hơi.

  • Vai trò của phương pháp VOF và mô hình LES trong mô phỏng dòng đa pha là gì? Phương pháp VOF theo dõi chính xác sự biến thiên của mặt phân cách lỏng - khí qua tỷ phần thể tích. Kết hợp cùng mô hình LES giải quyết các xoáy rối quy mô lớn, hệ thống tái hiện chân thực hiện tượng phân tách pha, cuộn xoáy và bóp nghẹt túi khí với sai số hình học dưới 6%.

  • Hiện tượng đóng kín và sụp đổ khoang khí gây ra những tác động cơ học nào? Tại vùng đuôi khoang khí, quá trình sụp đổ bọt khí sinh ra dòng chảy ngược và sóng xung kích áp suất cục bộ lên tới hàng trăm bar. Áp lực cực lớn này có thể gây xói mòn, nứt vỡ kết cấu kim loại hoặc được ứng dụng có kiểm soát trong cơ cấu kích nổ của vũ khí ngầm.

Kết luận

  • Làm chủ hoàn chỉnh quy trình thiết lập bài toán, phân rã lưới và mô phỏng số dòng chảy đa pha chuyển pha trên nền tảng mã nguồn mở OpenFOAM.
  • Chứng minh và định lượng thành công khả năng cắt giảm tới 90% lực ma sát nhớt nhờ lớp đệm khoang khí bao quanh vật thể chuyển động tốc độ cao.
  • Kiểm chứng độ tin cậy của mô hình số với sai số dưới 5% khi đối chiếu cùng dữ liệu thực nghiệm đo trường vận tốc PIV và ảnh chụp camera tốc độ cao.
  • Đặt nền móng kỹ thuật và công nghệ tính toán cho việc phát triển các phương tiện thủy động lực học hiệu năng cao phục vụ ngành kinh tế biển và quốc phòng tại Việt Nam.
  • Trong giai đoạn 2026–2028, các nghiên cứu tiếp theo cần mở rộng tích hợp tính nén được của chất lỏng ở dải vận tốc trên 1400 m/s; quý bạn đọc và kỹ sư hãy tham khảo ngay luận văn này để ứng dụng OpenFOAM vào các dự án nghiên cứu thủy khí động lực học chuyên sâu.