Báo cáo khoa học:
Mô hình số dòng chảy rối hai pha trong kênh dẫn
Tạp chí KHKT Nông nghiệp, Tập 1, số 2/2003

147

Mô hình số dòng chảy rối hai pha trong kênh dẫn
Numerical simulation on two-phase turbulent flow in channel
Hoàng Đức Liên
1
, Nguyễn Thị Thanh Loan
1

summary
This article presents a numerical model of two-phase non-isothermal turbulent
flow in channel basing on the model of Deitr & Philipov for the flow of steam, taking
in consideration the force of mutual influences between gas and particles.
Keywords: Numerical model, two-phase, non-isothermal turbulent flow. 1

11
u


11
v


0
PT động
lợng của pha
11
u


2
11
up

+

111
uv


z
R
2



111
ce +


(
)
2/
2
1111
chu +


(
)
2/
2
111
chv +


(
)
QvRuR
yz
+
+

222



z
R
2


PT CĐ ngang
của pha rắn

22
v


222
vu


2
21
v


y
R
2


PT năng
lợng của pha
22
e

,,
là
các thành phần lực tơng hỗ giữa các pha
và nhiệt lợng của một đơn vị thể tích của
pha rắn;
2

- mật độ thể tích tơng đối của
pha rắn. Khối lợng riêng của pha rắn có
thể xác định theo công thức:
k



22
=
.
Tiết diện ngang tính toán có dạng chữ
nhật nên sử dụng lới chữ nhật để tính
toán. Dựa theo theo đặc điểm của hỗn hợp
chất rắn và chất lỏng [3], [4] ta xác định
điều kiện biên của bài toán nh sau: tại
tiết diện ban đầu của dòng chảy có thể coi
sự phân bố của Entanpi h
0
, Entropi S là
không đổi và phơng của dòng chảy là
nh nhau. Tơng tự tại tiết diện ra khỏi
kênh kênh của dòng chảy, phân bố của áp
suất tĩnh sau tiết diện của lới đều. Trên

ở đây
2/1,2/1
++
nm

là diện tích tiết diện
biên của phần tử khối, các chỉ số bên dới
đợc sử dụng cho các giá trị của các thông
số tại thời điểm t và các chỉ số bên trên là
tại thời điểm (t+t).
Để hoàn chỉnh sơ đồ tính cần định rõ
phơng pháp tính các thông số tại biên của
các phần tử khối. Với sự góp mặt của các
quá trình trao đổi khối lợng, động lợng
và năng lợng giữa các pha, bài toán
không thể giải đợc bằng giải tích. Để
khắc phục vấn đề này việc tính toán các
giá trị các thông số của pha khí tại các
phần tử khối tại thời điểm sau đợc chia
làm hai bớc. Trong bớc đầu đầu ta chỉ
tính toán giá trị của các thông số tại biên
của phần tử khối mà không kể đến ảnh
hởng của sự tơng tác giữa các pha. Tác
động của chúng chỉ đợc xét đến trong hệ
phơng trình cơ bản dới dạng sai phân
(2) trong bớc hai.
Các thông số của pha rắn tại biên của
phần tử khối cũng đợc xác định qua lời
giải của bài toán dòng chảy một chiều
không ổn định của các phần tử khí với sự

ff
t
zcybzcybzcybzcyb
t
aa

(2)Hoàng Đức Liên, Nguyễn Thị Thanh Loan

149

diện ban đầu và giả thiết không có sự
tơng tác giữa các pha. Trong mô hình của
lời giải có sử dụng phép lặp tuyến tính khi
tính toán các thông số của khí và hạt pha
rắn tại biên của phần tử khối.
4. kết quả tính toán
Kết quả tính toán cho một số trờng
hợp đợc cho trên hình vẽ có các điều kiện
nh sau:
1) Vận tốc của pha khí (chỉ có theo
phơng của dòng chảy): 10 m/s;
2) Vận tốc của các hạt chất lỏng: tơng
tự nh pha khí;
3) áp suất d tại tiết diện vào: 0 Pa;
4) áp suất chân không tại tiết diện ra:
20 Pa.


khí (Hình 1). Còn sự phân bố áp suất
(Hình 3), Mật độ của pha rắn (Hình 4) và
Lực tơng tác giữa các pha rắn và khí
(Hình 5) chỉ thấy có sự biến đổi mạnh ở
đoạn đầu của dòng chảy. Điều đó chứng tỏ
kết quả tính toán số có sự phù hợp khá tốt
so với thực tế của dòng chảy.

5. Kết luận
Mô hình tính toán của bài toán trên là
mô hình số cấp thấp do cha kể đến ứng
suất rối. Kết quả tính toán cho thấy sự
phân bố không ổn định của các thông số
vật lý của dòng chảy từ tiết diện này sang
tiết diện khác theo chiều chuyển động.
Công việc nghiên cứu mô hình tiếp theo sẽ
bao gồm các nội dung:
Lựa chọn phơng pháp thích hợp cho
việc tính toán các thông số vật lý tại lớp
biên;


Numerical model of two-phase flow,
Scientific conference EMF01, Sofia, Vol.
III, pp. 42-47, Bulgarian.