Báo Cáo Nghiên Cứu: Mô Phỏng Số Dòng Xoáy Xâm Thực Hơi và Khí Không Hòa Tan Tính Đến Hiệu Ứng Nhiệt Động Lực Học

Tóm tắt nghiên cứu (Executive Summary)

  • Câu hỏi nghiên cứu chính (Research Question): Làm thế nào để mô hình hóa và dự đoán chính xác động học của dòng xoáy xâm thực (cavitation) nhiều pha trong môi trường nước khi chịu tác động đồng thời của cấu trúc xoáy dòng chảy, hiệu ứng nhiệt động lực học và sự hiện diện của khí không hòa tan?
  • Phương pháp luận (Methodology Snapshot): Nghiên cứu phát triển một chương trình mô phỏng số 2D dựa trên mô hình chất lưu nhiều pha đồng nhất (Homogeneous Equilibrium/Compressible Model) giải hệ phương trình Navier-Stokes đa pha bằng phương pháp sai phân hữu hạn (FDM) với sơ đồ sai phân ngược hướng cấp 2 Harten-Yee TVD kết hợp kỹ thuật tách thời gian Strang splitting. Đề tài đề xuất mô hình ngưỡng áp suất chuyển pha mới (SVM - Shear Strain & Vorticity Model) tích hợp hàm biến dạng cắt và độ xoáy cục bộ, được lập trình song song hóa OpenMP trên nền tảng FORTRAN/Linux.
  • Kết quả then chốt (Key Findings): Mô hình SVM cải thiện vượt bậc độ chính xác khi dự đoán phân bố hệ số áp suất ($-C_p$) và quá trình hình thành - tách - sụp đổ của bọt khí xâm thực trên cánh thủy lực NACA0015, vật thể đối xứng đầu trụ và vòi phun hội tụ - phân kỳ Abuaf. Nghiên cứu chỉ ra rằng khí không hòa tan thúc đẩy xâm thực xuất hiện sớm hơn về phía đầu mép cánh, trong khi hiệu ứng nhiệt động lực học kìm hãm kích thước khối bọt khí phụ thuộc mật thiết vào tỷ lệ thang thời gian đặc trưng ($\tau/\tau_{th}$).
  • Hàm ý thực tiễn (Implications): Cung cấp công cụ tính toán số tin cậy, chi phí thấp cho việc thiết kế tối ưu, kiểm nghiệm và nâng cao tuổi thọ của các thiết bị thủy lực công nghiệp, chân vịt tàu thủy, hệ thống đẩy hàng không vũ trụ và tua-bin thủy điện.

Bối cảnh và tầm quan trọng (Context & Significance)

1. Hiện trạng tri thức và thách thức công nghệ

Hiện tượng xâm thực (cavitation) xảy ra khi áp suất tĩnh cục bộ trong dòng chất lỏng giảm xuống dưới áp suất hơi bão hòa do sự gia tốc dòng chảy, dẫn đến việc hình thành các bọt khí/hơi. Khi các bọt khí này bị cuốn về vùng có áp suất cao hơn, chúng sụp đổ (implosion) mãnh liệt trong thang thời gian vi giây, sinh ra các vi tia nước (micro-jets) và sóng xung kích có áp suất cục bộ lên tới hàng nghìn megapascal cùng nhiệt độ hàng nghìn Kelvin. Hiện tượng này gây rung động mạnh, tiếng ồn thủy âm nghiêm trọng, làm sụt giảm hiệu suất và phá hủy xói mòn bề mặt cánh bơm, tua-bin, chân vịt tàu thủy và van điều tiết đường ống.

2. Khoảng trống nghiên cứu (Research Gap)

