Chúng ta đã biết là trong môi trường chất lỏng tĩnh ( không có sự chuyển động tương đối giữa các lớp chất lỏng). Tác dụng lên chất lỏng chỉ có thành phần lực khối.
Để xây dựng phương trình Euler tĩnh, chúng ta chọn hệ trục tọa độ gắn chặt với bình bể chất lỏng. Trong môi trường chất lỏng tĩnh lấy một điểm M bất kỳ có tọa độ x, y ,z và có áp suất p (Hình 1).
Tách ra một khối chất lỏng có dạng hình hộp chứ nhật, có các cạnh song song với các trục tọa độ, có độ dài tương ứng là dx, dy, dz. Điểm M là một đỉnh của hình hộp chữ nhật. Xét điều kiện cân bằng của khối chất lỏng được tách ra. Giả sử khối chất lỏng chịu sự tác động lực khối tổng hợp từ 3 thành phần theo phương các trục Ox, Oy, Oz. Lực khối đơn vị theo các trục tương ứng là X, Y, Z. Khi đó lực khối thành phần theo từng phương trục tọa độ bằng tích khối lượng của khối chất lỏng nhân với lực khối đơn vị tương ứng.
Áp suất trong lòng khối chất lỏng là một hàm không gian, tức là p=p(x,y,z). Theo tính chất của áp suất thủy tĩnh tại điểm M giá trị áp suất bằng nhau theo mọi phương, giả sử tại M có áp suất là pM. Khi di chuyển điểm M tới điểm N. Theo phương Ox tọa độ điểm M thay đổi một lượng dx. Kéo theo đó hàm p(x,y,z) biến thiên một lượng bằng vi phần từng phần (∂p/∂x)dx, vì vậy áp suất tại điểm N là pM+(∂p/∂x)dx.
Ở đó (∂p/∂x) – được gọi là gradient áp suất điểm M theo phương Ox.
Xét các cặp điểm tương ứng trên 2 mặt giới hạn vuông góc với Ox của khối chất lỏng, ví dụ M' và N', ta đều thấy rằng giữa các cặp điểm này khác nhau một lượng áp suất.
Tác dụng lên khối chất lỏng theo phương Ox có 2 thành phần lực: lực khối và lực chênh áp giữa 2 mặt giới hạn của khối hình hộp.
Lực khối theo phương Ox là:
Fkx=Xρdxdydz
(1.1)
Hệ 3 phương trình (1.1) là phương trình Euler tĩnh ở dạng vi phần từng phần.
Cộng theo vế 3 phương trình trên ta được:
Ta thấy tích phân toàn phần áp suất:
Suy ra ta có:
dp=ρ(Xdx+Ydy+Zdz) (1.2)
Phương trình (1.2) là phương trình Euler tĩnh ở dạng toàn phần.
Xét một số trường hợp đơn giản.
Trường hợp 1 :
Môi trường chất lỏng đứng yên chỉ chịu tác dụng của trọng lực.
Khi đó lực khối đơn vị Z= –g ; Y=X=0.
Ta có phương trình vi phân cân bằng :
dp= –ρgdz
Lấy tích phân ta có
p= –ρgz+C
C – hằng số tích phân
Trường hợp 2:
Chuyển động quay tròn đều quanh trục đứng của bình chứa chất lỏng.
Bình chứa chất lỏng mở quay đều với vận tốc góc ω quanh trục đứng của bình. Quan sát hiện tượng thấy chất lỏng cũng quay theo bình với cùng vận tốc ω; mặt thoáng của chất lỏng trong bình ở trung tâm hạ thấp xuống và ở mép thì nâng cao lên. Mặt thoáng của chất lỏng trong bình trở thành mặt dạng tròn xoay. Tác dụng lên khối chất lỏng trong trường hợp này có hai lực khối – trọng lực và lực ly tâm, tương ứng với 2 lực khối đơn vị là g và ω2r.
Áp dụng phương trình (1.2) để xác định qui luật thay đổi áp suất trong trường hợp này. Ta có:
Z= –g; X= ω2rcos(r,x) = ω2x; Y= ω2rcos(r,y)=ω2y; với x2+y2=r2
Thế vào phương trình (1.2) ta có :
dp=ρ(ω2xdx+ ω2ydy – gdz)
Lấy tích phân 2 vế thu được :
p= ρ ω2(x2+y2)/2 – ρgz +C
p= ρ ω2r2/2– ρgz +C
C – hằng số tích phân
Khi r=0, z=h, p=pat ta có : pat= – ρgh+C
Suy ra C=pat+ ρgh
Vậy p= ρ ω2r2/2 – ρgz+ pat+ ρgh
p= pat+ ρ ω2r2/2 + ρg(h – z)
Để xác định phương trình mặt thoảng của chất lỏng trong bình ta có
p= pat
Suy ra : z=h+ ω2r2/(2g)
Bài viết hôm nay của mình tạm dừng tại đây.
