Giả thiết cơ bản của mô hình đồng nhất trong dòng hai pha là gì?

Mô hình đồng nhất dựa trên hai giả thiết tiên quyết: các pha chuyển động cùng một vận tốc tức thời (không có hiện tượng trượt pha, tương ứng với tỷ số trượt bằng một) và các pha đạt trạng thái cân bằng nhiệt động học hoàn toàn tại mọi điểm dọc theo dòng chảy.

Mô hình đồng nhất khác mô hình thông lượng trôi như thế nào?

Trong khi mô hình đồng nhất giả định hai pha chuyển động với cùng vận tốc tuyệt đối, mô hình thông lượng trôi của Zuber và Findlay cho phép tồn tại vận tốc trôi tương đối giữa các pha do lực đẩy nổi và sự phân bố nồng độ không đồng đều trên tiết diện ống.

Khi nào không nên sử dụng mô hình đồng nhất?

Không nên áp dụng mô hình đồng nhất cho các dòng chảy phân tầng, dòng lượn sóng hoặc dòng chảy thẳng đứng có sự chênh lệch lớn về tỷ trọng giữa các pha, bởi vì hiện tượng trượt pha lớn sẽ khiến mô hình đánh giá quá cao tỷ lệ thể tích pha khí và đánh giá thấp nghiêm trọng lượng chất lỏng đọng lại trong ống.

Mô hình đồng nhất (homogeneous model) là gì?

Mô hình đồng nhất là mô hình thủy động lực học xem xét hỗn hợp dòng chảy đa pha như một chất lưu giả đơn pha đồng nhất, trong đó các pha cấu thành được giả định chuyển động cùng vận tốc tức thời và đạt trạng thái cân bằng nhiệt động học. Trong kỹ thuật nhiệt, công nghệ hóa chất, khai thác dầu khí và thủy lực nhiệt hạt nhân, mô hình này là công cụ nền tảng để dự đoán lưu lượng, độ sụt áp và hệ số truyền nhiệt. Bài viết này trình bày chi tiết bản chất vật lý, hệ phương trình bảo toàn vi phân, các phương pháp xác định thông số hiệu dụng, phân tích đối chiếu với các mô hình dòng chảy nâng cao, cùng các giới hạn ứng dụng thực tế.

Bản chất vật lý và các giả thiết cốt lõi

Trong tự nhiên và kỹ thuật, dòng chảy hai pha hoặc đa pha thường có cấu trúc giao diện phân tách pha vô cùng phức tạp, liên tục biến dạng và tương tác phi tuyến. Để đơn giản hóa bài toán mô phỏng mà vẫn giữ được các đặc trưng động lực học tổng thể, mô hình đồng nhất (homogeneous flow model) thực hiện phép lấy trung bình không gian và thời gian trên toàn bộ tiết diện dòng chảy. Thay vì giải các phương trình bảo toàn riêng biệt cho từng pha, mô hình này xem toàn bộ hỗn hợp như một môi trường liên tục đơn nhất với các tính chất vật lý tương đương.

Brennen (2005) tổng kết rằng mô hình dòng đồng nhất được xây dựng dựa trên hai giả thiết vật lý tiên quyết:

Giả thiết cân bằng động học (không có hiện tượng trượt pha): Các pha cấu thành, ví dụ như pha khí và pha lỏng, chuyển động với cùng một vectơ vận tốc tại mọi điểm trong trường dòng chảy. Nói cách khác, vận tốc tức thời của pha khí bằng vận tốc tức thời của pha lỏng, tương ứng với tỷ số trượt vận tốc giữa hai pha bằng 1.
Giả thiết cân bằng nhiệt động học: Nhiệt độ và áp suất giữa các pha đạt trạng thái cân bằng cục bộ tức thời tại mọi vị trí dọc theo chiều dài ống dẫn. Không có sự chênh lệch nhiệt độ hay áp suất qua bề mặt phân chia pha, nghĩa là quá trình truyền nhiệt và truyền khối lượng giữa các pha diễn ra với tốc độ đủ nhanh so với thời gian lưu chuyển của dòng chất lưu.

Khi cả hai giả thiết trên được thỏa mãn đồng thời, mô hình thường được gọi là mô hình cân bằng đồng nhất (homogeneous equilibrium model). Đây là công cụ xấp xỉ có độ phức tạp tính toán thấp nhất trong thủy động lực học dòng đa pha.