Mặc dù đã có nhiều mô hình mô phỏng CFD dòng xâm thực, phần lớn các phương pháp truyền thống tồn tại ba hạn chế lớn:

  • Bỏ qua tương tác giữa cấu trúc xoáy và áp suất chuyển pha: Áp suất tại tâm các xoáy cục bộ giảm sâu hơn nhiều so với dòng chảy bao quanh, nhưng các mô hình thương mại thường chỉ sử dụng áp suất hơi bão hòa tĩnh của môi chất ở trạng thái nghỉ làm ngưỡng chuyển pha.
  • Chưa tích hợp đầy đủ hiệu ứng nhiệt động lực học: Khi chất lỏng hóa hơi, nhiệt ẩn hóa hơi bị hấp thụ làm giảm nhiệt độ cục bộ tại giao diện pha (temperature depression), từ đó làm giảm áp suất bão hòa và kìm hãm sự phát triển của xâm thực. Hiệu ứng này đặc biệt quan trọng trong nước nóng hoặc môi chất lỏng đông lạnh (cryogenic fluids) trong động cơ tên lửa đẩy.
  • Đơn giản hóa vai trò của khí không hòa tan (Non-condensable Gas): Trong nước thực tế luôn tồn tại một lượng vi bọt khí không hòa tan, làm dịch chuyển áp suất ngưỡng chuyển pha sang $p = p_v + p_a$, khiến xâm thực khởi phát sớm hơn.

3. Tính thời điểm và tác động

Đề tài mã số QG.32 do TS. Lê Đình Anh (Trường Đại học Công nghệ – ĐHQGHN) chủ nhiệm là một nghiên cứu kịp thời, giải quyết trọn vẹn bài toán ghép nối đa vật lý (xoáy - nhiệt - khí không hòa tan) nhằm làm chủ công nghệ mô phỏng thủy khí động lực học tiên tiến tại Việt Nam.

Phương pháp tiếp cận và mô hình toán học (Methodology & Approach)

Nghiên cứu áp dụng quy trình mô phỏng số tiên tiến dựa trên phương pháp đa pha đồng nhất (Homogeneous Multiphase Model), tối ưu hóa hiệu năng tính toán nhưng vẫn bảo toàn đầy đủ các hiện tượng vật lý cốt lõi.

1. Hệ phương trình động lực học chất lưu đa pha

Dòng ba pha (nước lỏng – hơi bão hòa – khí không hòa tan) được đơn giản hóa thành mô hình hai pha nén được giả định ở trạng thái cân bằng động học và nhiệt động học:

  • Phương trình khối lượng riêng hỗn hợp $\rho_m$: $$\rho_m = (1 - \alpha)\rho_l + \alpha \rho_g$$ Trong đó mật độ pha lỏng $\rho_l$ tuân theo phương trình trạng thái Tamman, mật độ pha khí $\rho_g$ (gồm hơi bão hòa và khí không hòa tan) tuân theo phương trình khí lý tưởng: $$p = p_v + p_a = \rho_g R_g T$$
  • Phương trình năng lượng có xét nhiệt ẩn chuyển pha: Tích hợp số hạng nguồn nhiệt $L \cdot \dot{m}$ với $L = h_v - h_l$ là nhiệt ẩn hóa hơi, mô tả chính xác sự sụt giảm nhiệt độ tại vùng tạo bọt.
  • Mô hình chuyển pha Herz-Knudsen-Langmuir cải tiến: Xác định tốc độ hóa hơi $\dot{m}^+$ và ngưng tụ $\dot{m}^-$ dựa trên độ chênh lệch áp suất cục bộ và ngưỡng áp suất hiệu chỉnh.

2. Mô hình ngưỡng áp suất chuyển pha hiệu chỉnh (SVM)

Điểm cải tiến mang tính đột phá của nghiên cứu là xây dựng mô hình ngưỡng áp suất chuyển pha $p_v^$ tính đến tác động cục bộ của độ xoáy và biến dạng cắt: $$p_v^ = p_v(T) + p_{ff}$$ Với số hạng dao động áp suất dòng xoáy $p_{ff}$ được xác định theo điều kiện tương quan giữa độ lớn độ xoáy $\Omega = |\nabla \times \mathbf{u}|$ và độ lớn biến dạng cắt $S = \sqrt{2 \mathbf{S}:\mathbf{S}}$:

$$p_{ff} = \begin{cases} \mu S + C_k \rho k, & \text{khi } |\Omega| \le |S| \ \mu S + \mu \Omega + C_k \rho k, & \text{khi } |\Omega| > |S| \end{cases}$$

