intTypePromotion=1
zunia.vn Tuyển sinh 2024 dành cho Gen-Z zunia.vn zunia.vn
ADSENSE

Bài giảng Cơ học lưu chất - Chương 3: Động lực học lưu chất

Chia sẻ: Thiên Lăng Sở | Ngày: | Loại File: PDF | Số trang:43

31
lượt xem
1
download
 
  Download Vui lòng tải xuống để xem tài liệu đầy đủ

Bài giảng Cơ học lưu chất - Chương 3: Động lực học lưu chất cung cấp cho học viên những kiến thức về hai phương pháp mô tả chuyển động của lưu chất, một số khái niệm thường dùng, phân loại chuyển đông, gia tốc của phần tử lưu chất,... Mời các bạn cùng tham khảo chi tiết nội dung bài giảng!

Chủ đề:
Lưu

Nội dung Text: Bài giảng Cơ học lưu chất - Chương 3: Động lực học lưu chất

  1. 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) r M(x,y,z) ro r r O x H.3.0a r t=0 ⇒ ro (xo , yo, zo) r r r r ∀t ⇒ r (x, y, z) → r = f ( ro , t) (3.1) hay x = x(xo , yo, zo, t) y = y(xo , yo, zo, t) z = z(xo , yo, zo, t) www.datechengvn.com Copyright @datechengvn – 2010-2015
  2. 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 2r u= & a= 2 dt dt dx dy dz ux = , uy = , uz = (3.2.a) dt dt dt d 2x d2y d 2z ax = 2 , ay = 2 , az = 2 (3.2.b) dt dt dt + Đườ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), www.datechengvn.com Copyright @datechengvn – 2010-2015
  3. 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 www.datechengvn.com Copyright @datechengvn – 2010-2015
  4. 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 u1 u2 (C1) M (C2) H.3.0b H.3.0c H.3.0d www.datechengvn.com Copyright @datechengvn – 2010-2015
  5. 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 + m 2 ; R = A/P H.3.1d H.3.1c www.datechengvn.com Copyright @datechengvn – 2010-2015
  6. 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: - Lưu lượng thể tích : H.3.1f 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 : QM = ∫∫ ρ .u.dA (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 : www.datechengvn.com Copyright @datechengvn – 2010-2015
  7. 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 1 A A ∫∫ V= = u.dA (3.5c) 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 www.datechengvn.com Copyright @datechengvn – 2010-2015
  8. 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: ρ l 2V2 ρlV ρ l V V .l const. = const → Re = = μ Vl μ μ ν 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: www.datechengvn.com Copyright @datechengvn – 2010-2015
  9. 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 ∂a u = u(x, y, z) & a = a(x, y, z) ⇔ =0 & =0 ∂t ∂t + 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’ : u x = u x + u 'x Với, u x : 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 = T ∫ u .dt t x # 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). www.datechengvn.com Copyright @datechengvn – 2010-2015
  10. 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 H.3.3b 3.3.3 Phân loại theo không gian : + Chuyển động 3 chiều : - 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. www.datechengvn.com Copyright @datechengvn – 2010-2015
  11. 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. www.datechengvn.com Copyright @datechengvn – 2010-2015
  12. 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 Δy r 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ó : ∂u x ∂u ∂u ∂u u1x = ux + Δx+ x Δy+ x Δz+ x Δt ∂x ∂y ∂z ∂t Khi đó gia tốc của phần tử lưu chất được xác định như sau : du x u − ux ax = = limΔt →0 1x dt Δt www.datechengvn.com Copyright @datechengvn – 2010-2015
  13. Trường Đại Học Bách Khoa – ĐHQG TP. HCM PGS. TS. Lê Văn Dực ∂u x Δx ∂u x Δy ∂u x Δz ∂u x = limΔt →0 ( . + . + . + ) ∂x Δt ∂y Δt ∂z Δt ∂t ∂u x dx ∂u x dy ∂u x dz ∂u x ax = + . + . + ∂x dt ∂y dt ∂z dt ∂t du x ∂u x ∂u ∂u ∂u ax = = ux+ x .uy+ x .uz+ x (3.7a) dt ∂x ∂y ∂z ∂t Tương tự : du y ∂u y ∂u y ∂u y ∂u y ay = = .ux + .uy + .uz + (3.7b) dt ∂x ∂y ∂z ∂t du z ∂u z ∂u z ∂u ∂u az = = ux + uy + z uz + z (3.7c) dt ∂x ∂y ∂z ∂t 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 : ∂u z ∂u z ∂u ux + uy + z uz ∂x ∂y ∂z 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ộ : ∂u z ∂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 r r r a= = + u .( gr ad .u ) (3.8) dt ∂t www.datechengvn.com Copyright @datechengvn – 2010-2015
  14. Trường Đại Học Bách Khoa – ĐHQG TP. HCM PGS. TS. Lê Văn Dực Với toán tử : r r ∂ r ∂ r ∂ r gr ad = ∇ = i+ j+ k ∂x ∂y ∂z + 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). www.datechengvn.com Copyright @datechengvn – 2010-2015
  15. 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: X = ∫∫∫ k .ρ .dw (3.10a) 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) Và ở thời điểm t2 ⇒ Xt2 =XBt2 + XCt2 (b) 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) (c) 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 X ( A+ B ) t2 − X ( A+ B ) t1 X Ct 2 − X At2 = + Δt Δt Δt www.datechengvn.com Copyright @datechengvn – 2010-2015
  16. Trường Đại Học Bách Khoa – ĐHQG TP. HCM PGS. TS. Lê Văn Dực dX ΔX X ( A + B ) t 2 − X ( A + B ) t1 X Ct 2 − X At 2 = lim = lim + lim (d ) dt Δt → 0 Δ t Δt→ 0 Δt Δt→ 0 Δt (1) (2) ⎛ ∂X ⎞ ∂ (1) = ⎜ ⎟ thể tích kiểm soát = ∫∫∫w k .ρ .dw (e) ⎝ ∂t ⎠ ∂t 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à : rr rr dQ = un.dA = u .n.dA → dQM = ρ.un.dA = ρ. u .n.dA rr Do đó, giá trị đại lượng X của phần tử lưu chất đi qua dA là : k.ρ. u .n.dA → rr XCt2 = Δt. ∫∫ S2 k .ρ .u .n.dA (f) Và tương tự: rr XAt2 = -Δt. ∫∫ S1 k .ρ .u .n.dA (g) r rr (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) www.datechengvn.com Copyright @datechengvn – 2010-2015
  17. Trường Đại Học Bách Khoa – ĐHQG TP. HCM PGS. TS. Lê Văn Dực Từ đó ta suy ra: rr rr X Ct2 − X At2 Δt.∫∫ k .ρ .u .n.dA +Δt.∫∫ k .ρ .u .n.dA S2 S1 (2) = lim = lim Δt →0 Δt Δt → 0 Δt rr rr = ∫∫ S1+ 2 k .ρ .u.n.dA = ∫∫ k .ρ .u .n.dA S Do đó, phương trình (d) trở thành : dX ∂X rr hệ thống = + ∫∫ k .ρ .u .n.dA (3.10b) dt ∂t CV S Với: CV : thể tích kiểm soát S : diện tích kiểm sóat. Hay : d ∂X rr dt ∫∫∫ hệ thống (k.ρ.dw) = + ∫∫ k .ρ .u .n.dA (3.10c) ∂t 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 = ∫∫∫w k.ρ.dw = ∫∫∫w ρ.dw ⇒ k=1 dm ∂m rr → = + ∫∫ ρ .u .n.dA dt HT ∂t CV S ∂ρ rr → ∫∫∫ w ∂t .dw + ∫∫ ρ .u.n.dA = 0 S (3.11) Đâ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: rr r ∫∫ ρ .u.n.dA = ∫∫∫ div( ρ .u ).dw S w www.datechengvn.com Copyright @datechengvn – 2010-2015
  18. Trường Đại Học Bách Khoa – ĐHQG TP. HCM PGS. TS. Lê Văn Dực ∂ρ r → ∫∫∫ w ∂t .dw + ∫∫∫ div( ρ.u ).dw = 0 w ∂ρ r ∫∫∫ w [ ∂t + div(ρ. u )].dw = 0 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: ∂ρ r [ + div(ρ. u )] = 0 (3.12a) ∂t 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) - Trong hệ thống tọa độ Descartes (3.13) có dạng : ∂u x ∂u y ∂u z + + =0 (3.14) ∂x ∂y ∂z - Trong hệ thống tọa độ trụ (3.13) có dạng : 1 ∂ 1 ∂uθ ∂u z (r.u r ) + + =0 (3.15) r ∂r r ∂θ ∂z 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: rr ∫∫ S ρ. u .n.dA =0 (3.11a) H.3.9 www.datechengvn.com Copyright @datechengvn – 2010-2015
  19. 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 rr rr rr rr → ∫∫S ρ. u .n.dA = ∫∫A1 ρ. u .n.dA + ∫∫ A2 ρ. u .n.dA + ∫∫ Sb ρ. u .n.dA = 0 rr Trên Sb ⇒ u .n = 0 rr Trên A1 ⇒ u .n = -u1 rr Trên A2 ⇒ u .n = u2 → - ∫∫ ρ.u1.dA + ∫∫ ρ.u2.dA = 0 A1 A2 ∫∫ A1 ρ.u1.dA = ∫∫ A2 ρ.u2.dA Qm1 = Qm2 (3.16a) ρ1.V1.A1 = ρ2.V2.A2 = const (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) Hay V1A1 =V2.A2 (3.17b) www.datechengvn.com Copyright @datechengvn – 2010-2015
  20. 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.18a H.3.18b 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 vào = ∑Qra (3.18) 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 www.datechengvn.com Copyright @datechengvn – 2010-2015
ADSENSE

CÓ THỂ BẠN MUỐN DOWNLOAD

 

Đồng bộ tài khoản
2=>2