Hệ phương trình bảo toàn và mô hình hóa toán học

Dựa trên các giả thiết đồng nhất, hệ phương trình Navier-Stokes và phương trình năng lượng của dòng hai pha được rút gọn về dạng tương tự như dòng đơn pha, nhưng áp dụng cho các biến số trung bình của hỗn hợp:

Phương trình bảo toàn khối lượng liên tục cho hỗn hợp hai pha được biểu diễn dưới dạng vi phân vectơ:

\frac{\partial \rho_m}{\partial t} + \nabla \cdot (\rho_m \mathbf{u}_m) = 0

Trong phương trình này, ký hiệu $\rho_m$ biểu thị khối lượng riêng trung bình của hỗn hợp, ký hiệu $t$ là thời gian, và ký hiệu $\mathbf{u}_m$ là vectơ vận tốc khối lượng của hỗn hợp.

Phương trình bảo toàn động lượng của hỗn hợp được mô tả thông qua quan hệ cân bằng lực:

\frac{\partial (\rho_m \mathbf{u}_m)}{\partial t} + \nabla \cdot (\rho_m \mathbf{u}_m \mathbf{u}_m) = -\nabla p + \nabla \cdot \boldsymbol{\tau}_m + \rho_m \mathbf{g}

Trong đó, ký hiệu $p$ biểu thị áp suất thủy tĩnh cục bộ, ký hiệu $\boldsymbol{\tau}_m$ là tensor ứng suất nhớt hiệu dụng của hỗn hợp, và ký hiệu $\mathbf{g}$ là gia tốc trọng trường. Do giả thiết các pha chuyển động cùng vận tốc, số hạng ứng suất trượt do tương tác vận tốc giữa các pha hoàn toàn biến mất, giúp đơn giản hóa cấu trúc tensor ma sát nhớt.

Phương trình bảo toàn năng lượng cho dòng chảy một chiều dọc theo trục ống dẫn được biểu diễn qua entanpi hỗn hợp và mật độ dòng nhiệt truyền qua vách:

\frac{\partial (\rho_m h_m)}{\partial t} + \frac{\partial (\rho_m u_m h_m)}{\partial z} = \frac{\partial p}{\partial t} + u_m \frac{\partial p}{\partial z} - \frac{\partial q_m}{\partial z} + \frac{4 q_w}{D}

Ở đây, ký hiệu $h_m$ là entanpi riêng của hỗn hợp, ký hiệu $q_m$ là mật độ dòng nhiệt dẫn dọc trục, đại lượng $q_w$ là mật độ dòng nhiệt truyền từ vách thành ống vào dòng chất lưu, và ký hiệu $D$ là đường kính trong của ống dẫn.

Xác định các thông số vật lý hiệu dụng của hỗn hợp

Để khép kín hệ phương trình vi phân trên, người kỹ sư cần xác định các tính chất vật lý hiệu dụng của chất lưu giả đơn pha dựa trên tỷ lệ thành phần của từng pha đơn lẻ. Các thông số này phụ thuộc chặt chẽ vào độ khô chất lượng hoặc tỷ lệ thể tích pha:

Độ khô chất lượng: Đại lượng vô thứ nguyên ký hiệu là $x$ , biểu thị tỷ số giữa lưu lượng khối lượng của pha khí trên tổng lưu lượng khối lượng của hỗn hợp.
Tỷ lệ thể tích pha hơi: Ký hiệu là $\alpha$ , biểu thị tỷ lệ phần trăm diện tích hoặc thể tích mà pha khí chiếm dụng trong ống dẫn. Trong điều kiện không có sự trượt vận tốc, mối liên hệ giữa tỷ lệ thể tích và độ khô chất lượng được tính toán trực tiếp thông qua khối lượng riêng của pha khí và pha lỏng:

\alpha = \frac{1}{1 + \left( \frac{1-x}{x} \right) \frac{\rho_g}{\rho_l}}

Trong công thức trên, ký hiệu $\rho_g$ là khối lượng riêng của pha khí, và ký hiệu $\rho_l$ là khối lượng riêng của pha lỏng. Khối lượng riêng trung bình của hỗn hợp hai pha được tính toán theo một trong hai cách tương đương:

