Bài giảng Cơ lưu chất - Chương 3: Động lực học lưu chất - Lê Văn Dực
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
Chương 3: ĐỘNG LỰC HỌC LƯU CHẤT
) Lưu chất có thể chuyển động bên trong biên rắn (trong kênh, trong ống, …) hoặc bao quanh
biên rắn (máy bay, tàu thuyền…), và được xem như một môi trường liên tục. Mỗi phần tử
lưu chất được xem như có kích thước (một điểm) và có khối lượng vô cùng nhỏ.
) Chương này sẽ đề cập đến hai phần: động học và động lực học. Trong phần động học,
những thông số của dòng chảy được quan tâm xem xét là vận tốc, gia tốc và sự biến thiên
của chúng theo thời gian. Bên cạnh đó, phương trình liên tục của lưu chất cũng được trình
bày ở nhiều dạng khác nhau.
) Trong phần động lực học, cơ sở lý thuyết sự chuyển động của phần tử lưu chất có xét đến
nguyên nhân gây ra chuyển động (lực) và việc xây dựng các phương trình vi phân chuyển
động được trình bày. Các nguyên lý biến thiên động lượng và bảo toàn năng lượng được áp
dụng để xây dựng các phương trình cơ bản động lực học như Euler, Navier-Stoke, phương
trình năng lượng và động lượng được đề cập. Sau đó việc áp dụng các phương trình này cho
đoạn dòng chảy lưu chất trọng lực, không nén được, chuyển động ổn định được trình bày.
3.1 Hai phương pháp mô tả chuyển đông của lưu chất :
3.1.1 Phương pháp Lagrange :
+ Xét một hệ thống trục tọa độ cố định OXYZ (hình H.3.0a).
+ Chuyển động của lưu chất được mô tả bằng vị trí của các phần tử lưu chất theo thời gian.
z
y
Mo(xo,yo,zo)
M(x,y,z)
r
r
ro
r
O
x
H.3.0a
r
t = 0
⇒
⇒
r (xo , yo, zo)
o
r
r
r
r
∀ t
r (x, y, z)
→
r = f (r , t)
(3.1)
o
hay
x = x(xo , yo, zo, t)
y = y(xo , yo, zo, t)
z = z(xo , yo, zo, t)
35
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
→ xo , yo, zo và t được gọi là những biến số Lagrange.
r
r
r
+ Nếu như biết r , ta có thể tính được vận tốc u (ux, uy, uz) và gia tốc a (ax, ay, az) của phần tử lưu
chất tại thời điểm nào đó như sau :
r
r
r
dr
r
d 2 r
u =
&
a =
dt 2
dt
dy
dt
dx
dz
dt
ux =
ax =
,
uy =
,
uz =
az =
(3.2.a)
(3.2.b)
dt
d 2 x
dt 2
d2 y
dt2
d 2 z
dt 2
,
ay =
,
+ Đường nối vị trí của phần tử lưu chất theo thời gian được gọi là quỹ đạo của chuyển động của
phần tử lưu chất đó. Phương trình (3.2a) chính là phương trình quỹ đạo của phần tử lưu chất.
+ Nhận xét :
- Phương pháp này được dùng rộng rãi trong cơ học chất rắn, rất thuận tiện nếu như số lượng
phần tử ít.
- Trong cơ học lưu chất phương pháp này khó thực hiện vì số lượng phần tử rất lớn.
- Ngòai ra, trong lưu chất, tồn tại hiện tượng khuếch tán phân tử → các phân tử lưu chất chỉ giữ
được đặc tính riêng trong một khoảng thời gian rất ngắn → việc xác định quỹ đạo rất khó khăn.
- Trong thực tế việc nghiên cứu từng phần tử riêng lẽ là không cần thiết vì việc giải các bài toán
này rất phức tạp nên phương pháp này ít được dùng.
- Phương pháp Lagrange chủ yếu được áp dụng trong một số trường hợp như để nghiên cứu hiện
tượng sóng biển, quan sát vết sau vật.
3.1.2 Phương pháp Euler :
+ Trong phương pháp này, các thông số động học của các phần tử chất lỏng cùng đi qua một điểm
nào đó trong miền chuyển động được quan tâm. Ví dụ trạm đo vận tốc đặt cố định trên sông.
+ Trong một hệ tọa độ xác định, chuyển động của lưu chất được mô tả bằng vận tốc của các phần
tử lưu chất tại mỗi vị trí theo thời gian.
+ Tại một điểm M(x,y,z) cố định trong không gian, ở thời điểm t có một phần tử lưu chất qua M có
r
r
vận tốc u (M(x,y,z), t). Ở thời điểm t + dt có một phần tử khác đi qua M có vận tốc u (M(x,y,z),
36
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
t+dt). Như vậy trong toàn môi trường lưu chất (∀ M ∈ miền chuyển động), ta xác định được
r
r
r
r
trường các vectơ vận tốc u = u (M, t), hay u = u (x,y,z,t), với các thành phần :
ux =ux(x,y,z,t)
uy =uy(x,y,z,t)
(3.3)
uz =uz(x,y,z,t)
x, y, z và t được gọi là các biến Euler.
+ Nhận xét :
- Phương pháp Euler đặc biệt có ưu thế trong việc tìm lời giải của các bài toán lien quan đến
chuyển động của lưu chất nhờ vào lý thuyết trường → có thể xác định các đặc trưng chuyển
động tại những vị trí cần thiết (ví dụ : dòng chảy bao quanh vật, xác định vận tốc, áp suất tại các
điểm trên bề mặt vật rắn).
- Ngòai ra, vì số lượng các phương trình vi phân hữu hạn và có thể được biến thành các phương
trình sai phân tuyến tính, nên có thể tìm được lời giải bằng phương pháp gần đúng, chẳng hạn
phương pháp sai phân và được giải số thông qua chương trình máy tính.
3.2 Một số khái niệm thường dùng :
1. Đường dòng (lưu tuyến) :
+ Đường dòng tại một thời điểm t, là đường
vạch ra trong lưu chất sao cho các vectơ vận
tốc của các phần tử lưu chất trên đường đó
đều tiếp xúc với nó.
+ Đường dòng được xác định từ vectơ vận tốc
của các phần tử lưu chất khác nhau tại một
thời điểm xác định (≠ quĩ đạo là vết của
một phần tử lưu chất theo thời gian).
+ Theo định nghĩa, tiếp tuyến của đường dòng
tại mỗi điểm trùng với vectơ vận tốc tại
r
r
r
r
r
điểm đó → u // d r hay u x d r = 0
Ở đây :
r
r
r
r
u = ux i + uy j + uz k
37
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
r
r
r
r
d r = dx. i + dy. j + dz. k
→ Phương trình vi phân của đường dòng :
dx dy dz
=
=
(3.4)
ux
uy
uz
+ Nhận xét :
- Hai đường dòng trong cùng một thời điểm không thể cắt nhau hoặc tiếp xúc nhau. Điều này
có thể được chứng minh bằng phương pháp phản chứng: ở một điểm, ở một thời điểm chỉ có
r
một véctơ vận tốc u duy nhất (hình H.3.0b) , hoặc nếu chúng tiếp xúc nhau thì sẽ không bảo
đảm được phương trình bảo tòan khối lượng (hình H3.0.c); hoặc có thể được cứung minh
như trong chương thế lưu, đối với bài toán phẳng (chuyển động 2 chiều), trên mỗi đường
dòng, hàm dòng ψ = const và lưu lượng Q chảy giữa hai đường dòng bằng (ψ1 - ψ2) → ở
một điểm không thể có hai giá trị ψ khác nhau. Tuy nhiên, ngọai trừ trường hợp tại các điểm
kỳ dị trong trường chuyển động. Ví dụ điểm dừng trong chuyển động đều quanh hình trụ
tròn đứng yên (H.3.0d).
(C2)
(C1)
M
r
r
u2
u1
(C1)
M
(C2)
H.3.0c
H.3.0b
H.3.0d
38
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
- Trong chuyển động ổn định, đường dòng trùng với quỹ đạo và hình dạng các đường dòng
không thay đổi theo thời gian.
2. Dòng nguyên tố (sợi lưu tuyến) và dòng chảy:
H.3.1a
H.3.1b
+ Dòng nguyên tố là tập hợp tất cả các đường dòng đi qua các điểm trên dA (vi phân diện tích
trong không gian chứa lưu chất chuyển động).
+ Dòng chảy là tập hợp tất cả các dòng nguyên tố đi qua diện tích A hữu hạn.
+ Ống dòng là bề mặt bao quanh dòng chảy ⇔ tập hợp tất cả đường dòng đi qua các điểm trên
chu vi bề mặt của A. Vì vận tốc của các phần tử lưu chất trên ống dòng luôn theo phương tiếp
tuyến nên lưu chất không thể đi xuyên qua ống dòng.
3. Mặt cắt ướt, chu vi ướt, bán kính thủy lực :
m = cotang(α)
; A=(b+m.h).h
P = b + 2.h 1+ m2 ; R = A/P
H.3.1d
H.3.1c
39
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
H.3.1e
+ Mặt cắt ướt A là mặt cắt thẳng góc với tất cả các đường dòng qua nó. Mặt cắt ướt có thể là mặt
phẳng (khi tất cả các đường dòng song song và thẳng), hoặc là mặt cong (H.3.1c & d).
+ Chu vi ướt χ (hoặc P) là phần chu vi của mặt cắt ướt tiếp xúc với thành rắn (H.3.1c & d).
A
+ Bán kính thủy lực (ướt) R =
χ
4. Lưu lượng:
Lưu lượng là lượng lưu chất đi qua mặt cắt ướt trong một đơn vị
thời gian:
H.3.1f
-
Lưu lượng thể tích :
Q = u dA.
(3.5a)
∫∫
A
Với A là diện tích ướt, và dA là vi phân diện tích trên A và u là vận tốc (pháp tuyến với dA) tại một
điểm trên dA.
Thứ nguyên của [Q] = L3/T, đơn vị là m3/s, l/s.
-
Lưu lượng khối lượng :
∫∫ρ u d.A.
QM =
(3.5b)
A
Thứ nguyên của [Q M] = M/T, đơn vị là kg/s, g/s.
5. Vận tốc trung bình mặt cắt ướt :
40
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
+ Vận tốc trung bình V đối với mặt cắt ướt A:
Q
A
1
V =
=
u.dA
(3.5c)
A ∫∫
A
+ Vận tốc trung bình V là giá trị tưởng tượng, khác với vận tốc thực u tại mỗi điểm.
+ Trong thực tế, dễ xác định vận tốc trung bình V bằng cách đo lưu lượng Q và diện tích ướt A.
3.3 Phân loại chuyển đông :
3.3.1 Phân loại theo ma sát nhớt :
+ Chuyển động lưu chất lý tưởng:
Trong thực tế mọi lưu chất đều có tính nhớt. Tuy nhiên, khi xét chuyển động của lưu chất mà
tính nhớt hay lực ma sát nhớt cản trở chuyển động của các phần tử lưu chất không đáng kể,
người ta có thể bỏ qua lực ma sát nhớt (để làm đơn giản bài toán) và khi đó lưu chất được xem
là lý tưởng. Ví dụ: chuyển động lưu chất qua biên rắn ⇒ lớp lưu chất sát biên: lưu chất thực; lớp
lưu chất xa biên: lưu chất lý tưởng (bỏ qua ma sát).
+ Chuyển động của lưu chất thực:
- Có hai loại : chuyển động tầng và chuyển động rối.
- Osborne Reynolds (1883) là người đầu tiên đã làm thí nghiệm để chứng minh sự tồn tại của
hai loại chuyển động này.
- Mô tả thiết bị thí nghiệm :
Thiết bị thí nghiệm gồm một bình lớn đựng nước, một bình nhỏ đựng nước màu (có khối
lượng riêng bằng khối lượng riêng của nước). Nước chảy ra ngoài nhờ ống thủy tinh có van
điều chỉnh ở cuối ống. Dòng nước màu cũng cho chảy vào trong ống (xem hình H.3.2).
- Mô tả hiện tượng :
muc mau
nuoc
dong muc mau
H. 3.2a Thiết bị thí nghiệm Reynolds
H.3.2b Các trạng thái chảy
41
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
9 Khi vận tốc nước trong ống nhỏ, dòng màu chảy thành sợi chỉ thẳng tắp ⇒ lưu tốc song
song thành lớp → Trạng thái chuyển động tầng.
9 Khi vận tốc nước tăng lên (mở van), dòng màu bắt đầu gợn sóng và sau đó hòa lẫn vào trong
nước ⇒ các phần tử chuyển động hỗn độn, trộn lẫn nhau, vận tốc dao động liên tục theo mọi
phương và thay đổi trị số ⇒ áp suất dao động → Trạng thái chuyển động rối.
9 Khi chuyển động của phần tử lưu chất bị nhiễu ⇒ lực quán tính có khuynh hướng làm đổi
hướng vận tốc, lực ma sát nhớt do các phần tử bao quanh có xu hướng giữ nó ở trạng thái cũ
⇒
h Trong chuyển động tầng: lực quán tính nhỏ, lực ma sát nhớt quan trọng đủ để giữ các
phần tử lưu chất theo một trật tự nhất định.
h Trong chuyển động rối: lực quán tính lớn, lực ma sát không còn đủ sức khống chế các
nhiễu.
- Tiêu chuẩn phân loại trạng thái chảy :
Tiêu chuẩn để phân biệt trạng thái chảy là số Reynolds (Re), số không thứ nguyên, là tỷ số
giữa lực lực quán tính và lực nhớt:
ρ l2V2
μ Vl
ρ l V
μ
ρ lV V.l
const.
=const
→
Re =
=
μ
ν
Với :
ρ : khối lượng riêng của lưu chất (kg/m3),
l : chiều dài đặc trưng của dòng chảy (m),
V : vận tốc đặc trưng của dòng chảy (m/s),
μ : hệ số nhớt động lực học của lưu chất (Pa.s),
μ
ρ
ν =
là hệ số nhớt động học của lưu chất (m2/s).
# Khi Re < Re tới hạn (Re th) ⇒ chuyển động tầng (hình H.3.2b)
# Khi Re > Re tới hạn (Re th) ⇒ chuyển động rối (hình H.3.2b)
- Dòng chảy trong ống tròn: Đối với dòng chảy trong ống tròn, người ta lấy chiều dài đặc
trưng l = D. Và số Reynold tới hạn như sau:
42
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
ρ.V.D V.D
Reth =
=
= 2300
(3.6)
μ
ν
3.3.2 Phân loại theo thời gian :
+ Chuyển động ổn định :
Chuyển động ổn định (thường trực, dừng) khi các tính chất của chuyển động không thay đổi
theo thời gian ⇒ đạo hàm riêng phần theo thời gian của các đại lượng mô tả chuyển động bằng
∂*
0 → (
= 0) →
∂t
∂u
∂t
∂a
∂t
u = u(x, y, z) & a = a(x, y, z)
⇔
= 0
&
= 0
+ Chuyển động không ổn định :
Chuyển động không ổn định (không thường trực, không dừng) khi các tính chất của chuyển
động thay đổi theo thời gian →
u = u(x, y, z, t)
&
a = a(x, y, z, t)
+ Chuyển động rối :
- Chuyển động rối thực chất là chuyển động không ổn định và rất phức tạp.
- Thường người ta thay thế chuyển động rối thực bởi dòng chảy rối trung bình thời gian. Vận
tốc chuyển động của dòng chảy rối thực được xem là tổng của 2 thành phần : vận tốc trung
bình thời gian ⎯u và thành phần vận tốc mạch động u’ :
ux = ux + u'x
Với, ux : vận tốc trung bình theo thời gian theo phương x.
u’x : vận tốc mạch động theo phương x
t+T
1
ux =
ux.dt
∫
T
t
# Nếu ⎯u không phụ thuộc thời gian ⇒ chuyển động rối ổn định (H.3.3b).
# Nếu ⎯u phụ thuộc thời gian ⇒ chuyển động rối không ổn định (H.3.3b).
43
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
ux
u'x
ux
t
T
t
t+T
Hình 3.3a Vận tốc tức thời ux
3.3.3 Phân loại theo không gian :
+ Chuyển động 3 chiều :
H.3.3b
- Mọi chuyển động trong thực tế đều xảy ra trong không gian 3 chiều ⇒ vận tốc, áp suất và
các yếu tố khác của dòng chảy thay đổi theo các phương khác nhau.
- Trong một số trường hợp đặc biệt, ta có thể đơn giản hóa chuyển động dòng chảy thành dòng
chảy 2 hoặc 1 chiều, nếu như sự sai số do đơn giản là không đáng kể.
+ Chuyển động 2 chiều :
Hình 3.4 Chuyển động 2 chiều qua đập tràn
Khi các thông số chuyển động thay đổi theo 2 phương vuông góc nằm trong một mặt phẳng và xem
như không thay đổi trong các mặt phẳng song song với nó. Ví dụ dòng chảy qua bờ tràn có mặt cắt
ngang không đổi, bề rộng dài vô tận.
+ Chuyển động 1 chiều :
- Khi các thông số của chuyển động như vận tốc, áp suất, độ cao mô tả chuyển động tại một thời
điểm nhất định chỉ thay đổi theo 1 phương và có giá trị như nhau tại mọi điểm trên cùng một
mặt cắt ướt.
- Chiều chuyển động được chọn thường dựa theo chiều của đường dòng trung tâm có giá trị vận
tốc, áp suất và độ cao là giá trị trung bình trên mặt cắt ướt thẳng góc với đường dòng tại điểm
ấy.
44
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
- Ví dụ : chuyển động của lưu chất trong ống tròn khi dùng các yếu tố trung bình để diễn tả. Thực
ra vận tốc bằng không ở sát thành rắn và cực đại ở tâm ống. Do đó, khi dùng vận tốc trung bình
để diễn tả dòng chảy, người ta đưa thêm các hệ số điều chỉnh khi cần thiết (tăng độ chính xác)
trong các phương trình năng lượng và động lượng.
+ Dòng đều :
- Dòng chảy 1 chiều, có sự phân bố vận tốc trên mặt cắt ướt dọc dòng chảy không thay đổi.
- Điều kiện để có dòng đều :
# Điều kiện về đường dẫn (tiết diện, độ nhám, độ dốc …): không đổi dọc dòng chảy.
# Điều kiện về lưu lượng: không đổi dọc dòng chảy.
Hình 3.5 Biểu đồ phân bố vận tốc trong ống tròn
3.3.4 Phân loại theo tính nén được :
+ Chuyển động của lưu chất không nén được :
Khi lưu chất chuyển động nếu khối lượng riêng của lưu chất phụ thuộc rất ít vào áp suất và
nhiệt độ. Ví dụ chất lỏng hoặc chất khí với vận tốc nhỏ có số Mach M < 0,3.
+ Chuyển động của lưu chất nén được :
Khi lưu chất chuyển động nếu khối lượng riêng của lưu chất phụ thuộc vào áp suất và nhiệt
độ. Ví dụ chất khí với vận tốc lớn, có số Mach M > 0,3.
+ Tiêu chuẩn phân loại chuyển động :
Số Mach (M= u/a) là tỷ số giữa vận tốc chuyển động của lưu chất (u) và vận tốc truyền âm
trong lưu chất (a) được dùng để phân loại chuyển động của chất khí.
# M < 0,3 ⇒ ρ thay đổi không đáng kể → Lưu chất không nén được.
# M > 0,3 ⇒ ρ thay đổi đáng kể ⇒ ảnh hưởng của tính nén được phải được xét đến →
Lưu chất nén được.
45
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
• Nếu M < 1 : chuyển động dưới âm thanh (subsonic)
• Nếu M=1 : chuyển động ngang âm thanh (sonic)
• Nếu M > 1 : chuyển động trên âm thanh (supersonic)
• Nếu M >> 1 : chuyển động siêu thanh (hypersonic)
3.4 Gia tốc của phần tử lưu chất :
+ Theo phương pháp Euler :
Vận tốc tại một điểm của phần tử lưu chất là một hàm số theo tọa độ điểm và thời gian :
r
r
r
u = u (x,y,z,t). Tại thời điểm t, tại điểm M(x,y,z), lưu chất có vận tốc u . Ở thời điểm t + Δt; phần
tử lưu chất ở vị trí M1(x+Δx, y+Δy, z+Δz) và có vận tốc là:
r
u1
y1
r
Δy
u
y
M
Δx
x
x1
Hình 3.6 Sơ đồ tính gia tốc – phương pháp Euler
r
r
u1 = u ( x+Δx, y+Δy, z+Δz, t+Δt).
Xét theo phương x :
u1x = ux(x+Δx, y+Δy, z+Δz, t+Δt)
Giả sử hàm ux là hàm liên tục theo 4 biến số x, y, z, t. Khai triển chuỗi Taylor bậc 1, ta có :
∂ux
∂y
∂ux
∂z
∂ux
∂t
∂ux
∂x
u1x = ux +
Δx+
Δy+
Δz+
Δt
Khi đó gia tốc của phần tử lưu chất được xác định như sau :
dux
dt
u1x − ux
Δt
ax =
= limΔt→0
46
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
∂ux
∂ux
∂ux
Δy
∂ux Δx
Δz
= limΔt→0
(
.
+
.
+
.
+
)
∂x Δt ∂y Δt ∂z Δt ∂t
∂ux ∂ux ∂ux
∂ux dx
dy
dz
ax =
ax =
+
.
+
.
+
∂x dt ∂y dt
∂z dt ∂t
∂ux
∂y
∂ux
∂z
∂ux
dux ∂ux
=
ux+
.uy+
.uz+
(3.7a)
dt
∂x
∂t
Tương tự :
duy ∂uy
∂uy
∂y
∂uy
∂uy
∂t
ay =
=
.ux +
ux +
.uy +
.uz +
(3.7b)
(3.7c)
dt
∂x
∂z
duz ∂uz
∂uz
∂y
∂uz
∂z
∂uz
∂t
az =
=
uy +
uz +
dt
∂x
gia tốc đối lưu
gia tốc cục bộ
Theo (3.7), ta thấy gia tốc của phần tử lưu chất gồm 2 thành phần :
- Thành phần gia tốc đối lưu :
∂uz
∂x
∂uz
∂y
∂uz
∂z
ux +
uy +
uz
biểu thị suất biến thiên của vận tốc của các phần tử lưu chất do sự khác biệt vị trí trong trường
chuyển động.
- Thành phần gia tốc cục bộ :
∂uz
∂t
biểu thị suất biến thiên cục bộ của vận tốc do chuyển động không ổn định. Gia tốc cục bộ bằng
0 khi chuyển động ổn định.
Phương trình (3.7) có thể viết dưới dạng vectơ như sau :
r
r
r du ∂u
a =
r
r r
=
+u .( grad.u )
(3.8)
dt
∂t
47
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
Với toán tử :
r
k
r
r
i +
r
j +
r
∂
∂x
∂
∂y
∂
∂z
grad =∇ =
+ Theo phương pháp Lagrange :
r
r
Vectơ vận tốc là hàm số của tọa độ ban đầu u = u (xo , yo, zo, t). Tọa độ ban đầu là các hằng số nên
đạo hàm toàn phần hay gia tốc của phần tử lưu chất cũng chính là đạo hàm riêng phần theo thời
gian, không có thành phần gia tốc đối lưu.
r
r
r du ∂u
a =
=
(3.9)
dt ∂t
3.5 Phương pháp thể tích kiểm soát – Lý thuyết vận tải Reynolds:
H.3.7
9 Thể tích kiểm soát (control volume): là một thể tích cố định, W đối với một hệ trục tọa độ
OXYZ trong trường lưu chất chuyển động, nằm trong mặt cong kín S, được gọi là mặt kiểm
sóat (control surface), qua đó sự biến thiên của thông số (X, extensive parameter) dòng chảy
theo thời gian được nghiên cứu.
9 Trong dòng chảy xét thể tích kiểm sóat W được bao bọc bởi mặt kiểm sóat S (hình H.3.7).
Trong đó đại lượng X là hệ thống được nghiên cứu (khối lượng, năng lượng hoặc động
lượng…). Và, gọi k là giá trị của đại lượng X chứa trong một đơn vị khối lượng (k, intensive
parameter).
48
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
Ta có công thức để tính X như sau:
k.ρ.dw
(3.10a)
X = ∫∫∫
W
- Tại thời điểm t1 ⇒ thể tích W1 → chứa X1.
- Tại thời điểm t2 = t1 + Δt ⇒ thể tích W1 di chuyển đến thể tích W2 → chứa X2.
- Gọi A là thể tích mà: A ⊂ W1 và A ⊄ W2
- Gọi B là thể tích mà: B = W1 ∩ W2
- Gọi C là thể tích mà: C ⊂ W2 và C ⊄ W1
) Như vậy có nghĩa rằng, sau thời gian Δt:
- Một số phần tử lưu chất di chuyển ra khỏi thể tích kiểm soát và chiếm thể tích là C.
- Một số phần tử lưu chất bên ngoài vào và lắp đầy thể tích A
) Như vậy ở thời điểm t1 đại lượng X của toàn hệ thống nằm trong thể tích kiểm soát là:
Xt1 =XAt1 + XBt1
(a)
(b)
(c)
Và ở thời điểm t2 ⇒
Xt2 =XBt2 + XCt2
Sự thay đổi đại lượng X của hệ thống trong khoảng thời gian Δt là :
ΔX = (XBt2 + XCt2)-( XAt1 + XBt1)
Thêm vào và trừ XAt2 , ta được:
ΔX = (XBt2 + XAt2)-( XAt1 + XBt1)+ XCt2 - XAt2
= X(A+B)t2 - X(A+B)t1 + XCt2 - XAt2
Sự biến thiên của X theo thời gian là :
X( A+B)t − X( A+B)t
XCt − X At
ΔX
Δt
2
1
2
2
=
+
Δt
Δt
49
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
X ( A+ B )t − X ( A+ B )t
X Ct − X At
dX
dt
ΔX
Δt
2
1
2
2
= lim
= lim
+ lim
Δt→0
(d )
Δt→0
Δt→0
Δt
(1)
Δt
(2)
∂X
∂t
∂
∂t
⎛
⎜
⎞
⎟
(1) =
thể tích kiểm soát =
k.ρ.dw
(e)
∫∫∫
w
⎝
⎠
H.3.8
Với: XCt2 là giá trị của đại lượng X của các phần tử lưu chất ra khỏi thể tích kiểm soát thông
qua diện tích S2.
XAt2 là giá trị của đại lượng X của các phần tử lưu chất đi vào thể tích kiểm soát thông qua diện
tích S1.
́
r
r
- Xét diện tích dA ⊂ S2 có vectơ chỉ phương làn , lưu tốc qua dA là u . Thể tích qua dA trong
đơn vị thời gian là :
r r
r r
dQ = un.dA = u.n.dA → dQM = ρ.un.dA = ρ.u.n.dA
r r
Do đó, giá trị đại lượng X của phần tử lưu chất đi qua dA là : k.ρ.u.n.dA
→
r r
XCt2 = Δt.
k.ρ.u.n.dA
(f)
∫∫
S2
Và tương tự:
At2 = -Δt.
r r
X
k.ρ.u.n.dA
(g)
∫∫
S1
r
r r
(XAt2 mang dấu trừ vì lưu chất đi vào mặt S1, n hướng ra ⇒ u.n luôn luôn âm)
50
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
Từ đó ta suy ra:
r r
r r
Δt. k.ρ.u.n.dA +Δt. k.ρ.u.n.dA
X
− X
Δt
∫∫
∫∫
Ct2
At2
S2
S1
(2) = lim
= lim
Δt→0
Δt→0
Δt
r r
r r
=
k.ρ.u.n.dA= k.ρ.u.n.dA
∫∫
∫∫
S1+2
S
Do đó, phương trình (d) trở thành :
dX
dt
∂X
∂t
r r
=
+
k.ρ.u.n.dA
(3.10b)
hệ thống
∫∫
CV
S
Với: CV : thể tích kiểm soát
S : diện tích kiểm sóat.
Hay :
∂X
∂t
r r
k.ρ.u.n.d)A
ddt ∫∫∫
∫∫
hệ thống (k.ρ.dw) =
+
c
0 (
1
CV
S
3.6 Phương trình liên tục:
+ Theo định luật bảo toàn khối lượng: khối lượng của hệ thống không thay đổi theo thời gian
dm
= 0
dt
Ta xét X = m =
k.ρ.dw =
ρ.dw ⇒ k=1
∫∫∫
∫∫∫
w
w
dm
∂m
∂t
r r
∫∫ ρ.u.n.dA
→
HT
=
+
CV
S
dt
∂ρ
∂t
r r
→
.dw +
0
(
)
1
=
1
∫∫∫
∫∫ ρ.u.n.dA
S
w
Đây là dạng tích phân của phương trình liên tục khối lượng.
+ Dùng phép biến đổi GAUSS để chuyển từ tích phân mặt sang tích phân thể tích:
r r
r
div(ρ.u).dw
∫∫ ρ.u.n.dA = ∫∫∫
S
w
51
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
∂ρ
∂t
r
→
.dw +
div(ρ.u).dw= 0
∫∫∫
∫∫∫
w
w
∂ρ
∂t
r
[
+ div(ρ.u )].dw = 0
∫∫∫
w
Tích phân này áp dụng cho một thể tích W bất kỳ ⇒ nên ta có thể chuyển sang dạng vi phân của
phương trình liên tục.
3.6.1 Phương trình liên tục dạng vi phân, lưu chất nén được chuyển động không ổn định:
∂ρ
∂t
r
[
+ div(ρ.u )] = 0
(3.12a)
Nếu là chuyển động ổn định ta có:
r
div(ρ.u ) = 0
(3.12b)
3.6.2 Phương trình liên tục dạng vi phân, lưu chất không nén được, chuyển động ổn định:
Vì ρ=const, nên từ (3.12b), ta suy ra phương trình liên tục có dạng :
r
div(u ) = 0
(3.13)
(3.14)
- Trong hệ thống tọa độ Descartes (3.13) có dạng :
∂uy
∂y
∂ux
∂x
∂uz
∂z
+
+
= 0
- Trong hệ thống tọa độ trụ (3.13) có dạng :
∂u
1
r ∂θ
∂uz
∂z
1 ∂
r ∂r
θ
(r.ur ) +
+
= 0
(3.15)
3.6.3 Phương trình liên tục, lưu chất không nén được, chuyển động ổn định, 1 chiều:
Vì chuyển động ổn định, nên từ (3.11), ta suy ra:
r r
ρ.u.n.dA=0
(3.11a)
∫∫
S
H.3.9
52
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
Xét thể tích kiểm soát là một đoạn dòng chảy giới hạn bởi hai mặt cắt ướt A1 và A2 và diện tích bao
quanh ống dòng Sb ⇒
S = A1 + A2 + Sb
r r
r r
ρ.u.n.dA+
r r
ρ.u.n.dA+
r r
ρ.u.n.dA= 0
→
ρ.u.n.dA=
∫∫
∫∫A
∫∫A
∫∫
S
Sb
1
2
r r
Trên Sb ⇒ u.n = 0
r r
Trên A1 ⇒ u.n = -u1
r r
Trên A2 ⇒ u.n = u2
→
-
ρ.u1.dA +
ρ.u1.dA =
ρ.u2.dA = 0
ρ.u2.dA
∫∫A
∫∫A
1
2
∫∫A
∫∫A
1
2
Qm1 = Qm2
ρ1.V1.A1 = ρ2.V2.A2 = const
(3.16a)
(3.16b)
Với ρ1 và V1 là khối lượng riêng và vận tốc trung bình của lưu chất tại mặt cắt ướt 1-1, tương tự
ρ2 và V2 áp dụng cho mặt cắt ướt 2-2.
+Kết luận :
- Đối với lưu chất chuyển động ổn định, 1 chiều lưu lượng khối qua các mặt cắt ướt đều bằng
nhau.
- Khi lưu chất không nén được ρ1 = ρ2 = const ⇒ khi đó phương trình liên tục viết cho đoạn dòng
chảy có dạng :
Q1 = Q2
(3.17a)
(3.17b)
Hay V1A1 =V2.A2
53
Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực
3.6.4 Phương trình liên tục lưu chất không nén được, chuyển động ổn định tại nút:
Q2
Q3
Q4
Q1
Q5
Q1 = Q2 + Q3
Q1 + Q2 + Q5 = Q3 + Q4
H.3.18b
H.3.18a
Mở rộng phương trình (3.17a) đối với dòng chảy một chiều, chuyển động ổn định, lưu chất không
nén được, xét tại nút I, có nhiều nhánh vào và ra, phương trình liên tục, trở thành:
Q
=
Q
(3.18)
∑
∑
vào
ra
3.7 Phương trình vi phân chuyển động của lưu chất lý tưởng (phương trình Euler) :
3.7.1 Phương trình Euler :
H.3.19a
54
Tải về để xem bản đầy đủ
43 trang
11880
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Cơ lưu chất - Chương 3: Động lực học lưu chất - Lê Văn Dực", để tải tài liệu gốc về máy hãy click vào nút Download ở trên
File đính kèm:
- bai_giang_co_luu_chat_chuong_3_dong_luc_hoc_luu_chat_le_van.pdf
