Chương 1: Chảy rối và phương trình dòng chảy

Tài liệu cho sinh viên chuyên ngành công trình thủy

Chương 1

CHẢY RỐI VÀ PHƯƠNG TRÌNH DÒNG CHẢY

1.1 Khái niệm về dòng chảy rối:

1.1.1 Thí nghiệm về chảy rối

Nếu cho chất lỏng chảy qua một ống thủy tinh trong suốt tại đầu vào ta đặt một ngòi mực lỏng ta thấy:

Hình vẽ 1-1 Thí nghiệm về dòng chảy rối

- Khi vận tốc dòng chảy trong ống còn bé, dòng mực chảy thành đường thẳng, hay nói cách khác chất lỏng không xáo trộn lẫn nhau giữa các lớp, ta gọi trạng thái này là chảy tầng (chất lỏng được phân thành các tầng riêng biệt)

- Khi tăng vận tốc dòng chảy, dòng mực bắt đầu dao động, tiếp tục tăng vận tốc và tới một lúc nào đó dòng mực sẽ hòa tan vào chất lỏng, ta gọi trạng thái này là chảy rối (các phần tử nước chuyển động rối loạn không theo quy luật nào)

1.1.2 Phân định trạng thái chảy rối

Dòng chảy trong sông thiên nhiên nói chung là dòng chảy rối, để phân định được trạng thái chảy một cách định lượng người ta dùng số Râynôn không thứ nguyên:

υ - hệ số nhớt động học của chất lỏng; R - bán kính thuỷ lực

Với sông vì chiều rộng rất lớn so với chiều sâu nên có thể coi: TtbB

Hình vẽ 1-2 Mặt cắt ngang lòng sông

Với dòng chảy có áp:

- Khi Re<2320 ta có trạng thái chảy tầng; - Khi Re>2320 ta có trạng thái chảy rối Với dòng chảy không áp:

- Khi Re<580 ta có trạng thái chảy tầng; - Khi Re>580 ta có trạng thái chảy rối

Để nghiên cứu chảy rối người ta dùng công cụ xác suất thống kê

1.1.3 Sự hình thành dòng chảy rối:

Sự hình thành dòng chảy rối là do ma sát của chất lỏng đối với thành ống hoặc đối với bề mặt của kênh dẫn

Hình vẽ 1-3 Nguyên nhân dòng chảy rối

Bề mặt tiếp xúc với chất lỏng luôn có một nhám nhất định, phần chất lỏng chảy sát biên tiếp xúc sẽ có vận tốc thay đổi tạo thành các xoáy

- Khi vận tốc trong ống còn bé các xoáy không đủ năng lượng để lan truyền vào trong lòng chất lỏng ta có trạng thái chảy tầng;

- Khi vận tốc trong ống đủ lớn, do chất lỏng có tính nhớt (ma sát trong) nên các xoáy lan truyền vào trong lòng chất lỏng, giao thoa với nhau làm cho các phần tử nước chuyển động rối loạn, ta có trạng thái chảy rối

1.2 Đặc tính chung của dòng chảy rối:

1.2.1 Mạch động lưu tốc:

Trong chảy rối vận tốc tại một điểm đối với dòng chảy ổn định vẫn thay đổi theo thời gian:

Hình vẽ 1-4 Vận tốc tức thời trong chảy rối

Nếu gọi u là vận tốc trung bình theo thời gian tại một điểm ta có:

u - vận tốc trung bình; '

u- vận tốc theo phương x; v- vận tốc theo phương y; w- vận tốc theo phương z

Đối với dòng chảy phẳng, chảy theo một chiều theo phương x, ta có: 0

v ; w =0

Hình vẽ 1-5 Hệ trục tọa độ của dòng chảy 3 chiều

Vận tốc của dòng chảy là: u=u+u'; v=v′; w=w′