\rho_m = \alpha \rho_g + (1 - \alpha) \rho_l

\frac{1}{\rho_m} = \frac{x}{\rho_g} + \frac{1-x}{\rho_l}

Việc xác định độ nhớt động lực học hiệu dụng của hỗn hợp hai pha phức tạp hơn đáng kể do hiện tượng tương tác giữa các bọt khí hoặc giọt lỏng với trường vận tốc cục bộ. Awad và Muzychka (2008) đã hệ thống hóa và đánh giá các mô hình độ nhớt hiệu dụng phổ biến nhất:

Mô hình McAdams: Giả định nghịch đảo độ nhớt hỗn hợp bằng tổng có trọng số nghịch đảo độ nhớt từng pha:
$\frac{1}{\mu_m} = \frac{x}{\mu_g} + \frac{1-x}{\mu_l}$
Trong đó ký hiệu $\mu_g$ là độ nhớt động lực học pha khí và $\mu_l$ là độ nhớt pha lỏng. Mô hình này thường được ưu tiên khi tính toán sụt áp ma sát trong các dòng chảy tầng.
Mô hình Cicchitti: Sử dụng phép tính trung bình tuyến tính trực tiếp theo độ khô chất lượng:
$\mu_m = x \mu_g + (1-x) \mu_l$
Mô hình này đơn giản và thường cho kết quả tiệm cận tốt trong các dòng sôi cưỡng bức áp suất cao.
Mô hình Dukler: Tính toán độ nhớt hiệu dụng dựa trên sự phân bố thể tích thực tế của từng pha trong dòng chảy, phù hợp cho các chế độ dòng chảy phân tán mịn.

Bảng so sánh mô hình đồng nhất với các mô hình dòng hai pha

Để đánh giá vị trí và tính phù hợp của mô hình đồng nhất, bảng dưới đây đối chiếu đặc điểm tính toán và phạm vi ứng dụng của ba mô hình dòng hai pha kinh điển:

Tiêu chí so sánh	Mô hình đồng nhất	Mô hình thông lượng trôi	Mô hình hai chất lưu
Giả thiết vận tốc pha	Hai pha cùng vận tốc tuyệt đối	Có vận tốc trôi tương đối giữa các pha	Mỗi pha có trường vận tốc độc lập
Cân bằng nhiệt động	Cân bằng tức thời về nhiệt và áp	Có thể xét mất cân bằng nhiệt cục bộ	Mất cân bằng nhiệt và cơ học hoàn toàn
Số lượng phương trình	Ba phương trình bảo toàn hỗn hợp	Bốn hoặc năm phương trình bảo toàn	Sáu hoặc bảy phương trình bảo toàn
Độ phức tạp tính toán	Rất thấp, tính toán tức thời	Trung bình, phù hợp mô phỏng hệ thống	Rất cao, đòi hỏi tài nguyên tính toán lớn
Độ chính xác dòng bọt mịn	Chính xác cao, sai số nhỏ	Chính xác cao	Rất cao nhưng tốn kém
Độ chính xác dòng phân tầng	Kém, đánh giá thấp sụt áp	Trung bình	Rất cao, bắt trọn mặt phân giới

Tính toán độ sụt áp do ma sát trong đường ống

Một trong những ứng dụng phổ biến nhất của mô hình đồng nhất là tính toán tổn thất áp suất do ma sát dọc theo chiều dài đường ống dẫn chất lưu hai pha. Gradient áp suất ma sát một chiều được biểu diễn thông qua hệ số ma sát Darcy-Weisbach của hỗn hợp:

-\left( \frac{dp}{dz} \right)_{F} = \frac{f_m G^2}{2 D \rho_m}

Trong biểu thức này, ký hiệu $G$ là lưu lượng khối lượng trên một đơn vị diện tích mặt cắt ngang, ký hiệu $D$ là đường kính trong của ống dẫn, và ký hiệu $f_m$ là hệ số ma sát Darcy-Weisbach của hỗn hợp. Hệ số $f_m$ được tính toán dựa trên số Reynolds của hỗn hợp đồng nhất:

Re_m = \frac{G D}{\mu_m}