3. Thuật toán giải số và tính toán song song

  • Phương pháp sai phân hữu hạn (FDM): Sử dụng sơ đồ sai phân đón gió cấp hai Harten-Yee TVD (Total Variation Diminishing) giúp bắt chính xác các mặt gián đoạn, sóng áp suất và giao diện bọt khí mà không gây dao động số học phi vật lý.
  • Kỹ thuật tách toán tử thời gian (Strang Splitting): Tách riêng phần đối lưu, khuếch tán và số hạng nguồn chuyển pha cứng (stiff source terms) để đảm bảo tính ổn định thời gian.
  • Mã nguồn tính toán: Viết hoàn toàn bằng ngôn ngữ FORTRAN, tối ưu hóa tính toán song song với thư viện OpenMP trên môi trường máy trạm đa lõi Linux.

Các phát hiện và kết quả chính (Key Findings & Validation)

Nghiên cứu đã tiến hành kiểm nghiệm toàn diện mô hình SVM trên 3 đối tượng chuẩn quốc tế với các chế độ dòng chảy phức tạp:

1. Dòng xâm thực không dừng qua biên dạng NACA0015 ($T_0 = 293\text{K}$)

  • Phân bố hệ số áp suất: Khi góc tấn $\text{AoA} = 8^\circ$ và số xâm thực $\sigma = 1.4$, hiện tượng xâm thực mây (cloud cavitation) không dừng xuất hiện. Mô hình SVM cho kết quả phân bố hệ số áp suất trung bình $-C_p$ bám sát đường thực nghiệm của Kawanami et al., khắc phục hoàn toàn hiện tượng ước tính thiếu áp suất cục bộ của các mô hình HEM cũ.
  • Cấu trúc bọt khí: Tái hiện chính xác chu kỳ hình thành lớp xâm thực dừng ở mép trước, sự xuất hiện của dòng chảy ngược (re-entrant jet), sự đứt gãy tạo thành mây xoáy xâm thực và quá trình sụp đổ về phía đuôi cánh.

2. Dòng xâm thực nhiệt độ cao trong vòi phun hội tụ – phân kỳ ($T_0 > 373\text{K}$)

  • Hiệu ứng sôi bùng (Flashing Flow): Kiểm chứng với bộ số liệu thực nghiệm kinh điển của Abuaf et al. (trường hợp BNL133 và BNL296). Áp suất giảm đột ngột tại họng vòi phun dẫn đến sự hóa hơi diện rộng trong vùng khuếch tán, đẩy vận tốc dòng chảy lên vận tốc âm thanh hai pha ($\text{Mach} \approx 0.6$).
  • Sụt giảm nhiệt độ (Temperature Depression): Sự bay hơi mãnh liệt lấy đi nhiệt ẩn, làm nhiệt độ nước tại thành vòi phun lối ra giảm mạnh.
  • Quy luật thang thời gian vô thứ nguyên ($\tau / \tau_{th}$): Nghiên cứu phát hiện rằng sự suy giảm thể tích khối bọt khí không chỉ phụ thuộc vào nhiệt độ nước $T_0$ mà chịu sự chi phối quyết định bởi tỷ lệ giữa thời gian đặc trưng dòng chảy $\tau = D_{th}/U_{th}$ và thời gian nhiệt đặc trưng $\tau_{th}$. Khi $\tau / \tau_{th}$ đủ lớn, hiệu ứng kìm hãm xâm thực do nhiệt mới biểu hiện rõ rệt.

3. Tác động tương hỗ giữa nhiệt độ và khí không hòa tan

  • Dịch chuyển vị trí khởi phát xâm thực: Khi tính đến phần trăm thể tích khí không hòa tan $\alpha_0 = 0.1%$ ở $T_0 = 343\text{K}$, áp suất riêng phần của khí $p_a^*$ tại mép trước cánh $(x/c < 0.2)$ tăng cao đột biến.
  • Tăng tốc chuyển pha: Sự hiện diện của khí không hòa tan làm hạ thấp rào cản năng lượng chuyển pha, khiến quá trình tạo bọt diễn ra sớm hơn và vị trí bám của mây xâm thực dịch chuyển sát hơn về phía mép đón (leading edge) so với mô hình không xét khí (HEM).