Các vận tốc mạch động có tính chất: trung bình theo thời gian của các lưu tốc mạch động sẽ bằng 0 hay: u'=0; v'=0; w'=0

1.2.2 Cường độ mạch động:

Để đánh giá độ lớn của vận tốc mạch động người ta lấy trung bình bình phương của vận tốc mạch động (sau khi bình phương mới lấy trung bình), gọi là cường độ mạch động

σ u'u max

uσ maxv'

Kªnh hë ngè

Hình vẽ 1-6 Đồ thị biến thiên của cường độ mạch động theo chiều sâu

y - Tọa dộ theo phương tẳng đứng; H - Chiều sâu dòng chảy

Ta thấy:

- Cường độ mạch động theo hướng dòng chảy phụ thuộc vào số Re;

- Cường độ mạch động theo hướng vuông góc dòng chảy không phụ thuộc vào Re

1.3 Phương trình chuyển động và phương trình liên tục của dòng chảy rối:

1.3.1 Phương trình chuyểnđộng của dòng chảy rối:

Đối với dòng chảy ổn định của chất lỏng lý tưởng ta có phương trình:

ρ

Trong đó:

1.3.2 Phương trình liên tục

Phương trình liên tục của dòng chảy rối có dạng: 0

1.3.3 Ứng suất tiếp trong chảy rối

Xét dòng chảy phẳng, đều, ổn định có phương dòng chảy trùng với trục x khi đó ta có:

Bỏ qua phương trình 2 và các đại lượng có liên quan đến trục y Do dòng chảy đều nên các đại lượng không biến đổi theo x

⇒ ' ( ' ' )=0 ∂

dòng chảy có phương trùng với trục x nên:

Vì dòng chảy ổn định nên =0.∂∂=∂∂=∂∂

Do các đại lượng chỉ phụ thuộc vào biến z nên không còn đạo hàm riêng, hệ phương trình Râynôn trở thành:

gI +υ = (nhân cả hai vế với ρ) )

Lấy tích phân theo z ta được:

Đại lượng

µ như đã biết chính là ứng suất tiếp do chảy tầng

τ = do đó đại lượng −ρu'xu'z do chảy rối ký hiệu là τr

τ = + - ứng suất tiếp toàn phần

Dựa vào công thức: τ =γIx(H −z) ta thấy τ biến thiên theo quy luật tuyến tính: Khi:

0=⇒=H τ

HIz =0⇒τ =τ0 =γ x

τ = γ.oH Jx

Hình 1-7 Sự biến thiên của ứng suất tiếp theo chiều sâu

Theo thực nghiệm biểu đồ phân bố τt và τr theo chiều sâu như sau:

- Với phần lớn dòng chảy tính từ mặt nước xuống đáy τr chiếm phần chính; - Phần sát dưới đáy τt chiếm phần lớn

τ r τ t

Hình 1-8 Sự biến thiên ứng tiếp của chảy tầng và chảy rối

Chính vì vậy khi nghiên cứu τ ở phần cách đáy ta coi τ = Khi nghiên cứu phần τrsát đáy thì τ = τt

1.3.4 Hệ số, chiều dài hỗn hợp:

Để thuận tiện cho tính toán người ta dùng các công thức tính τr có dạng của τt:

dzuExr =

- Đối với sông tự nhiên K=0,54

1.4 Dòng chảy trong sông

Xét dòng chảy trong sông do tác dụng của động học và độ dốc Chọn hệ toạ độ xyz sao cho trục x song song với mặt đáy trung bình của lòng sông có độ dốc i0 Trục z hướng lên trên như hình vẽ:

Hình 1-9 Sơ đồ dòng chảy trong sông

Lực khối tác dụng theo 3 trục như sau: 0

- Bề rộng sông lớn hơn nhiều lần chiều sâu nên gia tốc theo trục z bé hơn so với trục x và y:

có thể bỏ qua