Để liên kết mô hình đồng nhất với các quan sát thực nghiệm trong dòng phân tách, Chisholm (1967) đã xây dựng cơ sở lý thuyết chặt chẽ cho hệ số nhân độ sụt áp hai pha. Hệ số nhân này cho phép quy đổi độ sụt áp ma sát tính theo chất lỏng đơn pha sang độ sụt áp thực tế của dòng hai pha thông qua tỷ số giữa khối lượng riêng và độ nhớt của hai pha, thiết lập cầu nối toán học quan trọng giữa phương pháp tiếp cận đồng nhất và phương pháp bán thực nghiệm.

Ứng dụng kỹ thuật và ý nghĩa trong bối cảnh Việt Nam

Tại Việt Nam, mô hình đồng nhất được ứng dụng rộng rãi trong các tính toán kỹ thuật ban đầu, thiết kế sơ bộ và thẩm tra an toàn hệ thống công nghiệp:

Hệ thống thu gom và vận chuyển dầu khí đa pha ngoài khơi: Trong các tuyến đường ống ngầm đáy biển vận chuyển hỗn hợp dầu thô, khí đồng hành và nước vỉa từ các giàn khai thác thềm lục địa về bờ, các kỹ sư khai thác dầu khí sử dụng mô hình đồng nhất để nhanh chóng ước tính áp suất đầu giếng và kích thước đường kính ống dẫn trước khi triển khai các phần mềm mô phỏng động học chuyên sâu.
Thiết kế lò hơi công nghiệp và thiết bị trao đổi nhiệt: Trong các nhà máy nhiệt điện than và các cơ sở công nghiệp chế biến hóa chất, việc phân tích dòng sôi hai pha trong các giàn ống sinh hơi dựa trên mô hình đồng nhất giúp tính toán phân bố nhiệt độ và dự báo nguy cơ khủng hoảng trao đổi nhiệt.
Kỹ thuật lạnh và điều hòa không khí: Trong ngành công nghiệp điện lạnh thương mại và dân dụng tại Việt Nam, mô hình đồng nhất được ứng dụng để định cỡ van tiết lưu và đường ống dẫn môi chất lạnh hai pha gas-lỏng trong các hệ thống điều hòa không khí trung tâm.

Giới hạn vật lý và phạm vi áp dụng

Mặc dù mang lại sự tiện lợi vượt trội về mặt toán học, mô hình đồng nhất có những ranh giới vật lý rõ ràng mà người kỹ sư cần lưu ý khi ứng dụng:

Hiện tượng trượt vận tốc pha là nguồn gốc sai số lớn nhất của mô hình. Trong các dòng chảy thẳng đứng hướng lên hoặc các dòng chảy có sự chênh lệch khối lượng riêng quá lớn giữa pha khí và pha lỏng, pha khí nhẹ hơn sẽ di chuyển nhanh hơn đáng kể so với pha lỏng do tác động của lực đẩy nổi. Zuber và Findlay (1965) đã chứng minh bằng thực nghiệm rằng trong chế độ dòng bọt lớn hoặc dòng nút bọt thẳng đứng, tỷ số trượt vận tốc thực tế có thể lớn hơn rất nhiều so với giá trị 1. Do giả định hai pha chuyển động với cùng vận tốc, mô hình đồng nhất đánh giá quá cao tỷ lệ thể tích pha khí và đánh giá thấp tỷ lệ chất lỏng đọng lại trong ống, dẫn đến sai số khi tính toán tổng khối lượng chất lưu và độ sụt áp tĩnh dọc theo tuyến ống.

Mặt khác, trong các bài toán an toàn nhiệt thủy lực hạt nhân liên quan đến hiện tượng dòng chảy tới hạn (critical flow) khi xảy ra sự cố vỡ đường ống làm mát, Tian và cs. (2025) chỉ ra rằng mô hình cân bằng đồng nhất thường dự đoán cận dưới của lưu lượng xả tới hạn. Nguyên nhân là do trong các quá trình giảm áp cực nhanh, sự hóa hơi bị trễ do hiện tượng chậm sôi nhiệt động học, khiến chất lỏng duy trì trạng thái quá nhiệt thay vì đạt trạng thái cân bằng pha tức thời như mô hình đồng nhất giả định. Do đó, đối với các phân tích an toàn chi tiết hoặc thiết kế van an toàn, các mô hình hai chất lưu không cân bằng nhiệt động cần được sử dụng thay thế.