Dòng chảy đa pha (tiếng Anh: multiphase flow; hay chuyển động đa pha) là dòng chuyển động đồng thời của hai hoặc nhiều pha vật chất nhiệt động học khác nhau (các trạng thái rắn, lỏng, khí) hoặc các chất lưu không hòa tan vào nhau (như nhũ tương dầu - nước) cùng tồn tại bên trong một không gian hình học dẫn dòng (như đường ống công nghệ, buồng đốt động cơ, lò phản ứng hạt nhân hoặc vỉa chứa địa chất). Điểm đặc trưng cốt lõi phân biệt dòng đa pha với dòng đơn pha là sự hiện diện của mặt phân cách pha (phase interface) động có thể biến dạng, đứt gãy hoặc hợp nhất liên tục dưới tác động của lực căng bề mặt, lực cắt nhớt và lực quán tính. Do sự chênh lệch lớn về khối lượng riêng, độ nhớt và vận tốc chuyển động tương đối (vận tốc trượt - slip velocity) giữa các pha, dòng đa pha thể hiện những cấu trúc thủy động lực học phi tuyến vô cùng phức tạp. Bài viết này trình bày toàn diện về bản đồ chế độ dòng chảy, các mô hình toán học giải tích và mô phỏng số CFD (Computational Fluid Dynamics), cùng các ứng dụng then chốt trong công nghiệp dầu khí, năng lượng nhiệt điện hạt nhân và kỹ thuật hóa học.
Phân loại các hệ dòng đa pha
Dựa trên sự kết hợp các trạng thái pha vật chất:
- Dòng hai pha Khí - Lỏng (Gas - Liquid Flow): Phổ biến nhất trong công nghiệp; ví dụ sự sôi sinh hơi nước trong ống lò hơi, vận chuyển dầu thô kèm khí đồng hành trong đường ống biển.
- Dòng hai pha Lỏng - Rắn (Liquid - Solid / Slurry Flow): Vận chuyển bùn quặng khoáng sản, nạo vét phù sa lòng sông, dòng huyền phù bột giấy.
- Dòng hai pha Khí - Rắn (Gas - Solid / Pneumatic Transport): Vận chuyển hạt nông sản (hạt ngô, hạt nhựa) bằng khí nén, dòng bụi than phun vào buồng đốt lò nhiệt điện, tầng sôi công nghiệp (Fluidized Bed).
- Dòng ba pha hoặc bốn pha (Multiphase Complex): Dòng hỗn hợp Dầu - Nước - Khí - Cát trong khai thác mỏ dầu khí ngoài khơi.
Các chế độ dòng chảy hai pha Khí - Lỏng trong đường ống
Sự phân bố không gian hình học của hai pha khí và lỏng tạo nên các chế độ dòng chảy đặc trưng (Flow Regimes) được mô tả trên Bản đồ dòng chảy Taitel - Dukler:
| Chế độ dòng chảy | Mô tả cấu trúc hình thái học | Đặc tính áp lực và dòng chảy |
|---|---|---|
| Dòng bọt khí (Dispersed Bubble Flow) | Pha khí phân tán thành các bọt khí nhỏ rời rạc chuyển động lơ lửng trong pha lỏng liên tục chiếm ưu thế. | Dòng chảy êm thuận, độ trượt giữa hai pha nhỏ; xảy ra khi lưu lượng lỏng rất cao và lưu lượng khí thấp. |
| Dòng phân tầng (Stratified Flow) | Pha lỏng nặng hơn chảy ở đáy ống, pha khí nhẹ hơn chảy ở phần trên của ống, ngăn cách bởi mặt phân cách nhẵn. | Thường xuất hiện trong ống nằm ngang ở vận tốc dòng chảy thấp. Khi vận tốc khí tăng, mặt phân cách sẽ sinh sóng tạo Dòng lượn sóng (Wavy flow). |
| Dòng nút khí - lỏng (Slug / Plug Flow) | Sóng trên mặt chất lỏng phát triển chạm trần ống tạo thành các "nút chất lỏng" (Liquid slugs) xen kẽ các bọc khí lớn hình viên đạn (Taylor bubbles). | Chế độ nguy hiểm nhất trong công nghiệp: Gây ra các xung áp lực dao động va đập thủy lực mạnh (slug hammering), gây rung lắc đường ống và sốc phụ tải thiết bị tách. |
| Dòng màng bao quanh (Annular Flow) | Pha lỏng tạo thành một màng mỏng bao quanh toàn bộ chu vi thành ống, trong khi pha khí chuyển động với vận tốc cực cao ở lõi trung tâm mang theo các giọt lỏng phân tán. | Xảy ra khi lưu lượng khí rất lớn; hiệu quả truyền nhiệt bốc hơi qua màng rất cao. |
Các phương pháp mô hình hóa toán học và mô phỏng số CFD
Để tính toán và mô phỏng dòng chảy đa pha trên máy tính, ba phương pháp tiếp cận chính được ứng dụng:
1. Mô hình Đồng nhất (Homogeneous Equilibrium Model - HEM)
Xem hỗn hợp đa pha như một chất lưu đơn pha giả định có các đặc tính khối lượng riêng trung bình (\( ho_{ ext{mix}}\)) và độ nhớt trung bình (\(\mu_{ ext{mix}}\)), giả thiết hai pha chuyển động với cùng vận tốc (không có độ trượt). Đơn giản, giải nhanh nhưng độ chính xác thấp đối với các chế độ dòng phân tầng hoặc dòng nút.
2. Mô hình Hai chất lưu Euler - Euler (Two-Fluid Eulerian Model)
Thiết lập hệ phương trình bảo toàn khối lượng, động lượng và năng lượng riêng biệt cho từng pha, đồng thời liên kết hai pha thông qua các số hạng tương tác liên pha (lực cản ma sát liên pha Drag force, lực nâng Lift force, lực khối lượng ảo Virtual mass và quá trình truyền nhiệt/chuyển pha). Đây là mô hình mạnh mẽ nhất để mô phỏng dòng tầng sôi và dòng bọt khí nồng độ cao.
3. Phương pháp Thể tích Chất lưu (Volume of Fluid - VOF)
Phương pháp chuyên dụng theo dõi ranh giới mặt phân cách tự do giữa hai pha không hòa tan (như sóng biển va đập chân đế giàn khoan, bọt khí nổi lên trong cột chất lỏng). VOF giải một phương trình Navier-Stokes chung và tính toán lực căng bề mặt qua mô hình CSF (Continuum Surface Force) của Brackbill.
Ứng dụng công nghiệp trọng điểm
- Khai thác và Vận chuyển Dầu khí tại Việt Nam: Mạng lưới đường ống ngầm đáy biển dẫn dầu khí từ các mỏ Bạch Hổ, Rồng Đôi, Lan Tây (Bể Cửu Long và Nam Côn Sơn) vào bờ vận chuyển dòng đa pha phức tạp. Tập đoàn Dầu khí Việt Nam (PVN) sử dụng phần mềm mô phỏng dòng đa pha chuyên dụng OLGA để thiết kế hệ thống bẫy nút chất lỏng (Slug Catcher) và bơm chất ức chế hydrat (MEG, Methanol) chống đông tắc đường ống biển.
- An toàn Lò phản ứng Hạt nhân và Nhà máy Nhiệt điện: Dự đoán chính xác hiện tượng khủng hoảng truyền nhiệt do sôi (Critical Heat Flux - CHF) và hiện tượng màng hơi cách nhiệt (Dryout) để ngăn ngừa thảm họa nóng chảy thanh nhiên liệu hạt nhân.