- Cũng do bề rộng lớn hơn nhiều lần chiều sâu nên sự biến đổi của ứng suất tiếp theo z mạnh hơn nhiều so với x và y Điều này chỉ không tính đến khi dòng chảy bứt ra khỏi bờ tạo thành vũng nước xoáy trên mặt bằng

- Dòng chảy trong sông vận tốc không lớn nên áp lực phụ thêm do chảy rối nhỏ hơn so với áp lực thủy động:

- Trong phần lớn dòng chảy ứng suất tiếp do dòng chảy rối trội hơn độ nhớt nên ta có thể bỏ qua υ∆ux,υ∆uy,υ∆uz

Phương trình Râynôn rút gọn có dạng như sau:

Phương trình liên tục: 0=∂∂+∂∂+∂∂

Tích phân phương trình 3 trong hệ phương trình chuyển động ta có:

p =−γ ⇒=−γ do p chỉ phụ thuộc z

p =γ =γ∂

∂

Hệ phương trình chuyển động trở thành:

Phương trình liên tục: 0=∂∂+∂∂+∂∂

1.5 Sự phân bố lưu tốc trong kênh hở:

Từ khái niệm chảy rối có thể tìm ra các công thức phân bố lưu tốc Sự phân bố lưu tốc và sức cản trong sồn thiên nhiên chiếm 1 vị trí quan trọng trong thuỷ động lực học

1.5.1 Công thức phân bố lưu tốc bán thực nghiệm:

1.5.1.1 Sự phân bố lưu tốc trong vách nhẵn:

Từ khái niệm ứng suất tiếp chảy rối ta có công thức tính ứng suất tiếp trong chảy rối:

Gần đáy ứng suất tiếp bằng ứng suất tiếp tại đáy hay τ =τ0 =γHIthì ta có:

Tích phân hai vế ta có: zCK

Với C là hằng số tích phân Nếu đặt: 0= HI = gHI =u*

gọi là lưu tốc ma sát (hay lưu tốc động lực) thì:

Gọi umaxlà vận tốc lớn nhất trên mặt nước (z=H) ta có:

CHuK

⇒

u - Lưu tốc bình quân của thuỷ trực

Vị trí điểm có lưu tốc bằng lưu tốc trung bình của thuỷ trực u =ubq: Ta có:

Khi u =ubq thì: ln=1

hay ez

Công thức của Gôntrarốp:

;

;

Trong đó:

1.5.2 Các công thức phân bố lưu tốc kinh nghiệm:

Các công thức thực nghiệm được xây dựng bằng cách xấp xỉ các số liệu thực đo theo các cách khác nhau thùy thuộc từng tác giả

1.5.2.1 Công thức dạng đường Parabol - Công thức Bazin:

Trong đó:

m - Hệ số có giá trị 20 ÷ 24; H – Chiều sâu dòng chảy; I - Độ dốc mặt nước

1.5.2.2 Công thức dạng log – Công thức Nicuratdê:

= * 1 1 lnξ

Hoặc:

= 1 1 lnξ

Trong đó:

u* = - Lưu tốc ma sát đáy; C - Hệ số Sedi;

ξ - độ sâu tương đối

1.5.2.3 Công thức dạng đường Ellipse - Karausep:

Căn cứ vào tài liệu thực đo của rất nhiều sông ở Liên Xô, Karausep tìm ra công thức sau:

⎝⎛ −−=

Trong đó:

uMP= bq

M phụ thuộc vào C:

- Khi 10 ≤C ≤ 60 thì M=0,7C+6;

- Khi 60 ≤C ≤960 thì M=48

Chương 1 1-1

1.1 Khái niệm về dòng chảy rối: 1-11.2 Đặc tính chung của dòng chảy rối: 1-21.3 Phương trình chuyển động và phương trình liên tục của dòng chảy rối: 1-41.4 Dòng chảy trong sông .1-71.5 Sự phân bố lưu tốc trong kênh hở: 1-9