Bảng tổng hợp các công bố khoa học và sản phẩm đào tạo

Đề tài hoàn thành vượt mức các chỉ tiêu nghiệm thu được phê duyệt:

STT Tên công trình / Sản phẩm Tạp chí / Nhà xuất bản / Cơ sở Phân loại Tình trạng
1 Improvement of Mass Transfer Rate Modeling for Prediction of Cavitating Flow Journal of Applied Fluid Mechanics (JAFM), Vol. 15(2), pp. 332-341 (2022) ISI/Scopus Q2 Đã xuất bản
2 Assessment of Flow Fluctuation Pressure Models for Simulating the Cavitating Flow Technical Physics Letters, Vol. 48(4) (2022) ISI/Scopus Q3 Đã xuất bản
3 Numerical Simulation of Hot Water Flashing Flow in a Converging - Diverging Nozzle Lecture Notes in Mechanical Engineering (Springer), pp. 62-71 (2022) Scopus Q4 Đã xuất bản
4 Bộ mã nguồn phần mềm mô phỏng xâm thực 2D (FORTRAN + OpenMP) ĐHQG Hà Nội lưu trữ Sản phẩm KH&CN Đã nghiệm thu
5 Luận văn Thạc sĩ: Nghiên cứu đặc tính thủy lực và xâm thực dòng 2D tốc độ cao... Học viên Nguyễn Văn Thanh (Đại học Bách Khoa Hà Nội) Đào tạo Thạc sĩ Đã bảo vệ thành công

Đóng góp khoa học và giá trị ứng dụng (Scientific & Practical Contributions)

1. Đóng góp lý thuyết học thuật

  • Hoàn thiện khung lý thuyết mô phỏng dòng nhiều pha bằng cách thiết lập mối quan hệ định lượng giữa ten-xơ biến dạng cắt, ten-xơ xoáy cục bộ và áp suất bão hòa chuyển pha.
  • Làm sáng tỏ bản chất nhiệt động lực học trong các dòng chảy tốc độ cao có kèm chuyển pha nhanh (cavitation/flashing), cung cấp cơ sở dữ liệu số đối chứng chất lượng cao cho cộng đồng cơ học thủy khí quốc tế.

2. Giá trị thực tiễn và kỹ thuật

  • Công nghiệp đóng tàu & Hàng hải: Hỗ trợ tính toán phân bố tải trọng và dự báo vùng xâm thực trên cánh chân vịt, giảm nguy cơ thủng vỏ chân vịt và giảm phát xạ tiếng ồn dưới nước.
  • Năng lượng & Thủy điện: Giúp tối ưu hóa hình học cánh tua-bin Francis, Kaplan và máy bơm ly tâm làm việc ở điều kiện cột áp cao hoặc môi chất nước nóng trong nhà máy nhiệt điện/hạt nhân.
  • Công nghệ Hàng không Vũ trụ: Phân tích chính xác hiện tượng xâm thực nhiệt trong bơm tăng áp nhiên liệu lỏng (oxy lỏng, hydro lỏng, ni-tơ lỏng) của tên lửa đẩy.

Đối tượng hưởng lợi (Target Audience)

  • Các nhà nghiên cứu và học giả chuyên ngành Cơ học Thủy khí (CFD / Multiphase Flows): Tiếp cận mô hình toán học SVM mới, phương pháp giải số FDM TVD bậc cao và các phát hiện vật lý về hiệu ứng nhiệt.
  • Kỹ sư R&D trong ngành Chế tạo Máy Thủy lực và Đóng tàu: Ứng dụng các quy tắc thiết kế biên dạng cánh, vòi phun để giảm thiểu xói mòn và nâng cao hiệu suất làm việc của thiết bị.
  • Các cơ quan quản lý và đơn vị đào tạo sau đại học: Tài liệu tham khảo chuẩn mực cho các chương trình đào tạo Thạc sĩ, Tiến sĩ chuyên ngành Kỹ thuật Cơ khí Động lực, Kỹ thuật Hàng không Vũ trụ và Kỹ thuật Nhiệt.

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

