Nghiên cứu quá trình lan truyền, tương tác giữa các song trong một số hỗn hợp chất lỏng hai pha
Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (1.29 MB, 51 trang )
1
MỤC LỤC
LỜI CẢM ƠN 2
Danh sách các kí hiệu và chữ viết tắt 3
MỞ ĐẦU 4
CHƯƠNG I TỔNG QUAN 8
CHƯƠNG II CƠ SỞ LÝ THUYẾT 12
2.1. HỆ PHƯƠNG TRÌNH THUỶ- NHIỆT ĐỘNG LỰC HỌC CHO MÔI
TRƯỜNG HAI PHA LỎNG HƠI. 12
2.1.1. Xây dựng mô hình 12
2.1.2. Hệ phương trình cơ sở. 13
2.2. PHƯƠNG PHÁP GIẢI SỐ VÀ CHƯƠNG TRÌNH TÍNH TOÁN 26
2.2.1. Phương pháp giải số 26
2.2.2. Chương trình tính toán 29
CHƯƠNG III MỘT SỐ KẾT QUẢ NGHIÊN CỨU QUÁ TRÌNH LAN
TRUYỀN VÀ TƯƠNG TÁC SÓNG XUNG KÍCH TRONG CÁC HỖN HỢP
CHẤT LỎNG CHỨA BỌT HƠI. 31
3.1. QUÁ TRÌNH LAN TRUYỀN, TƯƠNG TÁC GIỮA SÓNG NGẮN
VÀ SÓNG NGẮN TRONG MỘT SỐ HỖN HỢP LỎNG - HƠI 31
3.1.1. Hỗn hợp nước chứa bọt hơi 31
3.1.2. Hỗn hợp freon21 chứa bọt hơi 36
3.1.3. Hỗn hợp nito lỏng chứa bọt hơi 38
3.2. QUÁ TRÌNH LAN TRUYỀN, TƯƠNG TÁC GIỮA SÓNG DÀI VÀ
SÓNG NGẮN TRONG HỖN HỢP NƯỚC SÔI CHỨA BỌT HƠI 42
KẾT LUẬN 47
KIẾN NGHỊ VỀ NHỮNG NGHIÊN CỨU TIẾP THEO 48
TÀI LIỆU THAM KHẢO 49
1Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
2
LỜI CẢM ƠN
Để hoàn thành luận văn Thạc sĩ của mình em xin được gửi lời cảm ơn
chân thành tới Ban giám hiệu, phòng Đào tạo, Khoa Toán – Tin, các thầy cô
giáo trường Đại học Khoa học – Đại học Thái Nguyên và các thầy đang công
tác tại Viện Toán học Việt Nam đã nhiệt tình truyền đạt cho em những kiến
thức quý báu trong suốt quá trình học tập và hoàn thành luận văn Thạc sĩ.
Đặc biệt, em xin được gửi lời biết ơn chân thành nhất tới PGS.TS.
Nguyễn Văn Tuấn, người trực tiếp chỉ bảo, hướng dẫn và động viên em trong
suốt quá trình nghiên cứu và hoàn thành luận văn.
Xin được gửi lời cảm ơn đến BGH, các bạn đồng nghiệp của tôi tại
trường PTDT Nội trú cấp II-III Bắc Quang, Hà Giang đã quan tâm, tạo mọi
điều kiện thuận lợi nhất để tôi hoàn thành khóa học. Xin cảm ơn các anh chị
em học viên lớp cao học K4A đã đoàn kết, đùm bọc và giúp đỡ nhau trong
toàn khóa học.
Cuối cùng xin được gửi lời biết ơn sâu sắc đến những người thân trong
gia đình tôi, những người luôn động viên, khuyến khích và giúp đỡ tôi trong
suốt quá trình học tập. Thành quả đạt được chính là món quà mà tôi muốn
dành tặng gia đình thân yêu của mình.
Thái Nguyên, tháng 11 năm 2012
Học viên
Nguyễn Đức Chung
2Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
3
Danh sách các kí hiệu và chữ viết tắt
B - Hằng số khí
22
,,
vp
ccc
- Nhiệt dung riêng, nhiệt dung riêng khi áp suất và
vận tốc không đổi
j - Cường độ chuyển pha
l - Nhiệt hoá hơi của chất lỏng
n - Số lượng bọt
p,p
0
- áp suất hỗn hợp, áp suất ban đầu
p
i
- áp suất của pha i
p
e
- Cường độ sóng xung kích ban đầu
q - Dòng nhiệt
R,a - Bán kính của bọt
R
0
, a
0
- Bán kính ban đầu của bọt
T,T
0
- nhiệt độ hỗn hợp, nhiệt độ ban đầu
T
i
- nhiệt độ pha i
t - Thời gian
v, v
i
- Vận tốc của hỗn hợp, vận tốc của pha thứ i
w
1
- Vận tốc hướng kính của chất lỏng chuyển động
xung quanh bọt.
a
w
1
- Vận tốc màng bọt
WP,,,
- Đại lượng không thứ nguyên của mật độ, nhiệt độ,
áp suất và vận tốc màng bọt
i
- Phần thể tích của pha i trong hỗn hợp
0i
- Phần thể tích ban đầu của pha i
- Hệ số sức căng bề mặt
- Toạ độ Lagrange
21
,
- Hệ số hiệu chỉnh của hỗn hợp
1
- Hệ số nhớt động lực của chất lỏng
i
,
- mật độ của pha i, mật độ của hỗn hợp
0
i
,
0
- mật độ thực của pha i, mật độ ban đầu của hỗn hợp
- tỷ số giữa nhiệt dung riêng của hơi khi áp suất
không đổi và thể tích không đổi
4321
,,, bbbb
- cường độ của các xung áp suất
3Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
4
MỞ ĐẦU
Chất lỏng không đồng nhất hay hỗn hợp chất lỏng nhiều pha là một
môi trường rất phức tạp. Môi trường này thường gặp nhiều trong các quá
trình tự nhiên, trong công nghệ hoá học, vật lý… và cũng được ứng dụng rất
nhiều trong các nghành công nghiệp năng lượng như : khai thác , vận chuyển
và chế biến dầu khí…Tuy nhiên, chỉ từ 1950 trở lại đây, việc nghiên cứu về
dòng hai pha khí lỏng mới được tiến hành một cách có hệ thống cả về lý
thuyết và thực nghiệm. Trong các hỗn hợp này , thì quá trình trao đổi nhiệt -
chất là một trong những hiện tượng quan trọng không thể tách rời, nhất là
trong trường hợp tồn tại sóng xung kích lan truyền trong hỗn hợp.
So với môi trường là chất lỏng đồng nhất một pha thì môi trường hỗn
hợp hai pha khác xa về tính chất vật lý, nó thể hiện ở chỗ trong hỗn hợp do
có sự kết hợp các tính chất phi tuyến vật lý mạnh, sự tán sắc và quá trình hao
tán năng lượng nên biểu đồ mô tả các sóng có nhiều dạng. Chính vì vậy, khi
thay đổi các điều kiện thuỷ động lực sẽ dẫn đến sự thay đổi cấu trúc về sóng,
các tính chất vật lý nhiệt, và các quá trình tương tác giữa các pha. Tính chất
đặc trưng của hỗn hợp chất lỏng chứa bọt hơi trong các quá trình động lực
học là sự xuất hiện biến dạng cục bộ của hỗn hợp khi thay đổi thể tích môi
trường do sự thay đổi thể tích của bọt. Khả năng mức độ co nén của bọt phụ
thuộc mạnh vào sự trao đổi nhiệt và khối lượng giữa pha lỏng và pha khí. Sự
xuất hiện đồng thời những năng lượng do biến dạng này của sẽ dẫn tới sóng
có cấu trúc khác nhau. Ngoài ra, sự truyền sóng áp suất trong những môi
trường như vậy cũng dẫn đến khả năng hoá hơi và ngưng tụ của pha khí, từ
đó sẽ dẫn đến sự thay đổi chủ yếu cấu trúc của môi trường.
Do hỗn hợp chất lỏng chứa bọt hơi có tính chất đặc biệt như trên,
hơn nữa đây là hỗn hợp xuất hiện rất nhiều trong các lĩnh vực công nghiệp
4Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
5
năng lượng, công nghệ hoá học và các quá trình tự nhiên cho nên, sự hiểu
biết về các hiện tượng có thể xuất hiện khi sóng áp suất lan truyền qua chất
lỏng chứa bọt và nhất là khi xảy ra quá trình tương tác giữa các sóng là rất
cần thiết.
Căn cứ vào tình hình phát triển của các nghiên cứu về các quá trình
lan truyền của sóng xung kích trong hỗn hợp các chất lỏng- bọt trong và
ngoài nước mục đích chính của đề tài này bao gồm các vấn đề sau:
- Nghiên cứu lý thuyết về quá trình sóng trong hỗn hợp hai pha bao
gồm pha phân tán là các phần tử dạng bọt hơi cùng với pha chính là chất lỏng
Newton được chứa trong một ống nằm ngang. Giả sử tại hai đầu ống tồn tại
hai sóng xung kích lan truyền vào trong hỗn hợp, chúng tương tác và lan
truyền ra ngược nhau. Thành lập hệ phương trình vi phân đạo hàm riêng mô
tả các quá trình này .
- Nghiên cứu, phân tích và đánh giá các quy luật phụ thuộc chung của
sóng xung kích trong các hỗn hợp chất lỏng chứa bọt hơi trong các quá trình
lan truyền và quá trình tương tác giữa các sóng trên cơ sở của hệ phương trình
thủy - nhiệt động lực học. Sau khi đã biến đổi về dạng phù hợp với việc sử
dụng phương pháp số, sử dụng thuật toán và chương trình tính phù hợp với
mô hình khảo sát. Chương trình đã được kiểm chứng bằng cách so sánh kết
quả của chương trình với kết quả thực nghiệm đã được công bố của các tác
giả khác về sóng tới.
- Sử dụng chương trình tính nói trên để nghiên cứu và phân tích sự phụ
thuộc của cường độ sóng áp suất vào: các điều kiện đầu, điều kiện biên,
cường độ của sóng xung kích ban đầu, thể tích của pha hơi và các tính chất
vật lý nhiệt của hỗn hợp, trong các quá trình lan truyền và tương tác của các
sóng xung kích trong một số hỗn hợp chất lỏng chứa bọt hơi.
5Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
6
- Đưa ra nhận xét và kết luận những vấn đề nghiên cứu, đồng thời đề
xuất ý kiến góp phần giảm thiểu những hiện tượng bất thường xảy ra để đảm
bảo an toàn trong hệ thống truyền dẫn thủy lực.
Nội dung của luận văn được chia làm ba chương :
Chương I: Tổng quan
Trình bày tổng quan về sự phát triển và xu hướng phát triển trong lĩnh
vực truyền sóng xung kích trong hỗn hợp lỏng – hơi, các hiện tượng xảy ra
khi sóng lan truyền trong hỗn hợp.
Chương II: Cơ sở lý thuyết
Dựa trên cơ sở hệ phương trình thủy - nhiệt động lực học của hỗn hợp
lỏng hơi mô tả hiệ tượng vật lý, sử dụng thuật toán và chương trình tính phù
hợp để nghiên cứu, phân tích và đánh giá quá trình lan truyền và tương tác
sóng của các sóng xung kích trong hỗn hợp.
Chương III: Một số kết quả nghiên cứu quá trình lan truyền, tương tác
sóng xung kích trong các hỗn hợp chất lỏng chứa bọt hơi.
Trình bày một số kết quả nghiên cứu về quá trình la truyền và tương
tác giữa sóng ngắn và sóng ngắn, giữa sóng dài và sóng ngắn. Từ các kết quả
nhận được sẽ nghiên cứu, đánh giá sự ảnh hưởng của các yếu tố cường độ
ban đầu của sóng, nồng độ thể tích pha hơi lên quá trình lan truyền và tương
tác sóng trong hỗn hợp.
So sánh các quá trình tương tác giữa các sóng trong một số hỗn hợp
chất lỏng chứa bọt hơi nhằm chỉ ra được sự ảnh hưởng của các điều kiện
đầu, điều kiện biên và các tính chất vật lý nhiệt của hỗn hợp lên các quá
trình lan truyền và tương tác sóng trong hỗn hợp.
6Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
7
Để thự hiện các nội dung trên, tôi đã tiến hành thực hiện đề tài
“Nghiên cứu quá trình lan truyền, tương tác giữa các song trong một số hỗn
hợp chất lỏng hai pha ”.
Bản luận văn này được thực hiện tại trường Đại học Khoa học, Đại
học Thái Nguyên và hoàn thành từ ý tưởng nghiên cứu, sự hướng dẫn tận
tình và giúp đỡ về mặt khoa học của PGS. TS. Nguyễn Văn Tuấn. Tôi xin
được bày tỏ lòng biết ơn chân thành tới thầy.
Tôi xin chân thành cảm ơn Trường PTDT Nội trú cấp II-III Bắc
Quang, Hà Giang đã tạo mọi điều kiện thuận lợi nhất trong suốt quá trình
học tập, nghiên cứu và hoàn thành luận văn tốt nghiệp.
7Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
8
CHƯƠNG I
TỔNG QUAN
Như đã nêu trong phần mở đầu, hỗn hợp hai pha chất lỏng chứa bọt khí
hoặc hơi có tính chất đặc biệt, đó là sự kết hợp của các tính chất phi tuyến vật
lý mạnh, phân tán và hao tán năng lượng, nên biểu đồ mô tả các sóng có thể
có nhiều dạng. Chính vì vậy khi thay đổi các điều kiện thuỷ động lực, sẽ dẫn
đến sự thay đổi các cấu trúc về sóng và các quá trình tuơng tác giữa các pha.
Cho nên sự hiểu biết về các hiện tượng bất thường có thể xuất hiện khi sóng
áp suất lan truyền qua chất lỏng chứa bọt hơi là rất cần thiết để giải quyết một
loạt các bài toán thực tế như để phân tích chế độ làm việc quá độ của các thiết
bị năng lượng, phân tích tình huống hư hỏng và đảm bảo an toàn khi khai thác
của các nhà máy điện nguyên tử, phân tích các hiện tượng xâm thực trong các
máy tuốc bin, trong các hệ thống truyền dẫn thuỷ lực, để ứng dụng trong công
nghiệp khai thác, vận chuyển và chế biến dầu khí, trong công nghệ hoá học,
Dẫn đến sự cấp thiết của việc nghiên cứu thuỷ động lực học về các quá trình
sóng trong môi trường chất lỏng có bọt.
Môi trường hỗn hợp của chất lỏng với bọt của khí hoà tan và ngưng tụ
(hay khí không hoà tan và không ngưng tụ) thú vị ở chỗ trong chúng được cấu
thành từ ba yếu tố chính : tính phi tuyến, sự tán sắc và quá trình hao tán năng
lượng [1, 21]. Bức tranh sóng có thể có nhiều dạng và nó dễ dàng thay đổi
bằng cách thay đổi các điều kiện thuỷ động lực, cấu trúc và tính chất vật lý
nhiệt của hỗn hợp với các quá trình tương tác của các pha. Điều đặc biệt của
hỗn hợp chất lỏng chứa bọt trong các quá trình thuỷ động lực là sự xuất hiện
của năng lượng biến dạng cục bộ của hỗn hợp khi thay đổi thể tích môi
trường, sự thay đổi này chủ yếu do sự thay đổi thể tích của bọt trong hỗn hợp
do tính chất dễ co lại hay giãn nở của khí (hoặc hơi) trong bọt [3],[6].
Sự xuất hiện đồng thời của năng lượng biến dạng và sự đàn hồi dẫn đến
8Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
9
sóng có cấu trúc dao động bởi sự thay đổi thể tích của bọt. Ngoài ra, sự lan
truyền sóng áp suất trong những môi trường như vậy dẫn đến sự hoà tan hay
ngưng tụ của pha khí (hoặc hơi) và như vậy dẫn đến sự thay đổi cơ bản về bản
chất và cấu trúc vật lí của hỗn hợp.
Do sự thay đổi tính chất vật lí của hỗn hợp khi có sóng xung kích lan
truyền qua, nên trong hỗn hợp chất lỏng chứa bọt khí (hoặc hơi) thường xảy
ra những hiện tượng thể hiện tính chất phi tuyến của hỗn hợp đó là hiện tượng
khuếch đại hay tắt dần của sóng xung kích khi nó lan truyền trong hỗn hợp.
Về sự khuếch đại của sóng xung kích lan truyền trong chất lỏng chứa bọt khí
đã được đề cập đến trong [10], [11]. Trong các công trình này tác giả đã trình
bày hệ phương trình nhiệt thuỷ động lực học, giải hệ phương trình bằng
phương pháp số và đưa ra một số kết quả về sự ảnh hưởng của quá trình trao
đổi nhiệt – khối lượng lên động lực học sóng của môi trường hai pha lỏng -
khí hoặc hơi có xét đến khả năng ngưng tụ của hơi hay khí hoà tan của pha
khí [16]. Trong trường hợp không đi sâu nghiên cứu về cấu trúc của sóng
xung kích, mà chỉ xem sự truyền sóng xung kích như sự truyền của mặt gián
đoạn. Tác giả đã trình bày một số kết quả nghiên cứu về sự tăng áp suất của
sóng xung kích trong hỗn hợp chất lỏng chứa bọt khí hoặc hơi khi sóng này
tác động và bị phản xạ bởi một tường cứng, còn hiện tượng tắt dần có thể xem
trong [12].
Sóng xung kích trong hỗn hợp chứa bọt khí hoà tan hay ngưng tụ, về lý
thuyết đã được nghiên cứu trong các công trình [13, 18, 19], còn bằng thực
nghiệm đã được nghiên cứu bởi các công trình [12, 15, 17]. Nghiên cứu lý
thuyết về sự di chuyển của bọt và về sự co lại hay giãn nở của bọt trong dòng
chảy hai pha có thể xem trong [14]. Môi trường lỏng - bọt còn được sử dụng
trong công nghiệp năng lượng nguyên tử. Phân tích các chế độ làm việc quá
độ của của các thiết bị năng lượng, phân tích các tình huống hư hỏng và đảm
bảo an toàn khi khai thác các nhà máy điện nguyên tử đã được đề cập đến
9Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
10
trong công trình [4].
Nhận biết được những tính chất phi tuyến mạnh của môi trường hai pha
lỏng - bọt, mà trong đó có thể xuất hiện nhiều hiện tượng bất thường khi có
sóng xung kích lan truyền qua như hiện tượng tăng áp suất hoặc giảm đột
ngột của áp suất trong hỗn hợp [5, 19]. Hoặc một hiện tượng bất thường nữa
được xảy ra, đó là hiện tượng khuếch đại cường độ của sóng áp suất khi tác
động vào hỗn hợp lỏng hơi. Về vấn đề này đẫ được trình bày trong [9], trong
[9] tác giả đã chỉ ra rõ ràng rằng, trong giai đoạn đầu khi sóng xung kích tác
động vào hỗn hợp thì cường độ của sóng được tăng mạnh, tuy nhiên sự
khuếch đại này cong phụ thuộc vào hỗn hợp mà nó tác động. Trong tài liệu
[8], tác giả cũng đã đề cập tới quá trình trao đổi nhiệt và khối lượng trong quá
trình tương tác sóng. Tại đó cường độ của sóng tăng mạnh và cường độ của
sóng áp suất phụ thuộc vào tính chất vật lý nhiệt của hỗn hợp. Đây cũng là
một vấn đề rất được quan tâm. Do môi trường này lại được ứng dụng rộng rãi
trong rất nhiều ngành thuộc các lĩnh vực công nghiệp khác nhau. Chính vì
vậy, đây là một môi trường đã tập trung được rất nhiều nhà khoa học trong
nước và trên thế giới quan tâm. Với tổng quan tóm tắt trình bày thực trạng của
vấn đề đối với mô hình thủy - nhiệt mô tả dòng chảy hai pha lỏng – bọt, có
sóng xung kích lan truyền trong đó, tác giả đã cố gắng tìm hiểu và trình bày
một số nghiên cứu về lý thuyết cũng như về thực nghiệm của quá trình lan
truyền của sóng xung kích trong hỗn hợp chất lỏng hai pha, phân tích các hiện
tượng, các hiệu ứng bất thường có thể xảy ra trong hỗn hợp và một số nguyên
nhân cơ bản gây ra các hiệu ứng đó. Qua kết quả nghiên cứu của các công
trình đã cho sự hiểu biết đúng đắn về các hiện tượng có thể xuất hiện, để giải
quyết hàng loạt các bài toán thực tế như phân tích các chế độ làm việc của các
trạm năng lượng, phân tích các điều kiện hư hỏng và đảm bảo an toàn khi
khai thác các nhà máy điện nguyên tử, trong khai thác vận chuyển và chế biến
dầu khí, trong công nghệ hoá học và các quá trình tự nhiên…
10Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
11
Các kết quả nghiên cứu về lý thuyết và thực nghiệm của các công trình
được trình bày trên đây đã quan tâm đến hầu hết các vấn đề, các hiện tượng
có thể xuất hiện trong hỗn hợp lỏng - bọt khi tồn tại sóng xung kích truyền
qua như sự tăng áp suất hay giảm áp suất trong hỗn hợp. Đã tìm hiểu quá trình
chuyển động của bọt trong hỗn hợp, đến sự tách rời hay sự kết hợp lại của bọt
trong hỗn hợp vì chính bọt là nguyên nhân gây ra tính phi tuyến mạnh của
môi trường. Những vấn đề đặt ra trên đây có thể còn chưa được quan tâm
hoặc quan tâm chưa đúng mức.
Chính vì vậy, khi nhận biết được vấn đề này và hiểu được tầm quan trọng
của môi trường lỏng - bọt bởi nó được ứng dụng rất nhiều trong thực tế, nên
bản thân tôi đã cố gắng nghiên cứu, tìm hiểu và đã đưa ra được một số kết quả
nghiên cứu về quá trình lan truyền, tương tác của sóng xung kích trong hỗn
hợp lỏng - bọt đã nghiên cứu quá trình tương tác giữa sóng ngắn và sóng
ngắn, giữa sóng dài và sóng ngắn. Đã nghiên cứu, phân tích và tìm hiểu về sự
ảnh hưởng của các tham số đặc trưng như: các điều kiện ban đầu, các điều
kiện biên, các tính chất vật lý nhiệt của hỗn hợp, cường độ xung kích, phần
thể tích của pha hơi trong hỗn hợp lên quá trình lan truyền và tương tác kể
trên của sóng xung kích trong hỗn hợp lỏng - bọt. Các vấn đề đặt ra trên đây
cũng là nội dung chính sẽ được trình bày trong luận văn này.
11Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
12
CHƯƠNG II
CƠ SỞ LÝ THUYẾT
2.1. HỆ PHƯƠNG TRÌNH THUỶ- NHIỆT ĐỘNG LỰC HỌC CHO
MÔI TRƯỜNG HAI PHA LỎNG HƠI.
2.1.1. Xây dựng mô hình
Khảo sát hỗn hợp chất lỏng tĩnh chứa bọt hơi hình cầu được chứa trong
ống nằm ngang. Cùng với các giả thiết: độ nhớt và tính dẫn nhiệt chỉ quan
trọng trong quá trình tương tác giữa các pha, không quan trọng trong sự
truyền xung và năng lượng. Mật độ thực của chất lỏng chính là không thay
đổi = const. Trong quá trình tương tác không có sự phân chia bọt, bọt
không bị vỡ và vẫn có hình cầu.
Giả sử tồn tại hai sóng xung kích tác động vào hỗn hợp, lan truyền
trong hỗn hợp, tương tác nhau và lan truyền ra ngược nhau. Các sóng được sử
dụng là các sóng ngắn và sóng dài; sóng ngắn (còn được gọi là xung áp suất)
là sóng có quy luật của sự thay đổi nhanh áp suất, còn sóng dài là sóng có
cường độ luôn luôn không đổi. Nghiên cứu, phân tích và đánh giá các quy
luật phụ thuộc chung của các sóng xung kích trong các hỗn hợp chất lỏng
chứa bọt hơi trong hai quá trình: quá trình lan truyền, tương tác giữa sóng
ngắn và sóng ngắn, và quá trình lan truyền, tương tác giữa sóng ngắn và sóng
dài.
Quá trình sóng xung kích lan truyền vào, tương tác và lan truyền ra
được mô tả như sau:
a) Sóng xung kích lan truyền vào
p
p
12Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
13
b) Sóng xung kích tương tác nhau và lan truyền ra
2.1.2. Hệ phương trình cơ sở.
Hệ phương trình cơ học của hỗn hợp chất lỏng hai pha được xây dựng
trên cơ sở các định luật vật lý về bảo toàn khối lượng, xung lượng và năng
lượng bằng lý thuyết cơ học môi trường liên tục.[2]
2.1.2. 1. Hệ phương trình cơ học của hỗn hợp hai pha
Ta qui ước: chỉ số dưới i = 1 sẽ liên quan tới thông số pha chính ( chất
lỏng ), còn i = 2 liên quan tới pha phân tán ( pha hơi ). Để cho các phương
trình nhận được bớt phức tạp trong việc sử dụng phương pháp số ta giả thiết
như sau:
Các phần tử bọt của pha phân tán trong từng thể tích
V
tồn tại dạng
hình cầu có cùng bán kính R và mật độ thể tích của pha phân tán không quá
lớn :
1
2
2
Bỏ qua năng lượng và các ảnh hưởng khác của chuyển động hỗn loạn,
chuyển động quay, biến dạng của phần tử phân tán.
Bỏ qua tác dụng trực tiếp và sự va đập giữa các phần tử.
Không tồn tại quá trình chia nhỏ và sự hình thành các phần tử phân
tán mới.
Độ nhớt và tính dẫn nhiệt chỉ quan trọng trong quá trình tương tác giữa
các pha, không quan trọng trong sự truyền xung và năng lượng.
Vận tốc của các pha trong hỗn hợp là như nhau
.
21
vvv
13Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
14
Mật độ thực của chất lỏng chính là không thay đổi
const
0
1
.
Do giả thiết các phần tử của pha phân tán là hình cầu bán kính R, cho nên:
3
2 1 2
4
,1
3
Rn
.
Khi đó, mật độ pha
i
đặc trưng cho khối lượng pha trong một đơn vị thể
tích hỗn hợp được tính như sau:
,1
212
0
22
1
0
11
và mật độ của cả hỗn hợp :
2
0
21
0
121
Trong toạ độ Lagrange (
t,
), hỗn hợp chuyển động với vận tốc môi
trường v(
t,
), phương trình bảo toàn khối lượng pha, phương trình bảo toàn
số lượng bọt và phương trình bảo toàn xung lượng của cả hỗn hợp trong
trường hợp chuyển động liên tục một chiều có dạng sau: [10]
,4
2
0
11
nja
v
t
(2.1)
,4
2
0
22
nja
v
t
(2.2)
,4
3
4
20
2
3
jaa
t
(2.3)
,0
1
0
p
t
v
(2.4)
14Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
15
.0
0
vv
t
n
,
3
4
,,1,
3
22121
0
na
iii
.
2
2211
a
ppp
(2.5)
trong đó, t là thời gian, chỉ số i = 1, 2 liên quan tới tham số tương ứng của
chất lỏng và hơi nước;
0
,,,
iiii
p
là phần thể tích, áp suất, mật độ trung
bình và mật độ thực của pha thứ i ;
vp,,
là mật độ, áp suất và vận tốc
trung bình của hỗn hợp; n là số lượng bọt trong một đơn vị thể tích hỗn hợp;
j là cường độ của sự chuyển pha trong một đơn vị bề mặt diện tích giữa các
pha (j >0 tương ứng với trường hợp hoá hơi còn j < 0 ứng với trường hợp
ngưng tụ); là hệ số sức căng bề mặt. Từ đây về sau lấy đạo hàm riêng
t
khi cố định toạ độ Lagrăng có nghĩa là đạo hàm toàn phần theo thời gian.
Với giả thiết chất lỏng mang là không nén được, hơi nước tuân theo
phương trình trạng thái pha của khí hoàn hảo và tồn tại trạng thái bão hoà
trên bề mặt nó tuân theo phương trình Claperon-Clausiut chúng được cho
dưới dạng sau:
.1,,
0
1
0
2
0
2
2
2
2
2
0
22
0
1
l
T
dp
dT
TBpconst
(2.6)
trong đó T là nhiệt độ tuyệt đối, B là hằng số khí, l nhiệt dung riêng của hơi
nước bão hoà.
Trong hệ toạ độ Lagrăng - Ơle (
tr,,
) phương trình truyền nhiệt của chất
lỏng xung quanh bọt đối xứng cầu có dạng:
15Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
16
,
1
1
2
1
2
1
1
1
1
0
1
r
T
r
r
r
r
T
w
t
T
c
(2.7)
trong đó c
1
là nhiệt dung riêng của chất lỏng, w
1
là vận tốc hướng kính của
chuyển động nhỏ của chất lỏng xung quanh bọt thử nghiệm, được xác định
từ phương trình chuyển động nhỏ của chất lỏng xung quanh bọt thử nghiệm
cùng với tính không nén được của nó:
2
2
1
1
r
aw
w
a
. (2.8)
Điều kiện biên đối với phương trình (2.7) và (2.8) trong khuôn khổ
mô hình có thể được viết như sau:
r = a :
ar
aaa
r
T
qqqjlTT
1
112121
,,
,
.0:
1
3
1
2
r
T
ar
(2.9)
Trong trường hợp bọt hơi đồng đều như nhau thì:
,,,),(,),(
22
0
2
0
221
tppttTT
(2.10)
Dòng nhiệt
a
q
2
trong hơi nước từ giữa biên của các pha có thể nhận được
đường cong tích phân của phương trình dòng nhiệt của pha hơi từ 0 đến a
với tính toán tương ứng (3.5):
t
p
l
Tc
a
q
p
a
2
0
1
0
2
22
2
11
3
. (2.11)
trong đó c
p2
là nhiệt dung riêng của hơi khi áp suất không thay đổi.
16Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
17
Áp suất pha và kích cỡ của bọt liên quan tới điều kiện biến dạng
tương thích của mà nó được dùng trong phương trình tổng quát Rơlây-
Lamba (phương trình bảo toàn xung của chuyển động đối xứng cầu của chất
lỏng xung quanh bọt thử nghiệm ) với điều kiện ban đầu:
)1(
2
21
4
)1(5.1)1(
1
0
1
1
12
0
1
1
1
2
12
1
1
a
aa
a
jw
a
ppw
a
w
t
w
a
(2.12)
,)0()0(,
1
0
2
0
1
1 aa
wa
jj
w
t
a
2
22
3
1
2
2
2
2
3
1
2
1
0
1
3)2(
,
12
3
,)0(
aa
(2.13)
Hệ phương trình (2.1) đến (2.13) là hệ phương trình kín và khi với điều kiện
biên và điều kiện đầu thích hợp có thể sử dụng để khảo sát quá trình sóng
không dừng trong hỗn hợp chất lỏng chứa bọt hơi nước.
Để thuận lợi cho việc giải hệ phương trình trên bằng phương pháp số, ta có
thể biến đổi hệ phương trình này về dạng sau:
Kết hợp với (2.8), (2.9) ta viết phương trình bảo toàn khối lượng về dạng :
0
2
2
0
2
1
0
2
2
0
2
2
3
3
a
j
v
t
a
j
v
t
(2.14)
Đưa phương trìmh (2.2) về dạng :
v
j
att
0
222
0
2
0
2
2
3
(2.15)
17Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
18
Phương trình (2.3) có thể viết:
ja
t
a
a
t
a
20
2
2
0
2
3
44
3
4
kết hợp với (2.13):
0
1
1
j
w
t
a
a
, ta có:
0
1
0
2
1
0
2
0
2
3
j
wj
at
a
Thay vào (2.15), ta có:
v
j
at
j
wj
a
a
0
222
0
2
0
1
0
2
1
0
22
3
3
Suy ra:
v
j
w
at
a
0
2
0
1
0
2
1
0
2
2
0
2
2
3
1
(2.16)
với giả thiết:
1,,
2121
0
1
vvvconst
, ta có:
a
w
a
v
1
2
0
3
Do đó:
a
w
a
v
1
02
3
(2.17)
Thay phương trình (2.17) vào phương trình bảo toàn khối lượng hỗn hợp :
0
0
2
v
t
ta có:
18Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
19
aa
w
a
w
at
1
2
1
02
0
2
3
3
(2.18)
Lấy đạo hàm phương trình bảo toàn xung lượng hỗn hợp (2.4) theo
:
v
t
pp
pp
t
v
pp
t
v
p
t
v
0
2
0
2
2
0
0
2
0
2
2
0
2
0
2
2
0
2
0
11
0
11
0)
1
(
1
0)
1
(
Thế biểu thức (2.17) vào biểu thức trên, ta có:
a
w
at
pp
1
020
2
0
2
2
0
3
11
t
a
a
w
t
a
w
t
w
ata
w
w
t
a
a
t
w
a
w
at
w
at
pp
aaaa
a
a
aa
2
12
2
121
22
1
2
0
12
12
22
2
0
1
2
2
01
02
0
0
0
2
2
3
)(
)(
1
3
33
1
(2.19)
19Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
20
Phương trình (2.12) được viết lại như sau:
)1(
2
21
4
)1(
2
3
)1(
1
1
0
1
1
12
0
1
1
1
2
12
1
1
a
aa
a
jw
a
ppw
a
w
at
w
Biểu thức (2.16) sau khi sử dụng (2.17), có dạng:
aa
aa
w
a
j
a
w
a
w
a
j
w
at
1
0
2
22
0
1
2
1
2
1
02
0
2
0
1
0
2
1
0
2
2
0
2
2
333
3
3
1
0
1
11
2
12
0
1
1
2
33
j
w
a
w
j
w
a
aaa
(2.21)
Khi đó thế biểu thức (2.18), (2.21), (2.20) và
0
1
1
j
w
t
a
a
vào (2.19) và
biến đổi:
t
a
a
w
ta
w
t
w
ata
wpp
aaaa
2
12
2
121221
2
0
0
0
2
2
3
1
0
1
1
2
12
1
2
2
12
1
0
1
1
12
0
1
1
1
2
12
1
2
0
1
11
21
2
0
3
)1(
221
4
)1(
2
3
)1(
13
3
j
w
a
w
w
aa
w
jw
a
pp
w
a
w
aa
j
w
aa
w
a
a
a
a
a
aa
a
a
(2.20)
20Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
21
a
pp
jj
a
j
w
w
a
j
www
a
pp
jw
w
a
w
j
ww
a
a
a
aaa
a
aaaa
21
)1()1(
2
4
1
3
)1()1(3)1(
2
3
)1(3
)1(
3
)1()1(3
21
)1(
2
4
)1(
2
3
13
)1(
3
12
0
1
1
0
1
1
0
1
1
1
0
1
1
112211
2
1
1
2
2
2
0
1
0
1
11
2
112
12
0
1
1
0
1
1
1
1
2
12
0
1
1111
1
2
2
2
0
Đơn giản biểu thức này sẽ nhận được:
a
pp
a
w
w
a
pp
a
a
21
4
2
3
2
2
1
)1(
3
1
12
0
1
11
21
2
1
1
2
2
2
00
0
2
2
Hay phương trình này còn có thể biểu diễn dưới dạng:
pp
)(
2
2
(2.22)
trong đó:
0
0
1
)(
a
pp
a
w
w
a
a
a
21
4
2
3
2
2
1
)1(
3
12
0
1
11
21
2
1
1
2
2
2
0
Phương trình (2.22) có dạng elliptic mô tả sự truyền nhiễu áp suất trong hỗn
hợp hai pha lỏng khí khi chất lỏng là không nén được. Phương trình này cho
21Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
22
phép xác định được trường áp suất trong khoảng thời gian nhất định theo sự
thay đổi của các tham số cơ sở.
Để xác định áp suất pha hơi p
2
có thể sử dụng phương trình biểu diễn sự thay
đổi khối lượng của từng bọt (2.3) cùng với các điều kiện (2.7), (2.10) và
(2.12). Viết lại (2.3) dưới dạng:
ja
t
a
t
p
p
a
a
20
2
2
2
0
2
2
4)
3
(4
(2.23)
Từ phương trình trạng thái của pha hơi (2.6), ta có:
2
2
0
2
BT
p
đạo hàm hai vế của biểu thức này theo p
2
và sử dụng phương trình
Clapeyron-Clausius sẽ nhận được:
0
1
0
2
2
0
1
0
2
0
2
2
2
2
2
22
2
2
2
2
22
0
2
1
11
1
11
lBT
l
T
BT
p
BTp
T
BT
p
BTp
Thay vào (2.23) ta sẽ nhận được:
t
a
l
t
a
j
t
p
a
lBT
0
2
21
0
2
2
0
1
0
2
2
3
1
11
Sử dụng (2.11) ta nhận được biểu thức:
22Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
23
2
2
00
2
0
2 2 2 1
2
00
2 1 1
00
2
0
2 2 2 1
2
0 0 0
2 1 1 2
0
1
2
0
2
2
2
2
11
1 1 1
33
1 1 1
11
3
3
11
1
p
p
cT
pq
a a a
BT l t l l l t
cT
pq
aa
t BT l l l t l
q
a
tl
p
t
a
BT l
2
00
2
2
00
11
1
1
p
cT
ll
Bây giờ ta biến đổi biểu thức sau:
0
1
0
2
2
22
2
0
1
0
22
2
0
1
0
2
2
0
1
0
2
2
0
2
111
1
1
1
1
111
2
2
l
TcBT
l
BT
l
BT
p
l
Tc
llBT
p
p
Đặt
2
2
v
p
c
c
là tỷ số giữa nhiệt dung riêng của hơi khi áp suất không đổi và
thể tích không đổi. Mà
)1(
222
vvp
cBccB
, do đó biểu thức trên có
thể viết:
0
1
0
2
2
22
0
1
0
2
2
2
0
1
0
2
2
222
22
0
1
0
2
22
2
11
)1(
1
1
1
)1(
)1(
1
)1(
1
1
222
2
l
cT
l
cT
l
cT
p
l
TcTc
l
Tc
l
Tc
p
ppp
pv
vv
0
1
0
2
22
2
0
1
0
2
2
22
0
1
0
2
2
2
111)1(1
1
111)1(1
1
22
222
l
cT
l
cT
p
l
cT
l
cT
l
cT
p
pp
ppp
(2.2.4)
23Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
24
Đặt
0
1
0
2
22
111)1(1*
22
l
cT
l
cT
pp
Sử dụng điều kiện (2.9) vào biểu thức (2.24) ta nhận được phương trình sau:
t
a
r
T
a
p
t
p
a
1
0
2
122
*
3
(2.25)
Kết hợp phương trình truyền nhiệt của chất lỏng xung quanh bọt đối
xứng cầu với đạo hàm riêng, đặt
a
r
r *
thì có:
arr
a
a
r
rtt
t
rr
1
*
*
2
*
Khi đó phương trình truyền nhiệt (2.7) được viết lại như sau:
**
11
*
**
11
*
1
*
1
1
2
2
1
0
1
11
2
2
1
2
1
1
2
2
1
0
1
11
2
2
1
1
2
1
1
2
2
1
0
1
11
2
1
1
r
T
r
ra
rc
r
T
r
aw
a
ar
t
T
r
T
r
ra
rc
r
T
a
r
aw
r
T
a
ar
t
T
r
T
r
r
rc
r
T
r
w
t
T
a
a
a
mà
aa
w
jj
w
t
a
1
0
2
0
1
1
Do đó:
**
11
*
1
1
2
2
1
0
1
11
3
3
2
1
1
r
T
r
ra
rc
r
T
a
r
r
aw
t
T
a
(2.26)
24Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
25
Như vậy quá trình biến đổi ở trên đã có thể thu nhận được hệ phương trình
thuận tiện cho việc sử dung phương pháp số, đó là các phương trình: (2.6),
(2.12), (2.17), (2.18), (2.20), (2.22), (2.25), (2.26).
Điều kiện đầu và điều kiện biên
Nội dung của luận văn là nghiên cứu quá trình lan truyền và tương tác
sóng trong hỗn hợp hai pha được chứa trong ống nằm ngang với giả thiết tồn
tại hai sóng xung kích tác động vào hỗn hợp, lan truyền trong hỗn hợp,
tương tác nhau và lan truyền ra ngược nhau. Như vậy, điều kiện đầu và điều
kiện biên của bài toán được viết như sau:
a. Quá trình lan truyền, tương tác giữa sóng ngắn và sóng ngắn
1 0 2 0 0 1 2 0
0 1 2 0
2
0:
0
aa
t p p p p a a T T T
R
v v w w
(2.27)
0: ( )
: ( )
b
L
p p t
L p p t
b. Quá trình lan truyền, tương tác giữa sóng dài và sóng ngắn
1 0 2 0 0 1 2 0
0 1 2 0
2
0:
0
aa
t p p p p a a T T T
R
v v w w
(2.28)
0:
: ( )
e
b
pp
L p p t
Trong đó
p
e
là cương độ không đổi của sóng dài,
)(),( tptp
lb
là các xung
tam giác có thể biểu diễn như sau:
25Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