1. Điểm cốt lõi giúp mô hình SVM vượt trội hơn các mô hình xâm thực truyền thống là gì?

Mô hình SVM không cố định ngưỡng áp suất chuyển pha bằng giá trị áp suất bão hòa tĩnh $p_v(T)$, mà bổ sung số hạng áp suất động lực học $p_{ff}$ phụ thuộc trực tiếp vào độ lớn biến dạng cắt $S$ và độ xoáy $\Omega$. Điều này phản ánh chính xác hiện tượng áp suất sụt giảm cục bộ tại tâm các cấu trúc xoáy rối, giúp bắt kịp các hiện tượng xâm thực không dừng dạng mây xoáy mà mô hình cũ bỏ sót.

2. Tại sao hiệu ứng nhiệt động lực học lại có thể kìm hãm sự phát triển của bọt khí xâm thực?

Quá trình chuyển pha lỏng $\to$ hơi là quá trình thu nhiệt. Nhiệt ẩn hóa hơi được lấy trực tiếp từ lớp chất lỏng bao quanh bọt khí, khiến nhiệt độ tại bề mặt phân pha sụt giảm. Khi nhiệt độ giảm, áp suất hơi bão hòa tại chỗ giảm theo, làm suy yếu động lực chênh lệch áp suất gây hóa hơi, từ đó làm kích thước khối bọt khí bị co lại.

3. Khí không hòa tan trong nước ảnh hưởng như thế nào đến độ chính xác của mô phỏng?

Nước thực tế luôn chứa bọt khí không hòa tan. Khí này làm áp suất trong bọt khí trở thành $p = p_v + p_a$, khiến quá trình xâm thực khởi phát ở mức áp suất cao hơn và diễn ra sớm hơn (đặc biệt tại mép đón của cánh thủy lực). Bỏ qua khí không tan sẽ dẫn đến việc dự đoán sai lệch vị trí hình thành bọt khí đầu tiên.

4. Kết quả nghiên cứu 2D này có thể mở rộng cho bài toán thực tế 3D không?

Hoàn toàn có thể. Cấu trúc phương trình toán học và mô hình SVM được xây dựng dưới dạng giải tích ten-xơ tổng quát. Nhóm nghiên cứu đang phát triển mã nguồn mở rộng không gian 3 chiều (3D) để mô phỏng dòng xoáy đầu cánh (tip-vortex cavitation) trên chân vịt thực tế.

5. Mã nguồn FORTRAN của đề tài có ưu thế gì so với việc dùng phần mềm thương mại như ANSYS Fluent?

Mã nguồn phát triển riêng cho phép can thiệp trực tiếp vào cấu trúc hạt nhân của bộ giải số (explicit TVD scheme + Strang splitting), tùy biến mô hình chuyển pha tức thời mà không bị giới hạn bởi các mô hình đóng hộp (black-box) của phần mềm thương mại, đồng thời tối ưu hóa tính toán song song OpenMP trên hệ thống Linux hiệu năng cao.

Kết luận (Conclusion)

Đề tài nghiên cứu khoa học cấp ĐHQGHN QG.32 do TS. Lê Đình Anh chủ trì đã hoàn thành xuất sắc các mục tiêu đề ra, đem lại những đóng góp khoa học đột phá trong lĩnh vực mô phỏng số dòng chảy xâm thực đa pha. Việc tích hợp thành công mô hình ngưỡng áp suất chuyển pha SVM cùng hiệu ứng nhiệt động lực học và khí không hòa tan đã nâng cao đáng kể độ tin cậy của các dự báo CFD thủy lực.

Hướng phát triển tiếp theo của nghiên cứu sẽ tập trung vào việc mở rộng mô phỏng không gian 3 chiều (3D) và tích hợp kỹ thuật thích nghi lưới động (Adaptive Mesh Refinement) để giải quyết các bài toán công nghiệp phức tạp với độ chính xác và quy mô lớn hơn.